Electromagnetic excitation of the Z\( 1232) -resonance 



Vladimir Pascalutsa Marc Vanderhaeghen Shin Nan Yang ^ 

^Physics Department, The College of William & Mary, Williamsburg, VA 23187, USA 

^Theory Center, Jefferson Lab, 12000 Jefferson Ave, Newport News, VA 23606, USA 

^ECT*, Villa Tambosi, Str delle Tabarelle 286, 1-38050 Villazzano (Trento), Italy 

'^Department of Physics and Center for Theoretical Sciences, National Taiwan University, 

Taipei 10617, Taiwan 



Abstract 

We review the description of the lowest-energy nucleon excitation — the A(1232)-resonance. 
Much of the recent experimental effort has been focused on the precision measurements of 
the nucleon-to-A transition by means of electromagnetic probes. We confront the results of 
these measurements with the state-of-the-art calculations based on chiral effective-field the- 
ories (EFT), lattice QCD, large- A'^c relations, perturbative QCD, and QCD-inspired models. 
We also discuss the link of the nucleon-to-A form factors to generalized parton distributions 
(GPDs). Some of the theoretical approaches are reviewed in detail, in particular, recent dy- 
namical and unitary-isobar models of pion electroproduction, which are extensively used 
in the interpretation of experiments. Also, the novel extension of chiral EFTs to the energy 
domain of the A-resonance is described in detail. The two-photon exchange effects in the 
electroexcitation of the A-resonance are addressed here as well. 
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1 Introduction 



The A(1232) resonance — the first excited state of the nucleon — dominates the pion- 
production phenomena and plays a prominent role in the physics of the strong interaction. This 
resonance, first witnessed more than 50 years ago by Fermi and collaborators [1] in pion scat- 
tering off protons at the Chicago cyclotron (now Fermilab), is a distinguished member of the 
large family of excited baryons. It is the lightest baryon resonance, with a mass of 1232 MeV, 
less than 300 MeV heavier than the nucleon. Despite its relatively broad width of 120 MeV 
(corresponding with the lifetime of 10~^^s), the A is very well isolated from other resonances. 
It is almost an ideally elastic vriV resonance, 99% of the time it decays into the nucleon and 
pion, A Nn. The only other decay channel - electromagnetic, A N-j, contributes less 
than 1% to the total decay width (even though this channel will preoccupy the larger part of this 
review). The A(1232) has isospin 3/2 and as such comes in four different charge states: A++, 
A+, A°, and A~, all with (approximately) the same mass and width. The spin of the A(1232) 
is also 3/2, and it is the lightest known particle with such a spin. 

The A-resonance dominates many nuclear phenomena at energies above the pion-production 
threshold [2] (for recent reviews see [3,4]). In cosmology it is largely responsible for the "GZK 
cut-off" [5], which occurs due to the suppression of the high-energy cosmic ray flux by the cos- 
mic microwave background (CMB). Once the energy of the cosmic rays is sufficient to produce 
the A-resonance in the scattering off the CMB photons, the rate of observed rays drops dra- 
matically (for a recent calculation see [6]). This effect puts a cutoff on the primary cosmic ray 
energy at around 10^*^ eV for the rays coming from a distance larger than a few tens of Mpc. 

In a laboratory, the A's are produced in scattering the pion, photon, or electron beams off a 
nucleon target. High-precision measurements of the — > A transition by means of electromag- 
netic probes became possible with the advent of the new generation of electron beam facilities, 
such as LEGS, BATES, ELS A, MAMI, and Jefferson Lab. Many such experimental programs 
devoted to the study of electromagnetic properties of the A have been completed in the past few 
years. 

The current experimental effort has been accompanied by exciting developments on the the- 
oretical side, most recently in the fields of lattice QCD and chiral effective-field theories. These 
recent experimental and theoretical advances in understanding the A-resonance properties are 
the main subjects of this review. 

The electromagnetic — > A (or, in short 'jNA) transition is predominantly of the mag- 
netic dipole (Ml) type. A first understanding of the 7A^A transition can be obtained based on 
symmetries of Quantum Chromodynamics (QCD) and its large number-of-color (Nc) limit. The 
spin-flavor global symmetry of QCD, is utilized by many quark models and is exactly realized 
in the large- A/'c limit. In the quark-model picture, the ^ A transition is described by a spin 
flip of a quark in the s-wave state, which in the 7A^A case leads to the magnetic dipole (Ml) 
type of transition. Any d-wave admixture in the nucleon or the A wave functions allows also for 
the electric (£'2) and Coulomb (C2) quadrupole transitions. Therefore, by measuring the latter 
two transitions, one is able to assess the presence of the o?-wave components and hence quantify 
to which extent the nucleon or the A wave function deviates from the spherical shape. In this 
way one can hope to understand to which extent these particles are "deformed". 

The d-wave component of A's wave function can be separately assessed by measuring the 
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electric quadrupole moment of the A. However, this would be extremely difficult because of the 
tiny lifetime of the A. The small d-state probability of the A's wave function also enters into 
the A magnetic dipole moment, which is being measured, but the extraction of a small number 
from such a quantity is also very complicated. The A transition, on the other hand, was accu- 
rately measured in the pion photo- and electro-production reactions in the A-resonance energy 
region. The E2 and C2 transitions were found to be relatively small but non-zero at moderate 
momentum-transfers (Q^), the ratios Rem — E2/M1 and Rsm — C2/M1 are at the level of a 
few percent. 

Because the A excitation energy is only around 300 MeV and because the A-resonance al- 
most entirely decays into t^N, pions are expected to play a prominent role in the A properties. 
Early calculations of the 7A^A transition within chiral bag models revealed the importance of 
the pion cloud effects. Recall that pions are the Goldstone bosons of the spontaneously broken 
chiral symmetry of QCD. As the strength of the Goldstone boson interactions is proportional 
to their energy, at sufficiently low energy a perturbative expansion is possible. With the advent 
of the chiral effective field theory (xEFT) of QCD [7-9] and its extensions to the A-resonance 
region, it has become possible to study the nucleon and A-resonance properties in a profoundly 
different way. The advantages of such an approach are apparent: xEFT is a low-energy effective 
field theory of QCD and as such it provides a firm theoretical foundation, with all the relevant 
symmetries and scales of QCD built in consistently. The 7A'^A transition provides new chal- 
lenges for xEFT as it involves the interplay of two light mass scales : the pion mass and the 

— A mass difference. This topic will preoccupy a large part of this review, and we discuss up 
to which energy/momentum scales such an approach can be expected to hold. 

Considerable progress has recently been achieved as well in the lattice QCD simulations of 
hadronic properties. They hold the promise to compute non-perturbative properties of QCD 
from first principles. The present state-of-the-art results for hadronic structure quantities, such 
as the 7A^A transition form factors, are obtained for pion masses above 300 MeV. Therefore, 
they can only be confronted with experiment after an extrapolation down to the physical pion 
mass of 140 MeV. Such extrapolation can be obtained with the aid of xEFT, where the pion 
mass dependence is systematically calculable. The xEFT framework thus provides a connec- 
tion between lattice QCD calculations and experiment. 

The quark structure of the iV ^ A transition is accessible through the phenomenon of asymp- 
totic freedom of QCD at short distances. In a hard scattering process, such as 7* 7A where 
the virtual photon 7* transfers a large momentum, the QCD factorization theorems allow one 
to separate the perturbative and non-perturbative stages of the interaction. In this way, one ac- 
cesses in experiment the non-perturbative matrix elements parametrized in terms of new parton 
distributions, which are generalizations of the quark distributions from deep inelastic scattering 
experiments. One obtains in this way quark distribution information for the A transition, 
and the 7A^A form factors are obtained as first moments of such N ^ ^ generalized parton 
distributions. 

Traditionally, the resonance parameters are extracted from pion electroproduction experi- 
ments by using unitary isobar models, which in essence are unitarized tree-level calculations 
based on phenomenological Lagrangians. However, as discussed above, at low Q"^ the 7A^A- 
transition shows great sensitivity to the "pion cloud", which until recently could only be com- 
prehensively studied within dynamical models, which will also be reviewed here. 

The outline of this review is as follows. In Sect. 2, we introduce the definitions for the 7A'^A 
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and 7A A form factors and review the experimental status of tlie 7iVA transition at the real pho- 
ton point. We then discuss our present theoretical understanding of the 7A^A transition. Some 
physical insight is obtained from QCD inspired models based on quark degrees of freedom, 
pion degrees of freedom, as well as from the large Nc limit. The chiral symmetry of QCD is 
reviewed and the role of pions in the 7A'^A transition highlighted. We review the status of the 
lattice QCD calculations and their limitations. Subsequently, we discuss the quark structure of 
the ^ A transition. We introduce generalized parton distributions for the ^ A transition 
and discuss our phenomenological information of such distributions based on their relation to 
the 7A'^A form factors. Finally, we confront predictions made by perturbative QCD at very large 
momentum transfers with the available data. 

The theoretical formalism of dynamical models will be discussed in detail in Sect. 3, and 
compared with other approaches such as unitary isobar models and dispersion relations. 

In Sect. 4, we will review how xEFT can be applied to study the 7A/^A-transition at low mo- 
mentum transfers. It will be discussed that the xEFT may provide a theoretically consistent and 
phenomenologically viable framework for the extraction of the resonance parameters. In par- 
ticular we will review how the A can be introduced in a xEFT and discuss the field-theoretic 
intricacies due to the spin-3/2 nature of the A. We will apply the xEFT formalism to pion elec- 
troproduction in the A region and discuss applications to other complementary processes such 
as Compton scattering and radiative pion photoproduction in the A region. The convergence of 
the perturbative xEFT expansion will also be addressed. 

In Sect. 5, we compare predictions of both dynamical models and xEFT to the pion photo- 
and electroproduction observables. 

The information on the 7A^A transition as discussed above is obtained from pion electro- 
production assuming that the interaction is mediated by a single photon exchange. In Sect. 6, 
we give a brief description of our understanding of corrections to this process beyond the one- 
photon exchange, which may become relevant with increasing Q^. 

Finally in Sect. 7, we give our conclusions and spell out some open issues in this field. 

Several reviews exist in the literature on various aspects of the 7A^A transition. Among the 
more recent reviews, we refer the reader to Ref. [10] where the formalism to extract 7iVA 
form factors from pion photo- and electroproduction observables is outlined in detail. The ex- 
perimental status of the 'jNA, and 'jNN* transitions as obtained from meson photoproduction 
experiments is well described in Ref. [11]. The status of pion electroproduction experiments in 
the A region as well as a previous review of dynamical models can be found in Ref. [12]. 
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2 The electromagnetic — > A transition in QCD 



Both the nucleon and the A(1232) can be considered as quantum states in the rich and com- 
plex spectrum of the quark-gluon system. It is well established that the underlying theory which 
should describe this spectrum is Quantum Chromodynamics (QCD). However, direct calcula- 
tions in QCD of quantities such as the baryon spectrum, form factors, and parton distributions 
are extremely difficult, because they require non-perturbative methods. Only the lattice simula- 
tions have achieved some limited success to date in computing the hadron properties from first 
principles. Even then, such calculations are severely limited by the presently available comput- 
ing power. 

Fortunately, some insight into the properties of the spectrum can be obtained based on gen- 
eral principles of QCD, such as the global symmetries and large number-of -colors (Nc) limit 
of QCD. After introducing the definitions of 7*A^A and 7* A A form factors in Sect. 2.1, and 
reviewing the experimental status on the 7A^A transition at the real photon point in Sect. 2.2, 
we shall discuss in Sect. 2.3 the predictions for the 7iVA transition due to the spin-flavor sym- 
metry, utilized by many quark models and exactly realized in the large- A^'c limit. Furthermore, 
we shall compare predictions for the 7A'^A transition from varios "QCD-inspired models" based 
on quark and/or pion degrees of freedom. 

in Sect. 2.4 we shall review some large-iVc relations relevant to the 7 A^A transition. 

In Sect. 2.5, we shall discuss the chiral symmetry of QCD and the role of pions in the elec- 
tromagnetic N ^ A transition. The chiral effective field theory is able to provide predictions 
for the low momentum transfer behavior of the 7*A^A form factors. These predictions will be 
confronted with the dynamical models and experiment. 

In Sect. 2.6, we shall discuss the lattice QCD calculations of the 7A^A form factors. Present 
lattice QCD calculations are performed for quark masses sizably larger than their values in na- 
ture, corresponding with pion mass values of around 0.3 GeV or larger. It will be discussed how 
the chiral effective field theory can be useful in extrapolating the present lattice QCD calcula- 
tions to the physical pion mass. 

In Sect. 2.7, we shall review the quark structure of the — A transition and discuss gener- 
alized parton distributions (GPDs) for the electromagnetic ^ A transition. The GPDs can 
be accessed in hard exclusive processes such as deeply virtual Compton scattering where the 
hard probe ensures that the process occurs on a quark which is taken out of the initial nucleon 
and inserted into the final A. We shall see that the 7*A^A form factors are obtained as the first 
moment in the struck quark momentum fraction of such GPDs, which provide much richer in- 
formation on the A^ ^ A transition at the quark level. 

Finally, in Sect. 2.8, we shall consider the A^ — A transition in the formalism of perturbative 
QCD (pQCD). These considerations are valid only at very small distances, where quarks nearly 
do not interact — asymptotic freedom. In this limit, the 7*A'^A form factors correspond with a 
hard photon which hits a quark in the nucleon. The struck quark shares the large momentum 
with the other two (near coUinear) valence quarks in such a way that the final three quark state 
has A quantum numbers. We shall discuss the predictions made in this limit and confront them 
with the experimental status of the 7*A^A form factors at large momentum transfers. 
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Fig. 1. 'y*NA vertex (left panel) and vrA^A vertex (right panel). The four-momenta of the nucleon (A) 
and of the photon (pion) are given by p (p') and q (k) respectively. The four-vector index of the spin 3/2 
field is given by a, and fi is the four-vector index of the photon field. 

2. 1 Definitions and conventions 



Throughout this review we shall use the following conventions for the metric and 7-matrices: 
g'^-' = diag(l, -1,-1,-1), {7^7"} = 2^?'^^ 75 = WtV, ^0123 = +1- The ^NA and 
7AA form factors are introduced as follows. 

2.1.1 The 7*A^A vertex and form factors 

The general spin structure of the pion- and photon-induced — > A transitions, see Fig. 1, 
can be written as: 

(A(p') I e J^(0) I N{p)) ^ zu^ip') V'/^^uip) 

^ . [2 3e(MA + Mjv) 

" ' V 3 2M^[(Ma + Miv)2 + Q2] 

+ 9e{Q') {q" p'' - q ■ p' g'^^) 

+ gem (g" q' - q' g"-') H,} u{p), (2.2) 

where k is the pion and q is the photon 4-momentum, p (p') is the nucleon (the A) 4-momentum, 
Mtv (Ma) is the nucleon (the A) mass, u is the nucleon spinor, and represents the spin-3/2 A 
vector-spinor In Eqs. (2.1) and (2.2), the spin dependence in the N and A spinors is understood. 
The operator in Eq. (2.1), with a = 1, 2, 3 corresponding with the Cartesian pion fields vr", is 
the isospin 1/2 3/2 transition operator. The operator J'^ is the electromagnetic current opera- 
tor, and the factor y^2/3 in front of Eq. (2.2) corresponds with the isospin factor for the photon 
induced p A+ transition. Furthermore, ~ 92.4 MeV is the pion decay constant, Ha is a 
dimensionless constant representing the strength of the vriV — > A transition and thus related to 
the decay width of the A to nN, see Sect. 4.5.1. 

The strong transition form factor F^tata is normalized as F-^N/\{ml.) = 1, and can in principle 
depend as well on the invariant mass of the nucleon and A, i.e.,p^ andp'^. However, unless ex- 
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plicitly specified, we assume that nucleon and A are on the mass shell and therefore = Mjj, 
^u{p) = Mnu{p),'p'^ = Ml,fua{p') = MAU^ip'), andp>"(p') = = 7aM"(p')- 
We use the covariant normalization convention for the spinors, i.e. u{p)u{p) = 2Mjv, and 

In Eq. (2.2), the electromagnetic form factors gM, Qe, and gc represent the strength of the 
magnetic dipole, electric quadrupole, and Coulomb quadrupole ^ A transitions, respec- 
tively, as a function of the momentum transfer: Q'^ = — g^. Note that for the real photon case, 
Q'^ = 0, only the magnetic and electric transitions are possible, the Coulomb term drops out. 

These electromagnetic form factors relate to the more conventional magnetic dipole (G^), 
electric quadrupole (G*^) and Coulomb quadrupole (Gq) form factors of Jones and Scadron [13] 
as follows, ^ 



G*m{Q'^ 



■9m 



Q 



+ 



\{-Mi + Mf, + Q')gE + Q'gc 



Ge{Q^) = [h{-Ml + + gE + Q^gc 

GhiQ^) = 7^ + + _ 2Ml gE 

V+ 

where Q± is defined as: 



(2.3) 



(2.4) 



Equivalently, one can also parametrize the 7* A^A transition through three helicity amplitudes 

^1/2. ^3/2 and 5*1/2, which are defined in the A rest frame as illustrated in Fig. 2. These A rest 
frame helicity amplitudes are defined through the following matrix elements of the electromag- 
netic current operator: 



^3/2 



1 



e 1 



1/2 



S 



e 1 



1/2 



V2gZ (4MivMA)V2 



( A(0, +3/2) I J • eA=+i I N{-q, +1/2) ), 
( A(0, +1/2) I J • eA=+i I N{-q, -1/2) ), 
A(0, +l/2)|jO|iV(-g, +1/2)), 



(2.5) 



where the spin projections are along the 2; -axis (chosen along the virtual photon direction) and 

where the transverse photon polarization vector entering in A1/2 and is given by eA=+i = 
— l/\/2(l, i, 0). Furthermore in Eq. (2.6), e is the proton electric charge, related to the fine- 
structure constant as aem = e^/(47r) ~ 1/137, and ^a is the magnitude of the virtual photon 
three-momentum in the A rest frame: 



^ The form factors G\j, G*^, and used throughout this work correspond with the ones also used in 
Ref. [14]. Our sign convention leads to positive values of all three form factors at = 0. 
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+1^ -in 
7* ^AA/w\/^» • TV ^1/2 

A 



H/2 



+1 



-1/2 



7* ^AA/^AAA# ' TV -43/2 

A 



f3/2 







hl/2 
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7* ^AA/^AAAJ• ^ TV ^1/ 

A 

+1/2 

Fig. 2. 7*TV — > A helicity amplitudes in the A rest frame. The 7, TV, and A spin projections onto the 
z-axis, which is chosen along the virtual photon direction, are indicated on the figure. The corresponding 
helicity amplitudes are indicated on the right. 

= |q| = (2.6) 

The helicity amplitudes are functions of the photon virtuality Q^, and can be expressed in terms 
of the Jones-Scadron 7*TVA form factors, by using Eq. (2.2), as : 

^i/2 = -iV^{Q/-3G^,}, 

Sin-^^^O^c. (2.7) 
where TV is defined as : 

The above helicity amplitudes are expressed in units GeV~^/^, and reduce at = to the 
photo-couplings quoted by the Particle Data Group [15]. 

Experimentally, the 7*TVA helicity amplitudes are extracted from the Ml, E2, and C2 mul- 
tipoles for the 7*TV ttTV process at the resonance position, i.e. for ttTV cm. energy W = M^. 

(3/2) (3/2) (3/2) 

These pion electroproduction multipoles are denoted by Mij^ \ Ei_l and 5)+ (following 
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standard notation, see e.g. Ref. [16] and Sect. 3 for more details), where the subscript refers 
to the partial wave / = 1 in the nN system and "+" indicates that the total angular momen- 
tum is J = 1 + 1/2, being 3/2. The superscript (3/2) in the multipole notation refers to the 
total isospin 3/2. The 7*A'^A helicity amplitudes ^3/2, ^1/2 > and S1/2 are obtained from the 
7*A^ — > nN resonant multipoles as [17]: 



2 

Ay2 = -l {Mir+SE[f^}, 



S,/2 = -y/2S[f\ (2.9) 

where Mf!/^\ eI^/^\ and S[^_l^^ are obtained from the imaginary parts of the resonant multi- 
poles at the resonance position W — Ma as: 



Mlf\Q') = J^a^lmMif\Q',W = Ma), (2.10) 



where ^2/3 is an isospin factor, and similar relations define E^^^^-* and S[^_l'^'' . In these relations, 
a A is given by: 

/47rfcAMArAy/' 

OA = 77 , (2.11) 

where Fa denotes the A width, which we take as Fa = 0.115 GeV, and /ca denotes the magni- 
tude of the pion three-momentum in the nN cm. frame at the resonance position (W — Ma): 



|k| = 



2Ma 



1/2 



(MA + Mjv)'-m^ • (MA-M^)^-m^ 



1/2 



(2.12) 



which yields = 0.227 GeV. Furthermore, for = 0, ^a = 0.259 GeV and oa ^ 

1.29 GeVi/2. 

We can extract the Jones-Scadron 7*A^A form factors from the multipoles at the resonance 
position. Combining Eqs. (2.7) and (2.9), we obtain: 



E[f^ = -NG*E, 

-^ir^-A^^G^- (2.13) 

Of special interest are the multipole ratios at the position of the A resonance. The ratio of elec- 
tric quadrupole (E2) over magnetic dipole (Ml) is denoted by : Rem — E2/M1 (sometimes 
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also denoted by EMR), whereas the ratio of Coulomb quadrupole (C2) over magnetic dipole 
is denoted by Rsm — C2/M1 (sometimes also denoted by CMR): 



_ ^(f ) ^ _ .1^ ^^^^ 
Rem = EMR = _ ,„^ — — y=— — , 

5(3/2) /^c. 

Rsm ^ CMR ^ %^ = , " (2.14) 

Using the relations of Eq. (2.13), these ratios can be expressed in terms of the Jones-Scadron 
form factors as: 



From the values of the 7*A'^A form factors at Q"^ — 0, one can extract some interesting 
static quantities. For the dominant Ml transition, one can extract the static N ^ transition 
magnetic moment /jLn-^a from the value of at = as [18]: 



I^N^A = ^^GI,{0), (2.16) 

which is expressed in nuclear magnetons /jLn = e/(2Mjv). 

Furthermore, from the value of G*^ at Q'^ — 0, one can extract a static — > A quadrupole 
transition moment Qn-*a as [18]: 



where qA{0) is obtained from Eq. (2.6) for = 0, as qA{0) = {Ml - M^)/2Ma. 

Finally, we also like to note that in the literature, one sometimes uses so-called Ash form 
factors for the 7*A^A transition [19]. They are simply related to the Jones-Scadron form factors 
as: 



Gl,...M') = ^^%i^G^(Q^), (2.18) 
and analogously for the electric and Coulomb quadrupole form factors. 



2.1.2 The 7* A A vertex and form factors 

Consider the coupling of a photon to a A, Fig. 3. The matrix element of the electromagnetic 
current operator between spin 3/2 states can be decomposed into four multipole transitions: a 
Coulomb monopole (CO), a magnetic dipole (Ml), an electric quadrupole (E2) and a magnetic 
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Fig. 3. The 7* A A vertex. The four-momenta of the initial (final) A and of the photon are given by p (p') 
and q respectively. The four-vector indices of the initial (final) spin 3/2 fields are given by j3 (q), and // 
is the four-vector index of the photon field. 

octupole (M3). We firstly write a Lorentz-covariant decomposition which exhibits manifest 
electromagnetic gauge-invariance: 



+ 
+ 



2Ma 



i J y 



(2Ma) 



(2.19) 



where F* are the 7*AA form factors, ca is the electric charge in units of e (e.g., eA+ = +1), 
such that F*{0) = 1 . The relation to the multipole decomposition [20,2 1 ] can be written in terms 
of the magnetic dipole (/xa), electric quadrupole (Qa) and magnetic octupole (Oa) moments, 
given by: 



Qa 
Oa 



2Ma 

e 

Ml 



[eA + F;m 
CA - iF3*(0) 



2Mi 



eA + F;(0)-H^3(0)+i^:(0)) 



(2.20a) 
(2.20b) 
(2.20c) 



2.2 Experimental information at the real photon point 



In this section we briefly summarize the experimental information on the 7A^A transition at 

The 7A^A transition can be studied both in nucleon Compton scattering 7A^ — > 7A^ and in 
the pion photoproduction reaction 7A^ nN. The highest precision data on the 7A^A Ml 
and E2 amplitudes have been obtained in pion photoproduction experiments on a proton target, 
using linearly polarized photons, both by the MAMI/A2 Collaboration [22,23] and the LEGS 
Collaboration [24,25]. A detailed discussion of these 7A^ nN data will be presented in 
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Sect. 5.1. 

We quote here the real photon 7A^A amplitudes at the resonance position (W = Ma) for a 
proton target resulting from these high precision measurements. The MAMI/A2 Collaboration 
obtained the values [23]: 



Ai/2 = -{l31±l) [IQ-'^GeV-^/^], 
A3/2 = -(251±l) [lO-^GeV^^/^], 

Rem^- (2.5 ± Q.lstat. ± 0.2,^,0 %. (2.21) 
The LEGS Collaboration obtained the values [25]: 

= - (135.7 ± 1.3stat.+syst. ± 3.7weO [IQ-'GeV-'/^ 
A3/2 = - (266.9 ± 1.6stat.+syst. ± 7.8wd) [10-'GeV-'/% 

Rem = - (3.07 ± 0.26,tat.+syst. ± 0.24^odei) %■ (2.22) 



As summary of the above experiments, the Particle Data Group quotes as values [15] 



2 



2 = -(135 ±6) [lO-^GeV-^/^]^ 
2 = -(250 ±8) [lO-^GeV-^/^], 



Al/2 = - 

A3/2 = - , , , 

i?EM=-(2.5±0.5) %. (2.23) 



The experimental information on the Ml ^yNA transition can equivalently be expressed in 
terms of Gm(0) using Eq. (2.13) or equivalently in terms of the transition magnetic moment 
A*p^A+ using Eq. (2.16). Using the experimental values of Eq. (2.21), Ref. [26] extracted the 
values: 



Gm{0) = 3.02 ±0.03, 

A*p-.A+ = [3.46 ± 0.03] UN- (2.24) 

The value for the E2 ^NA transition can equivalently be expressed in terms of a quadrupole 
transition moment using Eq. (2.17). Using the experimental values of Eq. (2.21), Ref. [18] 
extracted: 



Qp_,^+ = - (0.0846 ± 0.0033) fm^ (2.25) 



^ Note that the central values of yli/2 and A2/2 quoted by PDG[15] yield Rem = —1.64 % when using 
Eq. (2.14), and therefore do not exactly correspond with the central value Rem = —2.5 % they quote. 
The photo-couplings are based on an analysis of a larger number of experiments, and within their error 
bars are compatible with the quoted value for Rem though. 
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2.3 Physical interpretation: quark and pion-cloud models 

2.3.1 Constituent quark models 

In a quark model, the nucleon appears as the ground state of a quantum-mechanical three- 
quark system in a confining potential. In such a picture, the ground state baryons (composed 
of the light up (u), down (d) and strange (s) quark flavors) are described by SU (6) spin-flavor 
symmetric wave functions, supplemented by an antisymmetric color wave function. 

In the Isgur-Karl model [27], the constituent quarks move in a harmonic oscillator type con- 
fining potential. For the ground state baryons, the three constituent quarks are in the Is oscillator 
ground state, corresponding with the [56]-plet of SU (6). The harmonic oscillator states can be 
represented by \B '^^'^^Lj)t, where B stands for either or A states, S specifies the spin, L the 
orbital angular momentum (L — S, P, D, ... in the common spectroscopic notation), and J the 
total angular momentum of the three-quark state. Furthermore, t{= S, M, A) refers to the sym- 
metry type : symmetric (S), mixed symmetric (M) or anti- symmetric (A) under exchange of 
the quarks in both the spin-flavor and space parts of the baryon wave function. In the Isgur-Karl 
model, the long-range confining potential is supplemented by an interquark force correspond- 
ing with one-gluon exchange. The one-gluon exchange leads to a color hyperfine interaction 
between quarks i and j of the form: 



(2.26) 

with as the strong coupling constant, Si(mj) the spin (mass) of quark i, and where rij(rij) 
specify the vector (distance) between quarks i and j. The first term in Eq. (2.26) corresponds 
with a zero-range spin-spin interaction, whereas the second term corresponds with a tensor 
force. The color hyperfine interaction of Eq. (2.26) breaks the SU{6) symmetry and leads to 
a mass splitting between iV(939) and A(1232), often referred to as the hyperfine splitting. It 
was found that it also predicts well the mass splittings between octet and decuplet baryons [28]. 
Furthermore, the tensor force in Eq. (2.26) will produce a £)-state (L — 2) admixture in the N 
and A ground states [29,30]. Because of this hyperfine interaction, the A^(939) and A(1232) 
states are described as superpositions of SU{Q) configurations. Including configurations up to 
the 2huj oscillator shell, they are given by (using the abovementioned spectroscopic notation): 

|7V(939)) = as \N ^8^/2)3 + a's \N -'S[/^)s + '^i/2)m + an \N ^Dy^)^.^]) 

I A(1232)) = hs I A %2)5 + h's | A ^,5^/2)5 + 6d | A ^D^,2)s + &d I A 'L>3/2)m. (2.28) 

By diagonalizing the hyperfine interaction and fitting the results to the baryon spectrum, Isgur 
et al. [30] obtained the following values for the wave-function coefficients: 

as ~ 0.93, a!s ^ -0.29, om ^ -0.23, an ^ -0.04, (2.29) 
hs ^ 0.97, h's ~ +0.20, ho ^ -0.10, h'^ ^ 0.07. (2.30) 
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NAN A 

Fig. 4. Schematic picture within a quark model of a Ml (left panel) and E2 (right panel) — > A 
transition induced by the interaction of a photon with a single quark in the nucleon. The Ml transition 
involves a S'-wave spatial wave function, whereas the E2 transition in this picture requires N and/or A 
wave functions to have a D-wave component (indicated by a non-spherical shape). 

From these values it is evident that the S'-wave component dominates the N and A wave func- 
tions in a constituent quark model. The A(1232) resonance is obtained from the nucleon by a 
spin flip of one of the quarks in the Is state so as to give a state of total spin 3/2. Therefore the 
electromagnetic ^ A transition is dominantly a Ml transition [31], see Fig. 4 (left panel). 
Using SU{6) flavor symmetry, i.e., setting as = bs = I and a'g = aM = an = b'g = 6d = 
b'n = in Eqs. (2.27) and (2.28), yields the relation between the magnetic moments of proton 
and p A+ as: 



/ip^A+ = — /^p = 2.63 /X AT ( SU(6)-symmetric quark model). (2.31) 

This is about 25 % lower than the experimental number Eq. (2.24). SU{6) breaking effects, due 
to the hyperfine interaction, decrease the quark model prediction even further to a value around 
2.3 /^AT. A very similar small number is also obtained in the MIT bag model calculation of 
Ref. [32]. The underestimate of /iiv^A in quark models is hinting at the importance of degrees 
of freedom beyond constituents quarks as will be discussed further on. 

The small non-zero values of bo and b'^ in Eq. (2.30) imply that the D-wave probability in 
the A (1232) ground state is around 1 %. As a result of such D-wave components, the N and A 
charge densities become non-spherical. For a static charge distribution, a measure of the non- 
sphericity (or deformation) is given by its quadrupole moment. Since the nucleon has spin 1/2, 
any intrinsic quadrupole moment of the nucleon cannot be directly measured though because 
angular momentum conservation forbids a non-zero matrix element of a (L = 2) quadrupole 
operator between spin 1/2 states. However this quadrupole deformation may reveal itself in an 
electromagnetically induced transition from the spin 1/2 to the spin 3/2 A state, see Fig. 4 
(right panel). In this way, the tensor force between quarks gives rise to non-zero values for the 
E2/M1 (Rem) and C2/M1 (Rsm) ratios. ^ 

Isgur et al. [30], having obtained the D-wave admixture reflected in Eqs. (2.27)-(2.30), find 
Rem — —0.41%. Similar small negative values, in the range —0.8 % < Rem < —0.3 %, were 
found in other non-relativistic quark model calculations [34-36]. Such a small value for Rem 
already indicates that any effect of deformation in the nucleon and/or A ground state is rather 
small and very sensitive to the details of the wave function. 

^ The relation between the tensor force, D-wave admixture, and the electromagnetic ^ A transition 
was already pointed out in the early paper of Glashow [33]. 
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Drechsel and Giannini [37] demonstrated that a calculation of the iV ^ A quadrupole tran- 
sition within the quark model can be very sensitive to a truncation in the quark model basis. 
Due to such truncation (e.g. to a 2hijj oscillator basis as in [30]) the charge-current operator 
for the quadrupole — > A transition violates gauge invariance which can lead to Rem ra- 
tios which differ by an order of magnitude when calculated either from the charge or current 
operator. Drechsel and Giannini showed that when performing a 2hw truncation in the quark 
model basis, a calculation based on the charge operator is more stable than a calculation based 
on the current operator. By using and A 3-quark wave functions which have a D-state prob- 
ability of around 1%, they obtained a value of Rem — —2% and a A quadrupole moment 
Qa ^ -0.078 fm^ [37,38]. 

The error induced due to the truncation in the quark model basis has been further investigated 
in the relativized quark model calculation of Capstick and Karl [39,40]. Using wave functions 
which have been expanded in a large harmonic oscillator basis up to 6hjj states, and are so- 
lutions of a relativized Hamiltonian, they obtained an even smaller negative value: Rem — 
-0.21%. 



2.3.2 Pion cloud models 

Even though the constituent quark model, despite its simplicity, is relatively successful in 
predicting the structure and spectrum of low-lying baryons, it under-predicts Hn-^a by more 
than 25 % and leads to values for the Rem ratio which are typically smaller than experiment. 
More generally, constituent quark models do not satisfy the symmetry properties of the QCD 
Lagrangian. In nature, the up and down (current) quarks are nearly massless. In the exact mass- 
less limit, the QCD Lagrangian is invariant under SU{2)l x SU{2)ji rotations of left (L) and 
right (R) handed quarks in flavor space. This chiral symmetry is spontaneously broken in nature 
leading to the appearance of massless Goldstone modes. For two flavors, there are three Gold- 
stone bosons — pions, which acquire a mass due to the explicit breaking of chiral symmetry by 
the current quark masses. 

Since pions are the lightest hadrons, they dominate the long-distance behavior of hadron 
wave functions and yield characteristic signatures in the low-momentum transfer behavior of 
hadronic form factors. Furthermore, as the A(1232) resonance nearly entirely decays into tiN 
pions are of particular relevance to the electromagnetic A transition. Therefore, a natural 
way to qualitatively improve on the above-mentioned constituent quark models is to include the 
pionic degrees of freedom. 

An early investigation of the 7iVA transition including pionic effects was performed by Kael- 
bermann and Eisenberg [41] within the chiral bag model. The chiral (or, cloudy) bag model 
improves the MIT bag model by introducing an elementary, perturbative pion which couples to 
quarks in the bag in such a way that chiral symmetry is restored [42]. 

Bermuth et al. [43] partially corrected and repeated the calculation of Ref. [41] for two for- 
mulations of the chiral bag model: one with pseudoscalar (PS) surface coupling and one with 
pseudovector (PV) volume coupling of the pions. It was found that using the charge opera- 
tor to calculate the quadrupole amplitude decreases the sensitivity to truncation effects which 
manifest themselves in a violation of gauge invariance in a model calculation. For a bag radius 
around ~ 1.0 fm, the photo-couplings in the chiral bag model are [43]: 



18 



^3/2 = -198 [10-=^ Rem = -1-8%, ( PS ), 

As/2 = -171 [10-=^ GeV-'/% Rem = -2.0%, ( PV ). 

A comparison with experimental data and the predictions of quark models, see Table 1, shows 
that the chiral bag model calculation for PS coupling yields an amplitude that is larger 
than quark models, but still smaller than experiment. For PV coupling, there is a seagull term 
with opposite sign which yields a partial cancellation and deteriorates the agreement with ex- 
periment. 

Lu et al. [44] in a later calculation within the chiral bag model improved upon the previous 
calculations by applying a correction for the center-of-mass motion of the bag. It was found that 
with a smaller bag radius, R ^ 0.8 fm, one is able to obtain a reasonable size for the magnitude 
of the helicity amplitudes, see Table 1. For such a small bag radius, the pionic effects are crucial 
as they account for around 75 % of the total strength of the amplitude ^3/2- 

For the Rem ratio, Lu et al. [44] obtained the value Rem — —0.03%, in disagreement with 
experiment. Both Bermuth et al. and Lu et al. found severe cancellations between the A1/2 and 
A3/2 amplitudes in Eq. (2.14), leading to a very small value for Rem, despite the 'reasonable' 
values for the helicity amplitudes. Bermuth et al. [43] therefore calculated Rem in the Siegert 
limit (|q| — > 0) where it can be obtained from Rsm using the charge operator, and obtained the 
value Rem — —2%, close to experiment. 

In the chiral bag model the quadrupole amplitude E2 has no contribution from the bare bag, 
which is a pure S'-state. Thus, the transition is mediated solely by the pion cloud. To get a more 
realistic value for Rem within the chiral bag model, one probably must either deform the bag 
or carry out a diagonalization like in the quark model (configuration mixing) [45]. 

The pion-cloud contributions also play an important role in the linear a-model and chiral 
chromodielectric model of Ref. [46], even though these models still under-predict the magni- 
tude of the helicity amplitudes, see Table 1. 

The Rem ratio has also been calculated in Skyrme models [47-49], which have only pio- 
nic degrees of freedom. In the Skyrme model, the nucleon appears as a soliton solution of an 
effective nonlinear meson field theory. Wirzba and Weise [47] performed a modified Skyrme 
model calculation, at leading order in the number of colors Nc, based on the chiral effective 
Lagrangian, which corresponds with the standard Skyrme model supplemented by stabilizing 
fourth and sixth order terms in the pion fields. This calculation obtained Rem values between 
—2.5% and —6%, depending on the coupling parameters of the stabilizing terms [47]. Certainly, 
the sign and order of magnitude of these results are consistent with the empirical result. Wal- 
liser and Holzwarth [49] included rotational corrections, which are of order 1/Nc, and lead to 
a quadrupole distortion of the classical soliton solution. Including such corrections, one finds a 
very good description of the photo-couplings (see Table 1) and obtains a ratio Rem — —2.3%, 
consistent with experiment. 

A model which has both quark and pion degrees of freedom and interpolates between a con- 
stituent quark model and the Skyrme model is the chiral quark soliton model (xQSM). This 
model is based on the interaction of quarks with Goldstone bosons resulting from the sponta- 
neous breaking of chiral symmetry. As for the Skyrme model, the xQSM is essentially based 
on the 1/iVc expansion. Its effective chiral action has been derived from the instanton model of 
the QCD vacuum [50], which provides a natural mechanism of chiral symmetry breaking and 
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enables one to generate dynamically the constituent quark mass. 

Although in reality the number of colors Nc = 3, the extreme limit of large Nc is known to 
yield useful insights. At large Nc the nucleon is heavy and can be viewed as "valence" quarks 
bound by a self-consistent pion field (the "soliton") whose energy coincides with the aggregate 
energy of the quarks of the negative-energy Dirac continuum [51]. A successful description of 
static properties of baryons, such as mass splittings, axial constants, magnetic moments, form 
factors, has been achieved (typically at the 30 % level or better, see Ref. [52] for a review of 
early results). After reproducing masses and decay constants in the mesonic sector, the only free 
parameter left to be fixed in the baryonic sector is the constituent quark mass The good agree- 
ment of the xQSM with the empirical situation is achieved for quark mass Mg ~ 420 MeV. 

The xQSM was also applied to the calculation of the 7A'^A transition [53,54]. In this model, 
the A is a bound state which corresponds to a soliton rotating in flavor space. The rotational en- 
ergy is responsible for the — A mass splitting. With a constituent quark mass Mq = 420 MeV, 
the empirical — A mass splitting is well reproduced. The 7A^A amplitudes are obtained tak- 
ing rotational {l/N^ corrections into account. In this way, Silva et al. [54] obtained as result: 
Rem — —2.1% (in the two-flavor case), fairly close to the experimental ratio, considering that 
in the xQSM calculation no parametrization adjustment has been made to the N A tran- 
sition. However, the value of the Ml 7A^A amplitude is largely underpredicted in the xQSM, 
which is also reflected in an underprediction of the magnitude of the photo-couplings, see Ta- 
ble 1. 

The above calculations incorporating the chiral symmetry of QCD (to lowest order in the pion 
fields) highlight the role of the pionic degrees of freedom in the 7A^A transition. A number of 
subsequent works have therefore revisited quark models, which break chiral symmetry, by in- 
cluding two-body exchange currents. As an example of such approach we consider the work of 
Buchmann, Hernandez and Faessler [55]. 

When the one-gluon exchange potential between quarks is complemented by a one-pion and 
one-sigma exchange potential between the quarks, the pion exchange gives rise to an additional 
tensor interaction between the quarks. The one-pion exchange then requires the presence of 
two-body exchange currents between the quarks due to current conservation. Within the non- 
relativistic framework, Buchmann et al. [55] found that the overall effect of the exchange cur- 
rents on the dominant Ml 7A^A transition are relatively small due to cancellation. As a result, 
the transition magnetic moment comes out to be underpredicted: yUp^A+ — 2.5 [/ijv], similar to 
the non-relativistic quark model without exchange currents. This is also reflected in the helicity 
amplitudes obtained in [55] and quoted in Table 1. They are smaller in the magnitude than the 
experimental values. 

However, the two-body exchange currents lead to non-vanishing 7A^A quadrupole (E2 or 
C2) amplitudes [55], even if the quark wave functions have no D-state admixture. In this pic- 
ture, the A is excited by flipping the spins of two quarks, see Fig. 5. According to Buchmann et 
al, this mechanism yields Rem — —3.5%. 

Also within their quark model, Buchmann et al. related the — > A and A+ quadrupole 
moments to the neutron charge radius as ^ : 



^ See also Refs. [56,57] for an analogous relation between the N —>■ A quadrupole form factor G^, and 
the neutron electric form factor. Possible higher order corrections to Eq. (2.32) allowed by QCD were 
recently investigated in Ref. [58] within a general parametrization method, see also Ref. [25]. 
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Fig. 5. Schematic picture within a quark model of a E2 N ^ A transition induced by the interaction of 
a photon with two quarks in the nucleon (correlated through, e.g., one-pion exchange). Both and A 
wave functions are S'-states. 

Qp^A+ = -^^l^ = ^l- (2.32) 

The origin of these relations in the model calculation of Ref. [55] is the dominance of two-body 
exchange currents in both Qp^/\+, Qa+ and r^. Using the experimental neutron charge radius, 
rl = -0.113(3) fm^ the relations of Eq. (2.32) yield: 



Qp^A+ = -0.08 fm^ Qa+ = -0.113 fml (2.33) 

The value of Qp-,A+ is close to the empirical determination of Eq. (2.25). 

Recently, the 7*A^A transition form factors were calculated in a Lorentz covariant chiral 
quark approach [62], where baryons are modeled as bound states of constituent quarks dressed 
by a cloud of pseudoscalar mesons. To calculate baryon matrix elements, a phenomenological 
hadron-quark vertex form factor was used, which parametrizes the distribution of quarks inside 
a given baryon. In Ref. [62], a Gaussian form was adopted for this hadron-quark form factor 
which involves a size parameter A^, which was used as a free parameter. Using a value ~ 
0.8 GeV, it was found that the contribution of the meson cloud to the static properties of light 
baryons is up to 20 % and, together with the relativistic corrections, helps to explain how the 
shortfall in the SU(6) prediction for hn^a is ameliorated. 

As a summary, we list in Table 1 the 7A^A photo-couplings A1/2 and A3/2 as well as the ratio 
Rem the various models discussed above. 



2.3.3 Intrinsic quadrupole moment and more about shape 

The relations of Eq. (2.32) were further interpreted by Buchmann and Henley [59] in terms of 
an intrinsic quadrupole moment of the nucleon and A states. The intrinsic quadrupole moment 
of a static charge distribution p(r) is given by: 



qo = y d^rp{r) (3^2 _ ^2^^ ^2.34) 

which is defined w.r.t. the body fixed frame. A charge distribution concentrated along the z- 
axis (symmetry axis of the system) corresponds with > (prolate deformation), whereas 
a charge distribution concentrated in the equatorial a;?/-plane corresponds with (5° < (oblate 
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A. ,n [10-3 GeV-^/^l 


ylo/9 flO-3 GeV-^/^l 


Rfm [%] 


experiment 
Ref. [15] 


—135 ± 6 


—250 ± 8 


—2.5 ± 0.5 


SU(6) symmetry 


-107 


-185 





non-relativistic quark models 

i?i:»-fo TOO '2A 
Keis. [zy, ju,j4— J / J 


1 (YX 
-Wji 


-1 /y 


— z to U 


relativized quark model 
Refs. [39,40] 


-108 


-186 


-0.2 


MIT bag model, Ref. [32] 


-102 


-176 





chiral (cloudy) bag models 
Kei. [4jj : (ro, K = 1 imj 
Kei. L^-jj . V, K — 1 imj 
Ref. [44] : (+ recoil, R = 0.8 fm) 


1 

-lUo 

Q1 

-yi 
-134 


1 oo 
1 71 

-1/1 
-233 


— i.O 

-0.03 


linear cr-model, Ref. [46] 


-107 


-199 


-1.8 


chromodielectric model, Ref. [46] 


-/U 


1 '2 1 


1 o 

-i.y 


Skyrme models 
Ret. [4 /J 
Ref. [49] : (+ l/N^ corrections) 


-136 


-259 


— D to —2.0 

-2.3 


chiral quark soUton model 
Ref. [54] : SU(2) flavor 


-70.5 


-133 


-2.1 


niiarlc modpl + tt n pxrhanpp 
Ref. [55] 


-91 


-182 


-3.5 


chiral quark model (Ag = 0.8 GeV) 

Ref. [62] 


-124.3 


-244.7 


-3.1 



Table 1 

Summary of the values of the 7A'^A helicity amphtudes A1/2 and at = 0, and the ratio 
REMiQ"^ = 0) in different models compared with experiment. 



deformation). This intrinsic quadrupole moment has to be distinguished from a measured (or 
spectroscopic) quadrupole moment Q. As an example, for a rigid rotor (which was considered 
within the context of the collective nuclear shell model [61]) these quantities are related as: 
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Fig. 6. Model picture of the intrinsic quadrupole deformation of the nucleon (left) and A (right) within 
the pion cloud model of Ref. [59]. In the the p-wave pion cloud is concentrated along the polar 
(symmetry) axis leading to a prolate deformation, whereas in the A, the pion cloud is concentrated in the 
equatorial plane, corresponding to an oblate deformation. Figure from Ref. [59]. 



Q 



3 J| - J( J + 1) 
(J + l)(2J + 3) 



(2.35) 



where Jz is the projection of the nucleon total spin J onto its symmetry axis (2;-axis in a body 
fixed frame). The difference between and Q represents the averaging of the nonspherical 
charge distribution due to its rotational motion as seen in the laboratory frame. One verifies 
from Eq. (2.35) that the multiplication factor is zero for a spin 1/2 particle, yielding Qp — 
for the proton. Eq. (2.35) does not preclude however that the proton has an intrinsic quadrupole 
moment Q^. 

Within a model where the nucleon consists of a spherically symmetric quark core surrounded 
by a pion with orbital angular momentum 1 = 1, Buchmann and Henley obtain [59]: 



QA+=rn = -Q?- (2.36) 

Thus, the proton and A+ have respectively a prolate and an oblate intrinsic deformation. In this 
hybrid (quark/pion-cloud) model, the pion cloud is fully responsible for the non-zero values of 
the quadrupole moments and Q%+, and hence for the non- spherical shape of these particles, 
see Fig. 6. 

Estimates such as Eq. (2.36) of intrinsic quadrupole moments are surely useful to reveal 
details of a given model calculation and gain physical insight. One should keep in mind, how- 
ever, that, even though a specific model such as the rigid rotor relates the intrinsic and total 
quadrupole moments, only the latter is directly related to observables. 



2.4 Large Nc limit 



Though the results from the different QCD inspired models discussed above may provide us 
with physical insight on the 7 A^^A transition, they are not a rigorous consequence of QCD. In the 
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following subsections, we will discuss what is known on the 7A^A transition from approaches 
which are a direct consequence of QCD in some limit, such as chiral effective field theory (chiral 
limit of small pion masses or momentum transfers) or lattice QCD simulations (continuum 
limit). 

Another first-principle technique is the 1/Nc expansion of QCD proposed by 't Hooft [63] 
and Witten [64] . The large- A^c limit is unique in that it provides a perturbative parameter at all 
energy scales. This expansion has proved quite useful in describing the properties of baryons, 
such as, ground-state and excited baryon masses, magnetic moments, electromagnetic decays, 
see e.g., Refs. [65,66] for reviews. 

In the large Nc limit, the baryons are infinitely heavy and can be treated as static. In this limit 
the baryon sector of QCD has an exact contracted SU{2Nf) spin-flavor symmetry, where Nf 
is the number of light quark flavors. The large Nc limit thus validates many of the quark model 
SU (6) spin-flavor symmetry results without making any model assumption (such as assuming 
non-relativistic quark dynamics for the baryon wave functions). For example, the ratio of proton 
to neutron magnetic moments in the large- A^c liniit is predicted to be given by Hp/ = — 3/2 in 
agreement with the naive quark model. Likewise at leading order in 1/ Nc, the N ^ A transition 
magnetic moment una is related to the isovector combination of proton and neutron magnetic 
moments as [67]: 

up to a correction of order 1 /N^. Using the empirical values for jip and /i„, one obtains in the 
large Nc limit ^n^a = 3.23 fi^, within 10 % of the experimental value of Eq. (2.24). 

The Rem ratio for the 'yNA transition is shown to be of order l/N^ [68]. Thus, the smallness 
of the 7A'^A Rem ratio is naturally explained in the large Nc limit. Using Eq. (2.14), this can 
equivalently be expressed as a large Nc prediction for the 7A^A helicity amplitudes: 



The large Nc limit also allows to obtain relations between the A and N ^ A quadrupole 
moments [69]: 



This result for Qa+ is different from the quark model ratio of Eq. (2.32) as explained in [69]. 
Using the phenomenological value of Eq. (2.25) for Qp-^A+, the large Nc relation of Eq. (2.39) 
yields for the A+ quadrupole moment: 




(2.38) 




(2.39) 



= - (0.048 ± 0.002) fm^ 



(2.40) 
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accurate up to corrections of order 

Buchmann, Hester and Lebed [69] derived recently another relation by making the additional 
assumption that the baryon charge radii and quadrupole moments arise from the same con- 
tributions. This assumption holds true when, e.g., the interaction between quarks arises from 
one-gluon or one-pion exchange, for which the ratio of spin-spin and tensor terms in the Hamil- 
tonian is fixed. Under this assumption, one finds a large- A^c relation between the N ^ A 
quadrupole moment and the neutron charge radius [69]: 



2 IN, + 5 



where both Qp^A+ and are of order l/N^. 

r.2 : 



The factor on the rhs of Eq. (2.41) after is unity both for the cases Nc — 3 and Nc — > oo, 
and deviates from unity by only about 1 % for all values in between. Therefore the large limit 
predicts to good approximation the same relation of Eq. (2.32), Qp^A+ = 1 /v^^n^ which was 
also derived in the quark model and pion cloud model, and which was found to be consistent 
with the empirical value of Eq. (2.25) for Qp^A+- 

More recently, Cohen and Lebed [70,71] have argued that the relations, such as in Eqs. (2.39) 
and (2.41), may be significantly altered in the real world where the pions are light. When taking 
iVc ^ oo at a fixed value of m^, the A^-A mass difference, which goes as l/N,, vanishes. 
Therefore, the "usual" large- Ai'c limit implies that one is in the region Ma — Mjy <^ m^, where, 
for instance, the A is stable. Cohen points out [70] that the large-iVg limit and the chiral limit 
(m^r — > 0) do not commute and for quantities which diverge in the chiral limit, such as the 
charge radii, one expects chiral corrections to dominate over the large- A^c predictions. 

Although the region — Mj^ -C is not accessible in nature, lattice QCD results, which 
are currently obtained in this regime (see Sect. 2.6), may provide an interesting testing ground 
for the above large predictions. 
Let us conclude the discussion of the large- A^'c relations by observing a new one, i.e. : 

Rsm^Rem, for g2 ^ and A^c ^ oo . (2.42) 

To derive this relation we follow the arguments given in the Appendix of Ref. [72]. Namely, we 
first make use of the expression of the ratios in terms of the little g's of Eq. (2.2), at Q"^ — 0: 



Rem = A ^Tnr~rjr^9E , (2.43a) 



R 



SM ■ 



9m - 2{Mn+M^)9E 



9e + ttTTSc 



2{Mn + MaY^ ami 



(2.43b) 



where A = M\ — is the A^-A mass difference. These are general expressions, which can, 
e.g., be obtained by substituting Eq. (2.3) into Eq. (2.15). Now, the little (yf's are the coupling 
constants from an effective Lagrangian (see Sect. 4.3) and, as can be seen from the similar 
operator structure of the corresponding terms, these constants have the same large- A^c scaling. 
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[a) 



Fig. 7. The vrA^ loop corrections to the 7pA+ vertex as calculated in xEFT. Diagram (a) : vr+n loop 
where the photon couples to the 7r+; diagram (b) : ir^p loop where the photon couples to the charge of 
the proton. 

Finally, we only need to recall that the baryon masses go as N^, while the A^-A mass difference 
as 

Mn, Ma ~ A/'c, ~ for A^'^ ^ oo, (2.44) 

to see that the gc term in Rsm is of higher 1/Nc order, and hence the relation Eq. (2.42) holds. 
From Eq. (2.43), one immediately verifies that both of these quantities are of order l/N^- 



2.5 Chiral effective field theory 

The A transition has also been studied within the chiral effective field theory (^EFT) 

expansions based on chiral Lagrangians with nucleon and A-isobar fields. In this framework 
the expansion for each of the 7*A^A-transition form factors begins with a low-energy constant 
(LEG), which then receives the chiral loop corrections of the type depicted in Fig. 7. 

A first such study was performed by Butler, Savage and Springer [83] in the framework of 
heavy-baryon ^PT [84]. At leading order they obtained a "chiral log" (i.e., Inm^r) enhancement 
of the E2 transition, which lead to relatively large and positive values of Rem- Their result thus 
showed that the chiral correction to the E2 transition diverges in the chiral limit. 

A more comprehensive study was subsequently carried out by Gellas et al. [85] using the 
so-called "small scale expansion" (SSE) [86], also called "e-expansion". In the SSE scheme 
the two light scales in the problem: the pion mass (m^r) and the A-resonance excitation en- 
ergy (A = Ma — Mn) are counted as having the same size. In addition, the heavy-baryon 
(semi-relativistic) expansion is performed. The calculation of Ref. [85] was performed to order 
p^, which is the next-to-leading order for Ml and leading order for the E2 and C2 transi- 
tions. More recently, Gail and Hemmert [87] have updated this work and analyzed some of the 
higher-order contributions. In overall, they show a good agreement with the low data for the 
7A^ A transition, as will be illustrated below. However, as a note of caution, one should point 
out that the Rem and Rsm ratios are reasonably well described only upon adding a nominally 
higher-order (p^) counter term with an unnaturally large LEG (Cg ~ — 17 GeV^^ in the notation 
of Ref. [87]). 

Two of us [88,89] have recently computed the 7*A^A form factors in a manifestly covariant 
(no heavy-baryon expansions) xEFT expansion. In the chiral limit, it was found that the C2 
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Fig. 8. The pion mass dependence of Rem and Rsm for three fixed values of Q^: (solid curves), 
0.025 GeV^ (dashed curves), and 0.1 GeV^ (dotted curves), resulting from the xEFT calculation of 
Refs. [88,89]. The figure shows the absence (presence) of the him^ divergence for Rem {Rsm) at 
Q2 = o. 

transition diverges for = 0, which is also in agreement with the result of Ref. [87] ^ . The 
coefficient of the chiral log turns out to be small. Only for tiny pion masses one starts to see the 
dominance of the chiral log in the C2 transition, see Fig. 8. 

Our xEFT study of the 7A^A transition [88,89] utilizes the "(5-expansion" scheme of Ref. [90]. 
In this scheme the two light scales 5 = Z\/A^sb and e = niT^/A^sB, with A^sb ~ 1 GeV the chi- 
ral symmetry breaking scale, are treated differently. In contrast, the SSE assumes e ~ 5, which 
leads to an unsatisfactory feature that the A-resonance contributions are always estimated to be 
of the same size as the nucleon contributions (hence, overestimating the A contribution at lower 
energies and underestimating some of them at the resonance energies). In the 5-expansion, one 
counts e = 5^, which is the closest integer-power relation between these parameters in the real 
world. 

At the physical pion mass, the 5-expansion provides an energy-dependent power-counting 
scheme designed to take into account the large variation of the A-resonance contributions with 
energy. As such it allows for an efficient calculation of observables in the resonance region. The 
relevant LECs can in this fashion be directly extracted from observables. In Refs. [88,89] the 

^ This result of both Refs. [87] and [88,89] is in disagreement with Butler et al. [83] where a chiral log 
enhancement was reported for E2 at = instead of C2. 
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Fig. 9. The low dependence of the quadmpole ratios EMR (left panel) and CMR (right panel) 
for the 7*p A transition. The data are from MAMI : Refs. Q [23], [77], S [78], T [79]; from 
Bates A [80]; and from CLAS □ [81]. The lattice QCD calculations with linear pion mass extrapolations 
are shown as x [82]. Two chiral EFT calculations are shown: the (^'-expansion result from Ref. [88,89] 
(red solid curves with error estimates), and the e-expansion result of Ref. [87] (black solid curves). The 
dynamical model predictions from SL [14] (green dashed-double-dotted curves) and DMT [91] (blue 
dashed-dotted curves) are shown alongside the phenomenological MAID2003 [92] (red dotted curves) 
and SAID [93] (black dashed-triple-dotted curves) models. The hypercentral (long dashed curves) [94] 
and relativistic (short dashed curves) [39] constituent quark models have been included. Figure from 
Ref. [79]. 



three LECs corresponding to the three 7A^A transitions are all extracted from a next-to-leading 
(]S[LO) calculation of pion photo- and electro-production observables (see Sect. 4.5.2). The re- 
sulting prediction of the dependence is then shown to be in a qualitative agreement with 
phenomenological extractions. 

Namely, in Fig. 9, we show the present experimental situation for the Rem and Rsm ratios at 
low and compare the data with xEFT calculations, dynamical and phenomenological model 
predictions, as well as quark model predictions. One sees from Fig. 9, that the Rem ratio stays 
small and negative, around —2 to —3 % up to momentum-transfer of around 0.5 GeV^. 

The Rsm ratio displays a steep slope at low (Siegert limit), and seems to level off to a 
negative value around —6 to —7 % above ~ 0.1 GeV^. Both the xEFT calculations are con- 
sistent with the low Q"^ dependence of the Rem and Rsm ratios. The 5-expansion results are 
quoted with a theoretical uncertainty band as will be explained in Sect. 4. Furthermore, one sees 
from Fig. 9 that the constituent quark models largely under-predict the Rem and Rsm ratios, 
as discussed above. The dynamical models of Refs. [14,91], which include pionic degrees of 
freedom, are in qualitative agreement with the empirical dependence. 

In Fig. 10, we compare Q^-dependence of the the pion cloud contributions to the 7* A^A form 
factors in dynamical models versus xEFT. In xEFT the values of the form factors G\j, G%, and 
G*Q at the real photon point are related to three LECs qm, Qe, and gc appearing in the chiral La- 
grangian (see Sect. 4), and determined from a fit to the data. Similarly, in the dynamical models 
the values at = are determined from a fit to the same data. However, the renormalization 
of the pion loop contributions is done in a dramatically different fashion. In xEFT one uses 
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Fig. 10. Comparison of the pion cloud in the dynamical model and chiral EFT calculations for the low 
behavior of the Jones Scadron form factors G\i, G*^, and (with value at = subtracted). 
The dynamical model predictions from SL [14] are shown for both the 'bare' version of the model 
(blue dashed curves) where the pion loop diagrams are not accounted for, and 'dressed' version of the 
model (black dashed dotted curves) which include the pion loop diagrams. The chiral EFT calculation of 
Ref. [88,89] are shown without (green dotted curves) and with (red solid curves) the pion loop diagrams. 

the dimensional regularization and absorbs infinities and other power-counting violating pieces 
into the LECs. In dynamical models, the loops are regularized by a cutoff, which is then fitted 
to the data. As the result the pion cloud contributions do not satisfy the chiral power-counting, 
and are usually larger than the ones obtained from xEFT. 

Therefore, to compare the dynamical model and xEFT results directly we consider the differ- 
ence G\.j{Q'^) — G'Jj(O), and analogously for G% and G^, such that the above renormalization- 
scheme dependence drops out. In Fig. 10, we compare the results of the dynamical model of 
Ref. [14] with the xEFT results of Ref. [88,89]. For both the dynamical model and xEFT we 
show the results with and without the pion-loop effects, which for these difference-quantities 
appear to be similarly renormalized. 

From this comparison, we conclude that the full result, including the pion-loop contributions, 
is very similar in both the dynamical model and the xEFT approach, for all the three form fac- 
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A(p') N(p) 

Fig. 11. The 7* A^A matrix element as evaluated in lattice QCD. The photon couples to one of the quarks 
in the nucleon at a fixed time ti to produce a A. 

tors. The pion cloud contribution to G\j appears to be larger in the dynamical model than in the 
XEFT calculation. This difference in size is mainly due to the photon coupling to the anomalous 
magnetic moments of the nucleons in the loops, which are included in the dynamical model, but 
not in xEFT since they are of higher than NLO considered in Refs. [88,89]. For the Coulomb 
quadrupole form factor Gq, one sees a very similar role of the pion cloud in both approaches, 
giving rise to a structure around Q"^ ^ 0.03 — 0.05 GeV^. For the electric quadrupole form 
factor G*^, on the other hand, one notices a qualitatively different behavior. Whereas the effect 
of the pion cloud is very small in the dynamical model calculation, a strong, nearly linear Q'^ 
dependence of the pion loops arises in xEFT. 

The role of the pion cloud will further be elucidated within the dynamical models in Sect. 3, 
and xEFT in Sect. 4. 

2.6 Lattice QCD and chiral extrapolation 

2.6.1 Lattice simulations 

Lattice QCD calculations of nucleon structure quantities have matured considerably in the 
recent past. They provide an ab initio calculation of quantities such as the 7*A^A transition 
form factors from the underlying theory of the strong interaction. The first such calculation was 
performed by Leinweber, Draper and Woloshyn [73]. 

The calculation of the 7*A^A form factors requires the evaluation of three-point functions, 
which involve the computation of a sequential propagator, see Fig. 11. Leinweber et al. [73] 
evaluated these in the so-called fixed current approach, which requires the current to have a 
fixed direction and to carry a fixed momentum. The initial and final states, on the other hand, can 
vary without requiring further inversions, which are the time-consuming part of the evaluation 
of three-point functions. For rather large quark masses, corresponding with pion masses in the 
range m,, ~ 0.65 — 0.95 GeV, they obtained [73]: Rem = (—3 ± 8)%. This initial result clearly 
indicates the need for high statistics to establish a non-zero value for a quantity such as Rem- 

Alexandrou, de Forcrand and Tsapalis [74] studied the deformation of baryons with spin 
higher than 1/2 in lattice QCD via three-density correlators. Such three-density correlators, 
when considered at a fixed time t, are function of the two relative distances ri and r2 between 
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Fig. 12. Three-density correlator of a baryon as defined in Eq. (2.45) and used in the lattice gauge cal- 
culations of Ref. [74]. The source (right) and sink (left) are taken at the same spatial location, and are 
separated by a Euclidean time distance T/2. The crosses indicate the density insertions (p"') into the 
up (down) quark lines respectively. 

the three quarks in a baryon B and are defined as: 



where the density insertion p'^ (q = u, d) into an up (u) or down (d) quark line is given by the 
normal order product: 



so that disconnected graphs are excluded. These correlation functions of quark densities, which, 
being expectation values of local operators, are gauge invariant, reduce in the non-relativistic 
limit to the wave function squared. The three-density correlator for a baryon is shown schemat- 
ically in^ Fig. 12. For simplicity, the baryon source and sink are taken at the same spatial 
location. Denoting the lattice extent in Euclidean time by T and taking periodic boundary condi- 
tions, yields a maximum time separation between source and sink of T/2. To isolate the ground 
state baryon, both the time distances between the source and the insertions (t) and between the 
sink and the insertions (T/2 — t) are taken to be as large as possible in such calculations. 

In Ref. [74], the correlator has been calculated in both quenched and unquenched lattice QCD 
for a in a spin -1-3/2 state. The quenched calculation was performed for pion masses in the 
range m.,, ~ 0.35 — 0.65 GeV, and no signal of deformation has been observed. This may not 
come as a surprise, because in a quenched calculation, where graphs with disconnected quark 
loops are neglected, pion cloud effects are only partially accounted for. It was discussed above 
in the context of phenomenological models that pion cloud effects are partly responsible for 
hadron deformation. This was checked by performing an unquenched lattice QCD calculation 
with two heavy dynamical quarks, using the Wilson Dirac operator in the same quark mass 
range as the quenched calculations [74]. The resulting correlator for the is displayed in a 

^ In addition to one-density insertions into each n-quark line, one can also have two-density insertions 
into the same u-quark line. The latter contribution was however checked to be small [74]. 



C(ri, r2, t) = I dh' {B I p"(r', t) p"(r' + n, t) p"(r' + r^, t)\B), 



(2.45) 



p%r,t) =: g(r,t)7%(r,t) : 



(2.46) 
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Fig. 13. Three-dimensional contour plot of the correlator (red) for the state with spin projection +3/2 
(the z-axis denotes the spin direction) as it follows from lattice QCD calculations [74] for two heavy 
dynamical quarks, corresponding with a pion mass m,r > 0.6 GeV. The contour (red), corresponding 
with a value (0.8) for the correlator which has been chosen to show large distances, indicates a slightly 
oblate deformation of the A+ in the +3/2 spin state. For comparison, the contour of a sphere is also 
shown (green). Figure from Ref. [74]. 

three-dimensional contour plot in Fig. 13. The unquenched lattice QCD calculation indicates a 
slightly oblate deformation for the in a spin +3/2 state, in agreement with the negative value 
of the A+ quadrupole moment obtained from the large Nc estimate of Eq. (2.40). It will be very 
interesting to check if this signal can be consolidated by higher statistics results using lighter 
quark masses and a larger lattice. 

A direct calculation of the 7*A^A transition form factors on the lattice has also been per- 
formed recently. Based on the method of Ref. [73], the Nicosia-MIT group [75] has performed 
a high- statistics calculations of the 7* A^A form factors. The Nicosia-JVIIT group has also imple- 
mented several improvements of the algorithm, such as smearing techniques to effectively filter 
the ground state, check on the volume dependence of the results, use of a larger lattice such as 
to simulate smaller quark masses. Both quenched and unquenched results have been obtained. 
The quenched calculation was performed for pion masses in the range 0.4 — 0.9 GeV. 

The unquenched calculation, using dynamical Wilson fermions, was performed in the range of 
^ 0.5 — 0.8 GeV. At such pion mass values, both the quenched and unquenched calculations 
give Rem values in the range from —1 to —5 %. For pion masses in the range 0.5 — 0.8 GeV 
unquenching effects were found to be within statistical errors, which means that pion cloud 
contributions, expected to drive the Rem value more negative are suppressed for these large 
quark (pion) masses. 

In order to compare the calculated Rem values with experiment, the Nicosia-IMIT group uses 
a linear fit in (corresponding with a linear fit in the quark mass) was used to extrapolate 
down to the physical pion mass. Such a naive extrapolation can not be justified, as will be dis- 
cussed below. Nonetheless, using such a linear extrapolation, the Rem value at Q'^ = 0.4 GeV^ 
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Fig. 14. Quenched lattice QCD results for the dependence of the (Ash) form factor 
^Mi ~ Ash/3' as defined in Eq. (2.18) for three values of m-,^ as indicated on the figure. The 
(blue) solid circles are the result of a linear extrapolation in m^. The experimental value at = (red 
square, PDG) corresponds with the value given in Eq. (2.24). The lower solid curve corresponds with a 
fit to the data, according to the MAID analysis [26]. Figure from Ref. [82]. 

was found to be [75] : 



which is to be compared with the experimental number of Ref. [81]: Rem = (—3.4 ± 0.4 ± 
0.4) %. One sees that the lattice QCD calculation clearly supports a negative value of Rem^ in 
basic agreement with experiment. 

In the subsequent worlc of the Nicosia-IVIIT group [82], the 7*A^A form factors were stud- 
ied within a fixed sinlc method in which the initial state, created at time zero, has the nucleon 
quantum numbers, and the final state, annihilated at a later time t2 has the A quantum numbers. 
The current can couple to a quarlt line at any time slice ti, see Fig. 11, carrying any possible 
value of the lattice momentum. Implementing further improvements, the method of Ref. [82] 
is superior to the fixed current approach discussed above, yielding a more accurate evaluation 
of all three 7*A^A form factors. The calculations in Ref. [82] were performed in the quenched 
approximation for pion masses: = 0.51, 0.45, and 0.37 GeV, which are also smaller than in 
the previous calculations. 

The quenched lattice QCD results of Ref. [82] for the magnetic 7*A^A form factor are shown 
in Fig. 15 together with a linear extrapolation in of the lattice results. The lattice results 
show a decrease of G*^ with decreasing pion mass. However the linear extrapolated values in 
are still lying significantly above the experimental results. The most likely explanation of 
this discrepancy is the inadequacy in the linear extrapolation. At low momentum transfer, where 
chiral effective field theory can be applied, one expects chiral logs to appear in the form factors, 
as will be discussed further on. 

In Fig. 15, the quenched lattice QCD results of Ref. [82] are shown for the 7*iVA ratios Rem 



Rem{Q^ = 0.4 GeV^) = (-4.4 ± 1.7) % 
Rem{Q^ = 0.4 GeV') = (-3.7 ± 1.2) % 



(quenched) , 
(unquenched) , 
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Fig. 15. Quenched lattice QCD results for the dependence of the ^*NA ratios Rem (left panels) and 
RsM (right panels), for three values of ttItt as indicated on the figure. The (blue) solid circles are the 
result of a linear extrapolation in m^. In the lower panels a comparison with experimental data is shown: 
real photon points from MAMI [23] and LEGS [25]; finite values are from BATES [95,80] (solid 
triangles), MAMI [78] (cross), and JLab/CLAS [81] (open boxes). Figure from Ref. [82]. 

and RsM, together with a linear extrapolation in m^. For the Rem ratio, the lattice results are 
accurate enough to show a negative value, which becomes more negative as one approaches 
the chiral limit. The linearly extrapolated lattice results seem to be in good agreement with the 
experiment, although the lattice results at larger Q"^ values show some scatter. The Rsm ratio is 
clearly negative over the whole range. For > 0.5 GeV^, the linearly extrapolated lattice 
results are in agreement with experiment. At lower on the other hand, the linearly extrap- 
olated lattice results fall increasingly short of the data, and cannot explain the large negative 
value of Rsm established in experiment. The present empirical results for Rsm at low have 
reached a high level of accuracy and have been cross-checked by several experiments at both 
BATES [95,315,80] and IVIAJMI [78,77], all obtain a rather large negative Rsm ratio of around - 
6% at low Q^. This puzzle was studied in Refs. [87-89] within the framework of xEFT, yielding 
some interesting results as is discussed in the following. 

2.6.2 Chiral extrapolations 

The extrapolation in the quark mass iriq is not straightforward, because the non-analytic de- 
pendencies, such as ^ym^ and Inm^, become important as one approaches the small physical 
value of niq. Therefore naive extrapolations often fail, while spectacular non-analytic effects 
are found in a number of different quantities, see e.g., Refs. [96,97]. The xEFT, discussed in 
the previous section, provides a framework to compute these non-analytic terms. We will first 
address the quark mass dependence of the nucleon and A masses within xEFT, which have been 
discussed extensively in the literature, see e.g. Refs. [98-105]. Subsequently, we will review the 
predictions for the pion mass dependence of the 7*A^A form factors. 

In the xEFT including A degrees of freedom, the pion mass dependence of and A masses 
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Fig. 16. The nucleon and A self-energy contributions in xEFT. Dashed lines represent pion propagators, 
solid lines represent nucleon propagators, whereas double lines represent A propagators. 

is calculated from the self-energy diagrams of Fig. 16. Near the chiral limit, i.e. for < A = 
Ma — Mjv, the pion mass dependence for the N and A masses goes as [98]: 

= - 4 c,^ ml - ^^gl ml + J^^h^ ^ Inm, + 0{mt), (2.47) 
Ma = Mf - 4 c,A ml - ^1^1 ml - ^^^,hl ^ In + 0(mt). (2.48) 

The non-analytic terms (proportional to ml, m^/Z\ In m^,...) on the rhs of Eqs. (2.47,2.48) are 
predictions of xEFT obtained from the loop diagrams of Fig. 16. These loop expressions de- 
pend on the coupling constants appearing in the lowest order chiral Lagrangian (see Sect. 4 
for details): qa = 1.267 is the axial coupling of the nucleon, Ha — 2.85 is the ttA^A cou- 
pling constant (in the notation of Ref. [105]), and Ha is the axial ttAA coupling constant, 
which is related with qa through the SU{6) relation, which coincides with the large Nc rela- 
tion: Ha = {9/5)gA — 2.28. The analytic terms on the rhs of Eqs. (2.47,2.48) are low-energy 
constants which have to be determined from experiment or from a fit to lattice QCD results. 
In particular M^ (M^) is the N (A) mass in the chiral limit, and the term proportional to cm 
(cia) is the quark mass contribution to the N (A) mass. For the nucleon, it is obtained from the 
experimental information on the pion-nucleon a-term. 

The above formulas of Eqs. (2.47,2.48) can be fitted to full lattice QCD results for the N and 
A masses. Care has to be taken however when fitting to quenched lattice QCD results, where sea 
quark loop effects are neglected. The quenched approximation also modifies the leading chiral 
expansion of baryon masses in the corresponding effective field theory, see Ref. [106]. For in- 
stance, in the quenched approximation, the pion mass dependence of N and A masses contain a 
term linear in m^, whereas such a term is absent in the full QCD expansions of Eqs. (2.47,2.48). 

In Ref. [100], the m^, dependence of N and A masses have been fitted to lattice results using 
both the quenched theory and full QCD. The non-analytic terms have been calculated from the 
one-loop diagrams of Fig. 16, including a phenomenological form factor at the nNN, ttNA and 
ttAA vertices to account for the finite size of the pion source. The parameters in the analytic 
terms were treated as free parameters and fit to the lattice results. Such a fit, is shown in Fig. 17. 
One clearly sees that this procedure is able to successfully account for the different behavior of 
the dependence of and A masses in quenched as compared to full QCD. One notices that 
the effect of unquenching is larger for the A as compared to the N. It will be interesting to test 
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Fig. 17. Pion-mass dependence of the N and A masses according to Ref. [100] compared with lattice 
QCD results. The soUd symbols are lattice results from the MILC Collaboration [107]. Solid symbols 
from bottom to top correspond with the nucleon unquenched, nucleon quenched, A unquenched, and A 
quenched lattice results. The open boxes are the fit of Ref. [100] to and A masses accounting for finite 
volume and lattice spacing artifacts, as described in the text. The curves describe the infinite-volume 
limit for dynamical (soUd curves) and quenched (dashed curves) N and A masses. 

the more singular behavior of the quenched result for the A as lattice results become available 
for smaller pion masses. 

A comparison of both quenched and full QCD also allows to draw interesting conclusions on 
the physical nature of the — A mass splitting, because the effect of the pion cloud is only 
partly accounted for in the quenched theory. Ref. [100] found that only about 50 MeV of the 
observed 300 IVIeV N — A mass splitting arises from the pion cloud, while the rest arises from 
short distance processes, such as gluon exchange. This short distance contribution is responsible 
for the finite — A mass splitting in the chiral limit (i.e. when — > 0). 

To account for higher order terms in the chiral expansion, the dependence of Mjv and 
Ma has been studied by two of us [105] in a manifestly covariant xEFT framework consistent 
with analyticity. The resulting relativistic loop corrections obey the chiral power-counting, after 
renormalizations of the available counter-terms are done. The relativistic expressions also con- 
tain the nominally higher-order terms, which are necessary to satisfy the analyticity constraint. 
In such approach, the analytic terms in the quark mass arising from the one-loop pion diagrams 
of Fig. 16 are partially resummed. We found [105] that the convergence of the chiral expansion 
for Mat and Ma can be improved without introducing additional parameters. In Fig. 18, the 
results for the m^r dependence of N and A-resonance masses in this approach are compared 
with full lattice QCD results. For both N and A, the chiral limit mass values M^°^ (Mi°^) and 
the parameters cin (cia) are fitted. As is seen from the figure, with this two-parameter form for 
Mjv and Ma, a good description of lattice results is obtained up to ~ 0.5 GeV^. 

Having discussed the dependence of Mjv and Ma, we turn to the dependence of the 
YNA transition form factors. The study of the rrij^ -dependence is crucial to connect to lattice 
QCD results, which at present can only be obtained for larger pion masses (typically > 0.3 
GeV) as discussed above in Fig. 14 for G^, and in Fig. 15 for Rem and Rsm- 
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Fig. 18. Pion-mass dependence of the N and A masses. The curves are two-parameter expressions for 
M]\j and according to the manifestly co variant xEFT framework of Ref. [105]. For both N and A, 
the chiral limit mass values M^^ (M^^) and the parameters cin (cia) are fitted. The red squares are 
lattice results from the MILC Collaboration [107] in full QCD. The stars represent the physical mass 
values. 

In Fig. 19, the -dependence of the magnetic 7*A^A-transition (Jones-Scadron) form fac- 
tor is shown in the relativistic xEFT framework of Refs. [88,89]. It is calculated from the 
one-loop diagrams of Fig. 7. Recall that the value of at = is determined by the low- 
energy constant gM- The -dependence then follows as a prediction of the NLO result, and 
Fig. 19 shows that this prediction is consistent with the experimental value at Q"^ = 0.127 GeV^ 
and physical pion mass. The m^-dependence of G\,j is also completely fixed at NLO, no new 
parameters appear. 

In Fig. 19, the result for 0%^ at = 0.127 GeV^ is shown both when the m^-dependence of 
the nucleon and A masses is included and when it is not. Accounting for the mTr-dependence 
in Mjv and Ma, as shown in Fig. 18, changes the result for G%f quite significantly. The xEFT 
calculation, with the m^r dependence of and Ma included, is in a qualitatively good agree- 
ment with the lattice data shown in the figure. The xEFT result also follows an approximately 
linear behavior in m^, although it falls about 10 - 15 % below the lattice data. This is just within 
the uncertainty of the NLO results. One should also keep in mind that the lattice simulations of 
Ref. [82] are not done in full QCD, but are "quenched", so discrepancies are not unexpected. 

In Fig. 20, the m^-dependence of the ratios Rem and Rsm is shown within the same rela- 
tivistic xEFT framework and compared to the lattice QCD calculations. As discussed above, 
the recent state-of-the-art lattice calculations of Rem and Rsm [82] use a linear, in the quark 
mass (nig oc m^), extrapolation to the physical point, thus assuming that the non-analytic rriq- 
dependencies are negligible. The thus obtained value for Rsm at the physical value displays 
a large discrepancy with the experimental result, as seen in Fig. 15. The relativistic xEFT calcu- 



37 



f\ I I I I I I I I I I I I I I I I I I I I I I I I I I I 

0.05 0.1 0.15 0.2 0.25 

m^'(GeV^ 

Fig. 19. The pion mass dependence of the real part of the Jones-Scadron 7*A^A form factor G\j for 
Q2 = and = o.l27 GeV^ in the relativistic ^EFT framework of Refs. [88,89]. The solid (dashed) 
curves are the NLO results for = 0.127 GeV^ (Q^ = 0) respectively, including the m,r dependence 
of Mn and A/a- The green dotted curve is the corresponding result for = 0.127 GeV^ where the 
dependence of Mj\[ and Ma is not included. The blue circle for = is a data point from 
MAMI [23], and the green square for = 0.127 GeV^ is a data point from BATES [95,80]. The three 
filled black diamonds at larger 771,7^ are lattice calculations [82] for values of 0.125, 0.137, and 0.144 
GeV^ respectively, whereas the open diamond near m^r ^ represents their extrapolation assuming 
linear dependence in m^. Figure from Ref. [89]. 

lation [88,89], on the other hand, shows that the non-analytic dependencies are not negligible. 
While at larger values of m^, where the A is stable, the ratios display a smooth dependence, 
at m^r = A there is an inflection point, and for < A the non-analytic effects are crucial, as 
was also observed for the A-resonance magnetic moment [108,109]. 

One also sees from Fig. 20 that, unlike the result for G\j, there is only little difference 
in the ratios between the xEFT calculations with the m^-dependence of Mn and Ma ac- 
counted for, and where this 171-,^ dependence of the masses is neglected. The xEFT framework 
of Refs. [88,89] thus shows that the assumption of a linear extrapolation in is not valid 
for Rem and Rsm- Once the non-analytic dependencies on the quark mass as they follow from 
XEFT are accounted for, there is no apparent discrepancy between the lattice results of Ref. [82] 
and the experimental results for Rsm- 

To test the difference between quenched and full lattice QCD results for the 7*A^A form 
factors, new lattice calculations, within full QCD, are underway [110]. The full QCD results 
obtained so far [110] using dynamical Wilson fermions are more noisy but in agreement with 
those obtained in the quenched theory. It will be interesting to test the predicted strong non- 
analytic effects in the Rem and Rsm ratios as shown in Fig. 20 once high statistics full lattice 
QCD results for pion masses smaller than 0.3 GeV become available. 
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Fig. 20. The pion mass dependence of Rem (upper panel) and Rsa4 (lower panel), at = 0.1 GeV^. 
The blue circle is a data point from MAMI [78], the green squares are data points from BATES [95,80]. 
The three filled black diamonds at larger are lattice calculations [82], whereas the open diamond 
near rriT^ ~ represents their extrapolation assuming linear dependence in m^. Red solid curves: NLO 
result when accounting for the dependence in and Ma; green dashed curves: NLO result of 
Refs. [88,89], where the ttItt -dependence of Mn and Ma was not accounted for. The error bands repre- 
sent the estimate of theoretical uncertainty for the NLO calculation. Figure from Ref. [89]. 

2. 7 Generalized parton distributions ( GPDs ) 

2.7.1 Definition ofN—>-N and N ^ A GPDs and sum rules 

So far we have discussed the ^ A transition as revealed with the help of the electro- 
magnetic probe. By measuring the response of the hadron to a virtual photon, one measures the 
matrix element of a well-defined quark-gluon operator (in this case the vector operator g7^g) 
over the hadronic state. This matrix element can be parametrized in terms of the 7*A^A transi- 
tion form factors, revealing the quark-gluon structure of the hadron. We are however not limited 
in nature to probes such as photons (or W, Z bosons for the axial transition). The phenomenon 
of asymptotic freedom of QCD, meaning that at short distances the interactions between quarks 
and gluons become weak, provides us with more sophisticated QCD operators to explore the 
structure of hadrons. Such operators can be accessed by selecting a small size configuration 
of quarks and gluons, provided by a hard reaction, such as deep inelastic scattering (DIS), or 
hard exclusive reactions such as deeply virtual Compton scattering (DVCS). We will be mostly 
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Fig. 21. The "handbag" diagram for the N ^ N and ^ A DVCS processes. Provided the vir- 
tuality of the initial photon (with momentum qh) is sufficiently large, the QCD factorization theorem 
allows to express the total amplitude as the convolution of a Compton process at the quark level and a 
non-perturbative amplitude parametrized in terms of generalized parton distributions (lower blob). The 
diagram with the photon lines crossed is also understood. 

interested here in DVCS reactions which are of the type 7* (qh) + N(p) — > 7(g') + B(p'), where 
the virtual photon momentum is the hard scale, and where the final state B stands for either 
the nucleon or the A state. The common important feature of such hard reactions is the pos- 
sibility to separate clearly the perturbative and nonperturbative stages of the interactions, this is 
the so-called factorization property. 

The all-order factorization theorem for the DVCS process on the nucleon has been proven 
in Refs. [111-113]. Qualitatively one can say that the hard reactions allow one to perform a 
"microsurgery" of a nucleon by removing in a controlled way a quark of one flavor and spin 
and implanting instead another quark (in general with a different flavor and spin) in the final 
baryon. It is illustrated in Fig. 21 for the case of the DVCS process with a A in the final state. 
The non-perturbative stage of such hard exclusive electroproduction processes is described by 
universal objects, so-called generalized parton distributions (GPDs) [1 14-1 16], see Refs. [1 17- 
120] for reviews and references. 

Before discussing the ^ A GPDs, we start by briefly reviewing the N ^ N GPDs. The 
nucleon structure information entering the nucleon DVCS process, can be parametrized at lead- 
ing twist-2 level, in terms of four (quark chirality conserving) GPDs . The GPDs depend on 
three variables: the quark longitudinal momentum fractions x and ^, and the momentum transfer 
= — to the hadron. The light-cone momentum ^ fraction x is defined by = xP+, where 
k is the quark loop momentum and P is the average nucleon momentum P = {p+p')/2, where 
p{p') are the initial (final) nucleon four-momenta respectively, see Fig. 21. The skewedness 
variable ^ is defined by g+ = —2^ P+, where q = p' — p is the overall momentum transfer in 
the process, and where 2^ —> xb/(1—xb/2) in the Bjorken limit: xb = Ql/{2p-qh) is the usual 
Bjorken scaling variable, with Ql = —q^ > the virtuality of the hard photon. Furthermore, 
the third variable entering the GPDs is given by the invariant = — g^, being the total squared 
momentum transfer to the nucleon. 

The DVCS process corresponds with the kinematics Q^, M^, so that at twist-2 level, 

terms proportional to Q"^ /Q\ or Mj^/Q^ are neglected in the amplitude. In a frame where the 



We do not consider chirally odd GPDs, which are also discussed in the literature, see e.g. the review 
of Ref. [119] 

^ We use the definition = (a° it a^) / for the light-cone components. 
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virtual photon momentum qf^ and the average nucleon momentum P^^ are collinear along the 
z-axis and in opposite directions, one can parametrize the non-perturbative object entering the 
nucleon DVCS process as (following Ji [1 15]) ^ : 



- I dy-e'^^'^y- {N{p'm-y/2)^-n^P{y/2)\N{p)) ^ ^ 
Ztt J y^=y±=o 

= H^{x, ^, Q') u{p) 7 • n u{p) + E%x, C, Q') u{p) iT'""^ u{p), (2.49) 

where ^ is the quark field of flavor q, u the nucleon spinor, and is a light-cone vector along 
the negative 2;-direction which can be expressed at twist-2 level ^° in terms of the external 
momenta as: 



= -2eP+^n, 5 = _2eP + g^, = (2.50) 

where q±_ — (0, q_L, 0) is the transverse component of the momentum transfer q, satisfying 
q± • n = qi_ ■ P = Furthermore the light-cone vector n, satisfying = is normalized in 
such a way that n- P = 1. The Ihs of Eq. (2.49) can be interpreted as a Fourier integral along the 
light-cone distance y^ of a quark-quark correlation function, representing the process where a 
quark is taken out of the initial nucleon (having momentum p) at the space-time point y/2, and is 
put back in the final nucleon (having momentum p ') at the space-time point —y/2. This process 
takes place at equal light-cone time (y"*" = 0) and at zero transverse separation (^j^ = 0) between 
the quarks. The resulting one-dimensional Fourier integral along the light-cone distance y~ is 
with respect to the quark light-cone momentum xP^. The rhs of Eq. (2.49) parametrizes this 
non-perturbative object in terms of the GPDs H'^ and E'^ for a quark of flavor q. The quark 
vector operator (7 • n) corresponds at the nucleon side to a vector transition (parametrized by 
the function H'^) and a tensor transition (parametrized by the function E'^). Analogously, there 
are two GPDs corresponding with a quark axial vector operator (7 • ny^), which are commonly 
denoted by the polarized GPDs H'^ and E'^. 

The variable x in the GPDs runs from —1 to 1. Therefore, the momentum fractions of the 
active quarks (x + for the initial quark and (x — for the final quark can either be positive 
or negative. Since positive (negative) momentum fractions correspond to quarks (antiquarks), it 
has been noted in [1 16] that in this way, one can identify two regions for the GPDs: when a; > ^ 
both partons represent quarks, whereas for x < both partons represent antiquarks. In these 
regions, the GPDs are the generalizations of the usual parton distributions from DIS. Actually, 
in the forward direction, the GPD H reduces to the quark (anti-quark) density distribution q{x) 
(q{x)) obtained from DIS: 



^ In all non-local expressions we always assume the gauge link: Pexp{ig J dx'^A^), ensuring the color 

gauge invariance. 

^°For kinematical expressions including correction terms proportional to Q'^/Qf^ and M^/Qf^, which 
are formally of higher twist, see e.g. Ref. [118]. 
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Hi{x,0,0) 



q(x), X > 0, 

(2.51) 

-q{—x), X < 0. 



The function E (and likewise the function E for the axial operator) are not measurable through 

DIS because the associated tensor in Eq. (2.49) vanishes in the forward limit (g ^ 0). Therefore, 
E is a new leading twist function, which is accessible by measuring hard exclusive electropro- 
duction reactions, such as DVCS. 

As the momentum fractions of initial and final quarks are different, one accesses quark mo- 
mentum correlations in the nucleon. Furthermore, in the region —^<x<^, one parton con- 
nected to the lower blob in Fig. 21 represents a quark and the other one an antiquark. In this 
region, the GPDs behave like a meson distribution amplitude and contain completely new in- 
formation about nucleon structure, because the region — ,^ < a; < ^ is absent in DIS, which 
corresponds to the limit ^ — > 0. 

Besides coinciding with the quark distributions at vanishing momentum transfer, the gener- 
alized parton distributions have interesting links with other nucleon structure quantities. The 
first moments of the GPDs are related to the elastic form factors of the nucleon through model 
independent sum rules. By integrating Eq. (2.49) over x, one obtains for any ^ the following 
relations for a particular quark flavor [115] : 



J^ldxH%x,C,Q') = F^iQ') , J^^'dxE%x,^,Q') = F^{Q') , (2.52) 

where (Q^) represents the elastic Dirac form factor for the quark flavor q in the nucleon. 
These quark form factors are expressed, using SU{2) isospin, as flavor combinations of the 
proton and neutron elastic form factors as: 



F^ = 2 Ff + + F^ , F^^ = 2 Ff + Ff + F{ , (2.53) 

where Ff and F" are the proton and neutron electromagnetic form factors respectively, with 
Ff (0) = 1 and Ff (0) = 0. In Eq. (2.53) Ff is the strangeness form factor of the nucleon 
(which is neglected in the practical calculations below). Relations similar to Eq. (2.53) hold for 
the Pauli form factors Fl . At = 0, the normalizations of the Dirac form factors are given by: 
Ff (0) = 2 (Ff (0) = 1) so as to yield the normalization of 2 (1) for the u ((i)-quark distributions 
in the proton. The normalizations of the Pauli form factor at = are given by Fl (0) = Kg 
(for q = u, d), where Kd can be expressed in terms of the proton (Kp) and neutron («„) 
anomalous magnetic moments as: 

Ku = 2Kp + Kn — +1.673, Kd = K,p + 2Kn — —2.033. (2.54) 

The sum rules of Eq. (2.52) also satisfy the condition that they are independent of ^, which is a 
consequence of Lorentz invariance . 



This is the simplest example of a so-called polynomiality condition when calculating moments of 
GPDs. 
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We next discuss the — ^ A matrix element for the vector twist-2 operator, which was worked 
out in Refs. [121,118]. There are four ^ A helicity amplitudes for the vector operator (as 
well as four for the axial-vector operator). However, the electromagnetic gauge invariance leads 
to only three electromagnetic form factors. Hence, there are four GPDs for the vector — > A 
transition, of which one has a vanishing first moment. In the following we shall neglect the 
GPD which has a vanishing first moment. This approximation is justified in the large Nc limit 
(discussed below) where this GPD is subdominant. The non-perturbative object entering the 
iV — > A DVCS process (lower blob in Fig. 21) can then be expressed as [1 18,121]: 



=y±=o 



(2.55) 



where u"{pa) is the Rarita-Schwinger spinor for the A-field, r3/2 is the third isospin generator 



for quarks, and y 2/3 is the isospin factor for the p A+ transition. 

Furthermore, in Eq. (2.55), the covariants A^^j;^'*^ are the magnetic dipole, electric quadrupole, 
and Coulomb quadrupole Jones-Scadron covariants [13] 



- 6(MA + Mjv) pA n a pKA (r,... 



3(Ma + M^) _ 



Here P = (pa+p) /2, q = Pa -p, and = M^. In Eq. (2.55), the GPDs Hm, He, and He for 
the N ^ A vector transition are linked with the three N ^ A vector current (Jones-Scadron) 
transition form factors G^, G*e, and introduced in Eq. (2.3) through the sum rules: 



£ dx Hm,e,c{x, e, Q') - 2 G*j^^E,c{Q'') ' (2.57) 

where the factor 2 arises because the electromagnetic form factors are conventionally defined 
with isospin generator T3/2 in the current operator in contrast to the operator T3 in Eq. (2.55) 
adopted in Refs. [121,118] to define GPDs. 

The above sum rules allow us to make a prediction for the iV ^ A form factors provided 
we have a model for the ^ A GPDs. Such a model will be discussed in the following 
sections. Conversely, the existing precise experimental information on the A^ — A vector form 
factors provides a strong constraint on the A?^ — > A GPDs. As discussed above, the GPDs are 
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however much richer observables and provide us with quark distribution information in the A 
resonance. They can be accessed by the N ^ A DVCS process as discussed in Ref. [122]. First 
experiments which are sensitive to the — > A GPDs have recently been reported [123]. 



2. 7.2 Model for the magnetic dipole N ^ A GPD 

Here we will be guided by the large Nc limit, which allows to connect the — > A GPD Hm, 
to the N ^ N iso vector GPDs. For the magnetic N ^ A transition, it was shown by Frankfurt 
etal. [121] that, in the large- A^c limit, the relevant N ^ A GPD Hm can be expressed in terms 
of the nucleon isovector GPD E'^ — E'^: 



Hm{x, e, Q') = {e-{x, e, Q') - E\x, Q') } , (2.58) 

where 3.70. Within the large A^c approach used in Ref. [121], the value 0*^(0) 

is given by G*m{0) = kv/ a/2, which is about 20% smaller than the experimental number . 
In order to give more realistic estimates, we will therefore use in the following calculations the 
phenomenological value G\j{Q) ^ 3.02 of Eq. (2.24). All other (sub-dominant) GPDs for the 
vector N ^ A transition vanish at leading order in the l/Nc expansion, consistent with the 
large Nc limit for the 7A^A form factors discussed in Subsect. 2.4. 
Using the large Nc estimate of Eq. (2.58), the sum rule Eq. (2.57) for can be written as: 




G*Mm = / ^ dx{E-(x,^,Q') - E'^(x,^,Q') 

_ G*M jO) 

Kv 



F^iQ') - F^iQ') , (2.59) 



where — F2 is the isovector combination of the proton (p) - neutron (n) Pauli form factors. 
Because the sum rule of Eq. (2.59) is independent of ^, we only need to constrain the GPD E'^ 
for ^ = in order to evaluate G*j^. The sum rule (2.59) was used in Ref. [124], using a model 
[126] in which the Gaussian ansatz for GPDs is modified at large Q"^ by terms having a power- 
law behavior. Refs. [127,128] used parametrizations which are motivated from the expected 
Regge behavior of the GPDs at small x and Q^. Guidal et al. [128] parametrized the function 
E'i{x, 0, (5^) at low (5^ through a Regge-type form (denoted by model Rl) as: 



Note the typo in the formula for Hm in Ref. [121]. Due to a different choice of isospin factors for the 
vector and axial vector transitions chosen there, one should correct Eq. (7) in Ref. [121] to be Hm = 
y^(^" _ E'^) instead of Hm = -^1 V2(£^" - E"^). 

In the large- A^c limit, the isovector combination H^ — H'^ is suppressed, therefore one could as well 
give as estimate iip^^+ ^ G^(0) ~ (/v,p — ~ 3.32 (where the magnetic moments are expressed 
in nuclear magnetons whereas is dimensionless), corresponding with Eq. (2.37) where Ma — Mjv, 
in the large- A^c limit. This estimate is accurate at the 10 % level. 
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0, = e«(x) x'^'^\ (2.60) 

The forward magnetic densities e'^{x) - unlike the usual forward parton densities q{x) - are 
unfortunately not known from experiment at present. The simplest idea is to take them propor- 
tional to the valence up-quark {uv{x)) and down-quark {dy{x)) densities as: 



e^{x) = -^Uy{x) and e'^{x) = K,ddy{x) , (2.61) 
which satisfy the normalization constraint of Eq. (2.52) at = 0: 



Kg = y dxe^^ix). (2.62) 

where and are defined in Eq. (2.54). One thus sees that the functions e^(x) encode the 
quark distribution information giving rise to the nucleon anomalous magnetic moments. 

As shown in Ref. [128], the Regge model R\ fits the nucleon Dirac (Fi) and Pauli (F2) 
form factor data for small momentum transfers < 0.5 GeV^. However, at larger the R\ 
model gives too strong suppression, and it consequently falls considerably short of the data for 
g2 > IGeV^. 

Experimentally, the proton helicity flip form factor F2(Q^) has a faster power fall-off at large 
than Fi{Q^). This means that the x ~ 1 behavior of the functions e(x) and q{x) should be 
different. To produce a faster decrease with Q^, the a; ~ 1 limit of the density e'^{x) should have 
extra powers of 1 — x compared to that of q{x). Aiming to avoid introducing too many free 
parameters, Guidal et al. [128] tried the next simplest ansatz for e^(x) by just multiplying the 
valence quark distributions by an additional factor (1 — x)'^", i.e. by taking: 

e«(x) = -^(1 - xy-Uy{x) and e\x) = -^(1 - xf-^d^ix) , (2.63) 
where the normalization factors Nu and Nd 



Nu ^ dx {I - xf- Uy{x) , Nd = dx {1 - xy^ dy{x) , (2.64) 

Jo Jo 

guarantee the condition of Eq. (2.62). In such modified Regge parametrization (denoted by R2), 
the GPD E'' entering the sum rule Eq. (2.59) for was parametrized in Ref. [128] as: 

EU^AQ') = ^{l-xr^qy(x)x^'^'-^^^\ (2.65) 

with q^ the valence quark distribution (q = u, d), and where the powers 1]^ and r]d have been 
determined from a fit to the nucleon form factor data as a = 1.105 GeV^^, = 1.713 and 
7]d = 0.566. Note that a value r]q = 2 corresponds to a l/Q^ asymptotic behavior of the ratio 
{Q')/F?{Q') at large Q\ 
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Fig. 22. The ^ A magnetic transition (Jones-Scadron) form factor G\i, relative to the dipole form 
(multiplied by a factor 3). The curves are calculated from a parametrization of the — > A magnetic 
GPD Hm of Ref. [128]. Blue dotted curve : Regge GPD parametrization according to Eq. (2.60). Red 
solid curve : modified Regge GPD parametrization according to Eq. (2.65). The data points are from 
JLab Hall C [130] (blue squares), and JLab CLAS : Refs. [81] (black triangles) and [131] (red circles). 
The data of Refs. [81,131] have been analyzed using the unitary isobar model of Ref. [132]. 



In the following estimates, the unpolarized valence quark distributions are taken at input scale 
/i^ = 1 GeV^ from the MRST2002 global NNLO fit [129] as: 



For the one-parameter Regge form Rl of Eq. (2.60) the same parameter value a = 1.105 
GeV^^ was chosen as in the model R2, which gives a good description of the proton charge 
radius. In Fig. 22, we show the results for the sum rule predictions for G\i using the GPD 
parametrizations of Eqs. (2.60) and (2.65). It is seen that both the Regge and modified Regge 
parametrizations yield a magnetic ^ A form factor which decreases faster than a dipole, 
in qualitative agreement with the data. The Rl Regge parametrization though gives too large 
suppression at larger Q^, as was also observed for the nucleon elastic form factors in Ref. [128]. 
For the modified Regge parametrization R2, it is seen that the sum rule prediction based on the 
large A^^^ estimate of Eq. (2.59) gives a good quantitative description of the data over the whole 
range without adjusting any parameters beyond the three parameters a' , and 7]^ which 
were determined from a fit to the nucleon elastic form factor data. 



= 0.262 x'^-^\l - (l + 3.83 + 37.65 x) , 
dy = 0.061 - (l + 49.05 + 8.65 x) . 



(2.66) 



(2.67) 
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2.7.5 Model for the electric quadrupole — > A GPD 

In this section, we make a very first attempt to model the electric quadrupole GPD, He- As 
in our modeling of Hm, we will also be guided by the large Nc limit. In the large Nc limit, 
Eq. (2.41) provides a relation between the N ^ A quadrupole moment and the neutron charge 
radius r^, which for Nc — 3 reduces to: 



Qp-.A+ = rl (2.68) 

Using Eq. (2.17), we can express Eq. (2.68) in a relation for G%{Qi), which reads to leading 
accuracy in the 1/A^c-expansion as: 



G,(0)^-i.J-L(<^. (,69, 

For small values of Q^, we can express the neutron electric form factor as G^iQ'^) ~ — 
Therefore, a natural extension of the large Nc relation of Eq. (2.69) to finite Q'^ is given by: 



Ge{Q ) - ^ ^ (2.70) 

The prediction which follows from the large A^,. motived expression of Eq. (2.70) is tested 
in Fig. 23 by comparing the Q'^ dependence of the neutron electric form factor G^ and the 
N ^ A Rem ratio. For Ge we use the modified Regge parametrization R2 of [128], which is 
seen to give a fairly good description of the available double polarization data. The Rem ratio 
is calculated using the large Nc relations Eq. (2.70) for G*e, and Eq. (2.59) for G^, as discussed 
in Sect. 2.7.2. These relations express the ^ A form factors in terms of nucleon elastic 
form factors. By using the three parameter R2 Regge form for the nucleon elastic form factors, 
we obtain in this way a prediction for Rem without adjusting any additional parameter. One 
sees that this yields a Rem ratio which has both the right size and displays a relatively flat 
behavior, up to a Q'^ value of about 1.5 GeV^, in surprisingly good agreement with the data. 

The form factor G*e is obtained from the first moment of the electric quadrupole ^ A 
GPD He through the sum rule of Eq. (2.57). We can therefore use Eq. (2.70) to propose a 
relation between the A^ — > A GPD He and the neutron electric GPD combination, which is 
consistent with this form factor sum rule, as: 



where the neutron GPDs are obtained in terms of the u- and d-quark flavor GPDs as: H^"'^ — 
-l/3i/" + 2/3H'^, and E'^''^ = -1/3E" + 2/3E'^. 
An interesting topic for future work is to perform different model calculations for Hm, He, 
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Fig. 23. GPD calculations for the dependence of the neutron electric form factor G% (upper panel) 
in comparison with the — > A Rem ratio (lower panel). For the neutron electric form factor, the 
modified Regge GPD parametrization R2 of Ref. [128] is used. For the Rem ratio, the large Nc relations 
Eq. (2.70) for G*^, and Eq. (2.59) for G\j are used, which express the ^ A form factors in terms of 
nucleon elastic form factors, which in turn are calculated using the modified Regge GPD parametrization 
R2. The data points for 0% are from MAMI [133-136] (red circles), NIKHEF [137] (blue square), and 
JLab [138-140] (black triangles). The data points for Rem are from BATES at = o.l27 [80] (blue 
square); MAMI (red circles): = q [23], = o.06 [79], = o.2 [141]; JLab CLAS [81] (black 
triangles); and JLab Hall A [142] (blue star). 

and He, as well as provide lattice QCD predictions for its moments, in order to cross-check the 
above estimates for Hm and He- 



2. 7.4 GPDs and transverse structure ofhadrons 

The interplay between the x and Q^-dependence of the GPDs contains new nucleon structure 
information beyond the information encoded in forward parton distributions depending only on 
X, or form factors depending only on Q"^. It has been shown that by a Fourier transform of the 
-dependence of GPDs, it is conceivable to access the distributions of parton in the transverse 
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Fig. 24. The magnetic dipole — > A GPD Hm as function of the quark momentum fraction x and 
the quark position b in the transverse plane (perpendicular to the average direction of the fast moving 
baryons), where b stands for either {h±)x or {h±)y. The calculation is based on the modified Regge 
parametrization R2 of Eq. (2.65). 

plane, see Ref. [143], and to provide a 3-dimensional picture of the nucleon [144]. 

For ^ = 0, one can define the impact parameter versions of GPDs which are obtained through 
a Fourier integral in transverse momentum q±. For the GPD E'^, this reads as: 



E'^ix, b^) ^ 1 1^ e^b^.q, _q2 )^ (2.72) 

and analogous definitions for the other GPDs. These impact parameter GPDs have the physical 
meaning of measuring the probability to find a quark which carries longitudinal momentum 
fraction x at a transverse position (relative to the transverse center-of-momentum) in a nu- 
cleon, seeRefs. [143,145]. 

When translating the GPD parametrization R2 of Eq. (2.65), into the impact parameter space, 
we obtain for the GPD E: 



^ g -bxVh4Q'(l-x)lnx] 



E'^ix, b^) = (1 - xr^ . . r-^. (2.73) 

Nq An [—a (1 — x) mx\ 

Using the large relation Eq. (2.58) we can then express the impact parameter version of the 
magnetic GPD Hm as: 
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Fig. 25. Perturbative QCD picture for the nucleon elastic and — > A transition electromagnetic form 
factors. The highly virtual photon resolves the leading three-quark Fock states of N and A, described by 
a distribution amplitude. The large momentum is transferred between the quarks through two successive 
gluon exchanges (only one of several possible lowest-order diagrams is shown). 

Hm{x, bx) = 2^mW fE^x, bx) - E'^ix, bx)| . (2.74) 



In Fig. 24, we display the impact parameter magnetic GPD HMix,h±) using the Regge 
parametrization R2 of Eq. (2.73) for E'^{x,h±). It is clearly seen from this image that for 
large values of x, the quark distributions are concentrated at small values of b^, reflecting the 
distribution of valence quarks in the core of the N and A. On the other hand, at small values of 
x, the distribution in transverse position extends much further out. 



2. 8 Perturbative QCD (pQCD ) 



2.8.1 pQCD predictions for helicity amplitudes and form factors 

The electro-excitation of the A provides a famous test for perturbative QCD, where scaling 
and selection rules for dominant helicity amplitudes were derived and are expected to be valid 
at sufficiently high momentum transfers [146]. A photon of sufficient high virtuality will see 
a nucleon (or A) consisting of three massless quarks moving coUinear with the nucleon. When 
measuring an elastic nucleon form factor or a ^ A transition form factor, the final state 
consists again of three massless coUinear quarks. In order for this (unlikely process) to happen, 
the large momentum of the virtual photon has to be transferred among the three quarks through 
two hard gluon exchanges as illustrated in Fig. 25. This hard scattering mechanism is gener- 
ated by valence quark configurations with small transverse size and finite light cone momentum 
fractions of the total hadron momentum carried by each valence quark. The hard amplitude can 
be written in a factorized form [147-149,146], as a product of a perturbatively calculable hard 
scattering amplitude and two distribution amplitudes describing how the large longitudinal mo- 
mentum of the initial and final hadrons is shared between their constituents. Because each gluon 
in such hard scattering process carries a virtuality proportional to Q^, this leads to the pQCD 
prediction that the helicity conserving nucleon Dirac form factor Fi should fall as l/Q^ (mod- 
ulo hiQ^ factors) at sufficiently high Q^. Processes such as in Fig. 25, where the interactions 
among the quarks proceed via gluon or photon exchange, both of which are vector interactions. 
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Fig. 26. 'y*N ^ A helicity amplitudes in the Breit frame, where N, A, and 7* momenta are collinear 
with the N momentum, and where the N and A tliree-momenta have the same magnitude but oppo- 
site sign. The 7, N, and A helicities (A^, Xn, and A a respectively) are indicated on the figure. The 
corresponding helicity form factors in the notation of Ref. [150] are indicated on the right. 

conserve the quark helicity in the limit when the quark masses or off-shell effects can be ne- 
glected. In contrast to the helicity conserving form factor Fi, the nucleon Pauli form factor F2 
involves a helicity flip between the initial and final nucleons. Hence it requires one helicity flip 
at the quark level, which is suppressed at large Q^. Therefore, for collinear quarks, i.e. moving 
in a light-cone wave function state with orbital angular momentum projection 1^ = (along the 
direction of the fast moving hadron), the asymptotic prediction for F2 leads to a l/Q^ fall-off 
at high Q"^. 

To derive the analogous relations for the 7*A^A transition form factors, it is most conve- 
nient to work in a Breit frame, in which the photon, A^, and A momenta are collinear with the 
incoming N. The outgoing A momentum has opposite direction but the same magnitude three- 
momentum. In this way, three independent N ^ A helicity form factors were introduced in 
Ref. [150] as illustrated in Fig. 26. They are defined as (for m = +, — , 0): 



G.. = ^ ( A, Aa = A, - - I e^T^ . J^^ \ N, = +-), (2.75) 

where is the electromagnetic current operator, and the factor l/(2M7v) is chosen to make 
Gm dimensionless. The 7, A^, and A helicities are denoted by A^, A^v, and Aa respectively. The 
transverse polarization vectors (for a photon moving in the 2-direction) are e''' = (0, =Fl5 —i, 0) / V2, 
and the longitudinal polarization vector satisfies e|°^ • e*^'^-''^ = 1, e^^^ ■ e^^^'' = 0, and ej^°^ -q^ = 0. 
One can express the N A helicity form factors G+, and Go in terms of the Jones- 
Scadron form factors 0%^, G% and as: 
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(2.76) 



Using the hadron helicity-conserving property of QCD at high Q'^, it is then easy to derive 
the asymptotic behavior of the 7*A'^A transition form factors. Since G+ is the only heUcity 
ampUtude with the same hehcity between initial N and final A, it will be the leading amplitude. 
At large , two-gluon exchange between coUinear quarks yields a 1 /Q^ behavior for [150]. 
The amplitudes Go and G_, requiring helicity flips, are asymptotically zero relative to G+. For 
each helicity flip, one expects one additional (m/Q) power suppression in pQCD, where m is 
some quark mass scale. In this way. Go (G_), requiring one (two) helicity flips, yield a large Q'^ 
behavior as i^/Q^) respectively [150]. 

Using the relations of Eq. (2.76), one can also derive the large behavior of the Jones 
Scadron form factors as well as the Rem and Rsm ratios. For the magnetic 'j*NA form factor 
G%[, one obtains al/Q^ behavior for — > oo. Furthermore, because G_ is asymptotically zero 
relative to G+, Eq. (2.76) yields G*^ — > — G^ in the limit Q'^ — > oo. This yields equivalently 
Rem +1 for — > oo. Likewise, the l/Q^ behavior for Go yields Rsm — ^ constant for 
Q^^oo [150]. 

We can test how well the above baryon form factor pQCD scaling predictions are satisfied 
at the available momentum transfers, see Fig. 27. One firstly sees from Fig. 27 that the proton 
Dirac form factor, which has been measured up to about 30 GeV^, displays an approximate 
scaling above 10 GeV^. For the proton ratio Ff/Ff , the data up to 5.6 GeV^ show no sign 
of a l/Q^ behavior as predicted by pQCD. Instead, the data show that the ratio / Ff falls less 
fast than with increasing Q^. In Ref. [156], the assumption of quarks moving coUinearly, 
underlying the pQCD prediction, has been investigated. It has been shown in Ref. [156] that 
by including components in the nucleon light-cone wave functions with quark orbital angular 
momentum projection = 1, one obtains the behavior F2/F1 ln^((5^/A^)/Q^ at large Q^, 
with A a non-perturbative mass scale. Choosing A in the range 0.2 — 0.4 GeV, Ref. [156] 
found that the data for Ff /Ff support such double-logarithmic enhancement. A same analysis, 
including states of orbital angular momentum projection = 1 in the or A light-cone wave 
functions, which find their physical origin in the transverse momentum of the quarks in the N 
and A, was also performed for the 7*A^A Coulomb transition form factor Gq [157]. It was 
shown in Ref. [157] that due to the orbital motion of the partons, Rsm acquires a double- 
logarithmic correction ln^(Q^/A^) at large compared with the standard scaling analysis, 
according to which Rsm constant for Q'^ — > 00. The arguments of Refs. [156,157] still rely 
on pQCD and it remains to be seen by forthcoming data at higher Q'^ if this prediction already 
sets in in the few GeV^ region. 

The test of the scaling behavior of the magnetic 7* A A form factor G^ is also shown on 
Fig. 27. As pQCD predicts a asymptotic scaling behavior, the ratio G^ / Go (with G o the 
dipole form factor) should approach a constant. The data for G^ up to about 6 GeV^ again do 
not support this scaling behavior. One sees instead from Fig. 27 that the data for G^ seem to 
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Fig. 27. Test of baryon form factor scaling predictions. Top panel : proton Dirac form factor multiplied by 
Q^; middle panel : ratio of Pauli to Dirac proton form factors multiplied by Q^; lower panel : ratio of the 
— > A magnetic transition (Jones-Scadron) form factor G\[ multiplied by to the dipole form factor 
Gd- The data for Ff are from [151] (solid squares). Data for the ratio F2 /Ff are from [152,153] (solid 
circles), [154] (open triangles), and [155] (solid triangles). Data for G\j- as in Fig. 22. The red curve is 
the calculation based on the three parameter modified Regge GPD parametrization R2 of Ref. [128]. 

support a similar behavior as for F2. The pQCD prediction for the Rem ratio which should 
approach +100% at high also fails dramatically. The Rem ratio is measured to be minus a 
few percent out to — 6 GeV^ [131], with no clear indication of a zero crossing. Also the 
RsM ratio which has also been measured up to — 6 GeV^ [131], does not seem to settle to 
a constant at large Q^. 



2.8.2 The road to " asymptopia" 

Although at high enough Q^, the pQCD scaling predictions should set in, the available data 
for the nucleon and ^ A electromagnetic form factors show that one is still far away from 
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this regime. This has been further investigated in several theoretical approaches. 

In Refs.[158,159], it has been argued that the above described hard scattering mechanism is 
suppressed at accessible momentum transfers relative to the Feynman mechanism [160], also 
called soft mechanism. The soft mechanism involves only one active quark, and the form fac- 
tor is obtained as an overlap of initial and final hadron wave functions. The hard scattering 
mechanism on the other hand, involving three active quarks, requires the exchange of two glu- 
ons each of which brings in a suppression factor «s/7r ^ 0.1. One therefore expects the hard 
scattering mechanism for Ff or to be numerically suppressed by a factor 1/100 compared 
to the soft term. Even though the soft mechanism is suppressed asymptotically by a power of 
relative to the hard scattering mechanism, it may well dominate at accessible values of 
Q^. In Refs.[158,159], the soft contribution to the 7*iVA form factors has been estimated us- 
ing a model based on local quark-hadron duality. In this approach it was found that the 7A^A 
transition is dominated by the magnetic form factor G^, while the electric quadrupole G% and 
Coulomb quadrupole form factors are small at accessible momentum transfers, in qualita- 
tive agreement with the data. 

In a more recent work [161], the soft contribution to the 7*A^A form factors was evaluated 
within the light-cone sum rule approach. In this approach, one analyzes a matrix element, in 
which the A is represented by an interpolating field rj^ of Ref. [162]. More specifically, one 
computes the correlation function of this interpolating field and the electromagnetic current 
operator given by the matrix element: 

T,u{P: q)=ij Ae^'-"(0 I T{r^^{Q)JM} I N{P)), {1.11) 

between the vacuum and a single-nucleon state \N[P)). It was found in Ref. [161] that the sum 
rule for G\j is dominated by contributions of subleading twist-4. They involve quark config- 
urations with a minus light-cone projection of one of the quark field operators, which can be 
interpreted as the importance of orbital angular momentum. The calculations of Ref. [161] are 
in agreement with the experimental observations that the Rem and Rsm ratios are small, al- 
though in the region of low < 2 GeV^, the result for is a factor of two below the data. 

In Sect. 2.7, we have shown that the nucleon elastic and A'^ ^ A transition form factors can 
be obtained from model independent GPD sum rules. These GPDs, represented by the lower 
blob in Fig. 21, are non-perturbative objects which include higher Fock components in the N 
and A wave functions. One can use a GPD parametrization to provide an estimate of the soft 
contributions, and expects this non-perturbative approach to be relevant in the low and interme- 
diate (5^ region for the form factors. This is shown in Fig. 27 (solid curves) from which one sees 
that the GPD Regge parametrization /?2, discussed in Sect. 2.7.2 is able to explain at the same 
time an approximate l/Q^ behavior for Ff and a behavior for -Ff/Ff which falls less steep 
than l/Q^. For G^, one sees from Fig. 27 that the GPD sum rule evaluation based on the large 
Nc relation of Eq. (2.59) is supported by the available data up to about 6 GeV^. Forthcoming 
experiments at the Jefferson Lab 12 GeV facility will extend the data for Ff/Ff and G^ to 
values around 15 GeV^. Such measurements will allow to quantify in detail the higher Fock 
components in the and A wave functions versus the simple three-quark Fock component, and 
pave the road to "asymptopia". 
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3 Phenomenology of pion photo- and electroproduction 

We will first review the dynamical model calculations which have been carried out by various 
groups and then turn to MAID, which is a variation of the effective Lagrangian method and has 
been very successful in describing the data. Lastly, we give a brief description on the current 
status of the dispersion-relation approach. 



3.1 Introductory remarks 

The electromagnetic production of pions is the main source of information about the elec- 
tromagnetic properties of nucleon resonances, such as the A(1232). Certainly the experimental 
measurements alone are not enough to obtain a quantitative insight, a theoretical interpretation 
of the resonance-excitation mechanism is necessary to extract, e.g., the strength of the 7A^A 
transition. One can distinguish three major theoretical tools which were developed in the last 
millennium for this purpose: 

(i) dispersion theory, first proposed by Chew et al. [163] for photoproduction, and by Fubini 
et al. [164] for electroproduction, has been successfully applied in the first multipole analy- 
ses [165,166] as well as in the modem ones [167-169]. The dispersion approach is based on 
general principles, such unitarity, analyticity, crossing symmetry, as well as relies on a phe- 
nomenological input from the nN scattering. A more detailed description of this approach is 
given in Subsect. 3.4. 

(ii) effective Lagrangian approach, where both nN scattering and pion production are calculated 
based on the same effective Lagrangian in terms of hadron fields. First attempts to develop 
this picture were made by Peccei [170], Olsson and Osypowsky [171]. More recently, this 
approach was extensively developed at RPI [172], Madrid [173], Gent [174], KVI [175], and 
Giessen [176]. It is important to emphasize that in this approach the effective Lagrangian 
is used only at tree level (nowadays usually by a unitarization procedure). This can be con- 
sidered as a weakness of this approach since, as is now known from xEFT calculations, the 
chiral loop corrections give rise to interesting and appreciable effects. 

(iii) dynamical models, where the tree-level effective Lagrangian is treated as a hadron-exchange 
potential of a quantum-mechanical scattering problem. First such models were developed 
by Tanabe and Ohta [177] and Yang [178]. Many more were developed over the past two 
decades, see it e.g., [179,180,14,181,91,182,183]. Three of these will be detailed below, 
where also the pros and cons of this approach in general will be addressed. 

In addition, the phenomenological multipole solutions SAID [184] and MAID [185] have 

proven to be useful in interpreting the experiment, as well as in obtaining the empirical infor- 
mation about individual amplitudes. Without these tools, it is extreamly difficult to extract the 
various amplitudes from experiment. In the rest of this subsection we discuss how it can be done 
in principle. 
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3.1.1 Measurement of pion photoproduction amplitudes 

Naively, without considering the discrete ambiguities, one would conclude that seven mea- 
surements are needed to determine the four pion photoproduction helicity amplitudes of H[s, [i = 
1, 4) (four magnitudes plus three phases) up to an arbitrary overall phase. The discrete ambigu- 
ities arise because the observables are all bilinear product of helicity amplitudes. To determine 
the amplitudes, mathematically speaking, amounts to solve a set of nonlinear equations. The 
solutions are not necessarily unique and the ambiguities could arise if the number of observ- 
ables measured are not enough. According to Ref. [186], nine measurements are required if all 
amplitudes would be determined without discrete ambiguities. Careful analysis by Chiang and 
Tabakin [187], however, reveals that only eight measurements would be sufficient to resolve 
all ambiguities, namely, four appropriately chosen double-spin observables, along with three 
single-spin observables and the unpolarized differential cross section. 

In pion photoproduction, there are 15 polarization observables (for detail see [10]). Among 
them three are three single polarization observables, namely, the polarized photon asymmetry 
S(^), polarized target asymmetry T{6) and the recoil nucleon polarization P{9). The rest are 
the 12 double polarization observables which describe reactions with polarized beam-polarized 
target, polarized beam-recoil nucleon polarization, and polarized target-recoil nucleon polar- 
ization. They are commonly denoted as E, F, G, H, Cx' , Cz' , Ox' , Oz' , Tx' , Tz' , Lx' , and Lz' . 
For definition see, e.g., Refs. [186,188]. 

We illustrate here only the differential cross section and one of the single polarization observ- 
ables, the polarized photon asymmetry, which has been extensively investigated experimentally 
in the study of E2/M1 mixing ratio for the N ^ A transition. 

The unpolarized differential cross section for the pion photoproduction can be expressed in 
terms of the helicity amplitudes Hi's (see Appendix) as follows: 



^ = ili V I H p 

The polarized photon asymmetry is defined as. 



^ = (3.2) 

where da± (da^) is the differential cross section for a linearly polarized photon with polariza- 
tion vector perpendicular (parallel) to the reaction plane. In terms of the helicity amplitudes, 
da± {dan ) given as 



da I 



■(I Hi + H^ |2 + I H2-Hs |2), 



dn 2|q| 



(3.3) 



56 



The expressions given above contain contributions from all partial waves. However, in the 
case of 7r° photoproduction in the A (1232) resonance region where most of the current exper- 
imental information are derived from, the dominant contributions come from s- and p-wave 
multipoles. Therefore, a truncated multipole approximation is often used in the analysis of the 
data for this reaction and only Eq^, Ei+, Mi+ and Mi_ multipoles are kept. Such an approxi- 
mation greatly simplifies the analysis and allows one to express the differential cross section in 
the following simple form : 



^ = - [Aj + B. cos 9 + cos' 9] (3.4) 
dVt q ' ' 

where the coefficients Aj^ Bj and Cj are bilinear functions of the s- and p- wave multipoles 
and J indicates the parallel (||), perpendicular (±), and unpolarized (0) components. These co- 
efficients are quadratic functions of the s- andp-wave amplitudes. The parametrization Eq. (3.4) 
can be used separately for da± and da\\ in the analysis of the photon asymmetry. In Refs. [22,24], 
it was found that da\i is very sensitive to the small Ei+ multipole due to the interference with 
large Mi+. We can see this from the following expressions : 

A|| = \Eo+\^+ I 3Ei+ + Mi+ - Mi_ p , 

S|| = 2Re [Eo+{m+ + Mi+ - Mi_)*] , (3.5) 
C\i^l2Re [£;i+(Mi+ - Mi_)*] . 

In Ref. [22,24] this sensitivity was used for the experimental determination of the Rem ratio of 
Eq. (2.14), because at the A resonance position one has. 



^--l^X- (3.6) 



3.1.2 Measurement ofpion electroproduction amplitudes 

Consider now the pion electroproduction^^ which kinematically is illustrated in Fig. 28. In 
this diagram the four-momenta of the initial and final electron are defined as, — (cj, Qil), If — 
(e/, Q/l); while the momenta of the other particles are given as, 

• in the lab frame: q = {u^^, Ql), p = (Mjv, 0), p' = (p'o^, Pl), and k = (a;k^, k^,); 

• in the center-of-mass (cm.) of the t^N system: q — [uj^, q), k — (a;k, k), p — [E, — q), and 
P' = {E\ -k). 

More kinematical details can be found in the Appendix. 



Until Sect. 6 we assume that the electroproduction process proceeds via the one-photon exchange. In 
this case, electroproduction differs from photoproduction only in that the incoming photon is now virtual. 
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Fig. 28. Kinematics for pion the electroproduction on the nucleon. 

The 5-fold differential cross section, in the case of the unpolarized target and beam, can be 
written as: 

r^, (3.7) 



which defines the virtual photon cross section. 



da,, dan 



dfl^ dCl^ 



dai 

da 



da- 



TT 



COS 207r + ^2£(1 + e) cos 0^ 



dn^ 



+ h\2e{l 



where 



— daTu ^.^ 



27r2 ei 



l-e 



(3.8) 



(3.9) 



s,h,9e,q,L,ei, and e/ denote the degree of transverse polarization of the virtual photon, the 
helicity of the incoming electron, the angle between the incident and final electrons, the three- 
momentum of the virtual photon, energy of the initial and final electron in the lab frame, respec- 
tively, r is the flux of the virtual photon and K — {s — M'^)/2M]sf, often called the 'equivalent 
photon energy'. 

In Eq. (3.7), ef,Qf denote the energy and solid-angle of the scattered electron in the labora- 
tory, respectively, and is the tilt angle between the electron scattering plane and the reaction 
plane, as shown in Fig. 28. d^l^, is the pion solid-angle differential measured in the cm. frame 
of the final pion and nucleon. Since F can be interpreted as the number of virtual photons per 
electron scattered into dej and dQ/„ the quantity on the r.h.s. of Eq. (3.7), da^/dQ^^ then repre- 
sents the cm. differential cross section for pion photoproduction by virtual photons. 

Also let us note that Eq. (3.8) differs slightly from Eq. (20) of Ref. [10] often used in the 



and hence in additional to transverse polarizations can have the longitudinal or scalar ones. 
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literature, where the degree of longitudinal polarization of the virtual photon el = {Q^ 
appears in some of the coefficients in Eq. (3.8). The difference can be viewed as a different def- 
inition for the longitudinal cross sections which simply modifies their relations with the nuclear 
response functions [189] by a kinematical constant, as will be specified below. 

The first two terms in Eq. (3.8) are the transverse (T) and longitudinal (L) cross sections. 
They do not depend on the pion azimuthal angle 0. The fourth and fifth terms describe the 
transverse-longitudinal interferences (TL and TL'). They contain an explicit factor sin^^ [190] 
and therefore vanish along the axis of the momentum transfer. The same is true for the third 
term, a transverse-transverse interference (TT), proportional to sin^ G. 

In the pion-nucleon cm. frame, it is useful to express these cross sections in terms of five 
response functions which depend only on three independent variables, i.e. it! = R{Q^, W, ©). 
The corresponding expressions are the following 



d<^T |k| d(TTT |k| 

d(^L |k| daxL |k| Q darv |k| Q 

d\ 2^ qw uj^ a\ 1^^ qw ai 2^ qw oOq 

with qw = (W^ — M^)/2W the "photon equivalent cm. energy", and cjq = {VV^ — — 
M'^)/2W the virtual photon energy in the cm. frame. The five response functions may be 
expressed in terms of the six independent CGLN amplitudes Fi,...,Fq. Explicit expressions 
can be found in Ref. [10]. In terms of the helicity amplitudes we can get simple expressions for 
the response functions: 

Rl = {\H,\^ + \He\'), 
Rtt = Re [H^Hl — H4H*) , 

Rtl = ^ Re[H;{H, - H^) + H;{H2 + H^)], 

Rtu = Im [{H, - H^)H; + {H2 + Hs)H;] . (3.11) 

Many pion electroproduction experiments in the A resonance region have been performed to 
study the electromagnetic — A transition. Most of them are on 7r° production reaction, with 
either unpolarized electron beam {h = 0) or longitudinal electron beam (/i ± 1). By measuring 
coincident cross sections at different azimuthal angles (0 = 0, in plane and 7^ out of 
plane) and different polar angles 9, it has been possible to obtain data on arr, cflt, ctlt' and un- 
polarized cross section ao — aT + sgl- Values for the Mi+, Rem, and Rsm have been inferred 
from the models which provide an overall agreement with the data. This will be discussed in 
more detail in Sect. 5. 

We mention in passing that, again naively, without considering the discrete ambiguities, one 
would expect that eleven measurements would be required to determine the six pion electro- 
production helicity amplitudes of H[s, {i = 1, 6) (six magnitudes plus five phases) up to 
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an arbitrary overall phase. However, in contrast to the case of photoproduction, the question of 
how many measurements are needed to determine all amplitudes without discrete ambiguities 
has not been properly addressed yet. 



3.2 Dynamical models 



Here we turn to the discussion of the dynamical models and as an example consider three of 
the more recent ones: the model of Sato and Lee (SL) [14], DMT (Dubna-Mainz-Taipei) model 
[91], and the DUO (dynamical Utrecht-Ohio) model [183]. 



3.2.1 General framework: unitarity, relativity, gauge invariance 

The dynamical model formulation for the pion electromagnetic production reaction starts 
from the following scattering equation 



for the following four processes 




U t 



(3.12) 



TlN - 

7* TV 



ttN, 
> ttN, 



7*7V ^ 7*iV, 



(3.13) 



where the transition matrices fs are related to the elements of the T-matrix by a kinematical fac- 
tor and t>'s are the driving potentials of the vrA^ scattering (vrvr), pion electromagnetic absorption 
(717), production (77r), and the nucleon Compton scattering (77). and are, respectively, 
the pion-nucleon and photon-nucleon two-particle propagators. 

To first order in e, the transition matrix element for 7*A'^ — > nN, in a dynamical model is 
then given as 



with 



(3.14) 



t-irN — Vt^n + Vt^nGqIt^n, (3.15) 

where we have replaced the subscript tttt used in (3.12) by t:N and introduced Co = Ctt 
which will be used hereafter. In this approximation, only the integral equation Eq. (3.15) for 
TiN scattering has to be solved and the rest is determined in a one-loop calculation. 

In a Bethe-Salpeter (BS) formulation, Eqs. (3.14) and (3.15) are four-dimensional equations 
and the driving terms v's represent the sums of all two-particle irreducible amplitudes. Two 
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approximations are commonly made to simplify the solution of the BS equation [191]. The first 
is to approximate the driving potentials by the tree diagrams of an effective chiral Lagrangian 
and the other is to replace Go by a propagator go which would reduce the dimensionality of 
the integral equations from four to three, go is chosen to have an invariant form so that the 
covariance of the original equation is preserved. The reduction in the dimensionality is achieved 
with a 5— function in go which imposes a constraint on the time component of the relative 
momentum variable, i.e., the relative energy. In addition, this new propagator must be chosen 
such that the resulting scattering amplitude has a correct elastic cut from the elastic threshold 
(m^r + Mtv)^ to oo in the complex s plane, as required by the unitarity condition. 

It is well known for example, see Ref. [192], that the choice of go is rather arbitrary and can 
be done in infinitely many different ways. In Ref. [191], a class of three dimensional equations 
of the following form, 

Go{k;P) = -^^ I fis,s')[ais,s')r+l^ + M^] 

[Zn)'^ J s — s + IS 

X 5^+\[rjN{s')P' + A;]2 - m%)5^+\[ri^{s')P' - k]^ - ml). (3.16) 

were employed to investigate the ttN scattering. In Eq. (3.16), P' — ^J^P defines the "off- 
shellness" of the intermediate states. The superscript (+) associated with 5-functions means that 
only the positive energy part is kept in defining the nucleon propagator, k = 'r]T,{s)p — riN{s)q 
is the relative momentum k with r]'s any function of s constrained by the condition r]T,{s) + 
Vn{s) = 1. To have a correct TriV elastic cut, the arbitrary functions f{s,s') and a{s,s') must 
satisfy the conditions 

/(s,s) = l, a{s,s) = 7]n{s). (3.17) 

It is easy to verify that for (m^ + Mjv)^ < s < oo, Eqs. (3.16) and (3.17) give the correct 
discontinuity of the propagator Go 

Disc[Go(k; P)] = -^{7^^(s/) + ^ + MN)5^+\[r)N(s)P + k]^ - M^) 

x5^+\[r]^{s)P -kf -ml). (3.18) 

The class of propagators given in Eq. (3.16) has the feature that both particles in the interme- 
diate states are put equally off-mass- shell such that the relative energy dependence in the inter- 
action is removed. Several three dimensional formulations developed in the literature, including 
those developed by Blankenbecler and Sugar [193], Kadyshevsky [194], Thompson [195], and 
Cooper and Jennings [196], can be derived from using Eqs. (3.16) and (3.17) and were studied 
in detail in [197] for TriV scattering. All of these schemes set ^n{s) = 6^(8) /{eN{s) + eT^{s)) 
and ?7^(s) = e^{s)/{eN{s) + e^{s)), where eN{s) = {s + M% - m\)ll^s and e^{s) = 
(s — M"^ + m^)/2y^ are the center of mass (cm.) energies of nucleon and pion, respectively. 
The resultant ttN interaction obtained with Cooper- Jennings reduction scheme has been exten- 
sively used in the DMT model calculation for the pion electromagnetic production reactions. 
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The functions a{s, s') and f{s, s') specific to the Cooper-Jennings reduction are 



a(s,s')^r]N{s), 

f{s, s ) = - — -g— 77-^. 



(3.19) 
(3.20) 



The resulting three-dimensional scattering equations then take the Lipmann-Schwinger form 



t-irN — Vt^n + VnNgotnN- 



(3.21) 



One of the most important features of the dynamical model approach is that it provides a 
unified theoretical framework to describe ttN scattering and pion production in a consistent 
way. In most cases, v^^n is derived from a tree approximation to an effective chiral Lagrangian 
which involves field operators of the nucleon, A, pion, rho meson, and a fictitious scalar meson 
a [191]. Pseudovector ttNN coupling is used as it is consistent with the leading order of chiral 
perturbation theory. It leads to a driving term which includes the direct and crossed and A 
terms, and the t-channel a- and p-exchange terms. 

Furthermore in the dynamical models, one needs to specify the 7A^A vertex, shown in Fig. 1 . 
It is described by Eq. (2.2) with qm, Qe, and gc, the magnetic dipole Ml, electric quadrupole 
E2, and Coulomb quadrupole C2 excitation strength of the A, respectively. They are then only 
three new parameters in the dynamical model, besides those determined by the t:N scattering 
and — > 7r7 reactions, to be determined from the electromagnetic pion production data. 

Since multiple pion rescattering in the final state is treated explicitly, unitarity in inherent in 
Eq. (3.21). In fact, if we take 



1 



E-Ho + is 



(3.22) 



with Hq the free Hamiltonian of the pion-nucleon system, then a multipole decomposition of 
Eq. (3.21) gives [178] 



q;E + is) = e'^°' cos 5a 



X 



77r 



{kE,q) 



Jo E- E^N{k') 



(3.23) 



where 5a and R^j^ are the ttN scattering phase shift and reaction matrix, in channel a, re- 
spectively; kE is the pion on-shell momentum and q —\ q | is the photon momentum. In the 

energy region where only two channels are open: ttN elastic scattering and single pion photo- 
or electroproduction, v^^,, f^jv and the reaction matrix i?^]v are real numbers. In this case we 
see explicitly from Eq. (3.23) that the phase of the t"^ is equal to the ttN scattering phase in the 
corresponding channel, e.g.. 
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t^,{kE,q;E)^\t^^{kE,q;E)\e 



(3.24) 



This is the well known Fermi-Watson theorem [198]. This theorem, which is a consequence the 
unitarity of the S-matrix, and time-reversal invariance, imposes a phase condition on the indi- 
vidual multipole amplitudes below the two-pion production threshold. It simplifies the analysis 
of the pion photo- and electroproduction processes in the A resonance region connecting the 
real and imaginary parts of the reaction amplitudes. 

Gauge invariance of the electromagnetic interaction requires the following current conserva- 
tion condition. 



as ty^ is proportional to T/j = £,x J'^ of Eq. (A.3). The first term in t^^r of Eq. (3.21) clearly 
satisfies the current conservation condition of Eq. (3.25) for on-mass-shell incoming and outgo- 
ing particles since it is obtained from the tree diagrams of a Lagrangian resulting from gauging a 
chiral effective Lagrangian for ttN scattering. However, there is a fundamental difficulty to im- 
pose the current conservation condition on the second term of Eq. (3.21), which represents 
the sum of the following ladder series, 

VjTvgotTTN = v^T,gQVT,N + v^T,gQVT,Ngov-„N H ■ (3.26) 

It is well-known that the sum of a set of diagrams is gauge invariant if photon is hooked to 
every line which carries charge. Consequently, v^t^qoIt^n is not a gauge invariant quantity. It is 
inherent in the approximation scheme which leads to Eq. (3.21). There are many recipes pro- 
posed [199,179,91,183] to make it to satisfy the current conservation condition of Eq. (3.25). 
However, they are all ad hoc in nature and hence not unique. 

The dynamical model approach as summarized in Eq. (3.21) hence contains four theoretical 
ingredients. The first one is the choice of f^^. It is commonly chosen to consist of tree dia- 
grams of a chiral effective Lagrangian [171], which include the Bom terms in pseudoscalar 
coupling, contribution from channel (p, uj) vector meson exchanges, and s— and it— channel 
A— exchanges. Different dynamical calculations differ mostly in the pNN and ujNN coupling 
constants used. The second ingredient concerns the choice of the three-dimensional propagator 
Qq. Another input is the model chosen for t^^jv since experimental TriV phase shifts constrain only 
the on-shell behavior while the physical pion production multipole amplitude t^^^ depends on 
the half-off-shell matrix elements of i^jv. The off-shell t^N rescattering effects have been shown 
[200] to play an important role to explain the threshold 7r° data. The last theoretical ingredient 
is the recipe employed to satisfy the current conservation condition of Eq. (3.25). 

3.2.2 Sato-Lee model 

The Sato-Lee (SL) model [14] made a strong effort to derive the driving terms of v-^^^ and v^^^ 
of Eq. (3.21) in a consistent manner. 
SL started from a model Lagrangian with N,A,7r, and p fields which would generate the tree 



(3.25) 
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diagrams of ChPT. The strong interaction Lagrangian is then extended to include the 7 field 
with "minimal substitution". The Lagrangians which describe the 'JttV interaction of Eq. (3.43) 
with V = (p, uj) and the 7A^A vertex of Fig. 1 are then added. The unique feature of SL's 
calculation is that they then use a unitary transformation method, called the SKO method [201], 
to derive from the above-mentioned Lagrangian an energy-independent effective Hamiltonian. 
The essence of the SKO method is to systematically eliminate the virtual processes from the 
considered Hamiltonian by using unitary transformation. The so-called "virtual processes" are 
the processes like N Ntt, N ^ Np, N ttA, and tt vrp, which can not take place in the 
free space because of the energy-momentum conservation. The effects of the virtual processes 
are included as effective operators in the resultant Hamiltonian. Another advantage of such a 
scheme is that it does not have to perform the renormalization for the nucleon as the N <-> Nn 
vertex has been transformed away since it is a "virtual" process. 
The final effective energy-independent Hamiltonian of SL model then takes the form 

Heff ^Ho + v^r, + < + {h^"L + h% + h.c.), (3.27) 

where v^j^ is the background ttN potential, and v!^^ describes the background — > ttN 
transition. contains Bom terms and channel p exchange, while v^^ is consisted of Bom 

terms in PV coupling and channel (p, a;) vector-meson exchange. h^^liA ^^'^ ^^jNA denote the 
ttN — > A^°) and 7*A'^ — > A^°) excitations of a bare A^^\ respectively. The matrix elements of 
hf}^^ and /i^^ta take the familiar forms, cfr. Sect. 2. 



<A'\h%\ci> 



1 1 



^"W + ^'^-lS.kjT.. (3.28) 



^jv(q) + Mjv 3 (Ma + 

X {iGl,{Q^)S X q • e + G*e(Q^)(S • etr • q + S • qtr • e) 



Ma 



-S • qo- ■ qeo)}, (3.29) 



where T and S are the isospin and spin | | transition operators, respectively. 
The driving terms v's in Eq. (3.21) in the SL model are then given by 



VnN = V^TV + V^N, '^jn = v!^^ + V^^, (3.30) 



where 



/,(0)t , (0) ;,(0)t , (0) 

A _ 'Vwa/Vtva a _ ''-TrAfA'SA^A /"2 -an 
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The energy denominators (E—M^^) inEq. (3.31)arise because SL have chosen gro in Eq. (3.21) 

to be the Schroedinger propagator, i.e, go = 1 / (E — Hq). 

Following Ref. [177], SL decomposed, with the use of two-potential formula, the resulting 
t-matrix as follows [202], 

t,^{E) = i^,{E) + t^,{E), tME) = i^AE) + i^AE): (3.32) 

where the background contributions if^^{E) and i^j^{E) are given by. 



i^,{E) = + go{E) i^AE), i^AE) = v^^ + v^^ go{E) i%{E). (3.33) 

This kind of background terms are further denoted as "non-resonant" background because they 
do not contain any resonance contributions from v^^ or v^j^ of Eq. (3.31). Note that in partial 

channel a, the on-shell matrix elements of both i^^{k,E, q', E + ie) and t^AE){kE-, kE] E + ie) 
have the same phase 5^ as can be proven in the same way as Eq. (3.23). 
The A contribution to t^N scattering i^^{E) in Eq. (3.32) takes the form 



i^AE) = hij,AE)9AKNA{E), (3.34) 
where h.^NA{E) describes the dressed ttA^^A^ vertex 

Kna{E) = h^^L + h^'L 9o{E) i^AE), 

hi^AE) = /^S/a + i^AE) 9o{E) hfS^. (3.35) 

It can be easily seen that the physical matrix elements of both h-j^NAiE) and h^NAiE), i.e., 
Kna{E + is, kg) and hlj^AE + £, k^;), have the phase 5a. Note that hlj^AE) ^ ^!rivA(-^)- 
gA{E) is the dressed A propagator, 

g^^^g^^-E^iE), (3.36) 
where the A self-energy Ea {E) is given by 



EA{E)^h^;l^gohi^AE). (3.37) 

The matrix element of T^a{E) is related to the dressed mass Ma{E) and the width Ta{E) of 
the physical A by 



< AW \^a{E + ie) I AW >= MAE) - i^^. (3.38) 
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Similarly, the A contribution toY^ reaction t^T^{E) can be expressed as 

i%{E) = hi^^{E)g^h^M^{E), (3.39) 
where 

h,NA{E) = h% + /iS^ go{E) i^,{E) 

= h^'L + Kna{E) go{E) v^,. (3.40) 

To solve the integral equation for i^j^ in Eq. (3.33), SL introduced dipole form factors for 
i'^jv and the ttA^A vertex. The form factors they introduced for the different pieces in v!^]^ are 
different. The off-energy-shell matrix elements of v^^ are needed to calculate t^^ of Eq. (3.33) 
and hyNA{E) of Eq. (3.40). This is dictated by the unitary transformation they used as it re- 
quires that the time components time component of the momentum in the propagator of any 
— > ttN amplitude is evaluated by the external momenta associated with the strong interac- 
tion vertices. 

For the electroproduction, Sato and Lee used the following substitution for the current oper- 
ator 



JM ^ j/^ _ i—n'^, (3.41) 
n ■ q 

with n — (0, 0, 0, 1) for q — (a;, 0, 0, q) to preserve gauge invariance in their calculation. This 
amounts to adding interaction currents Jz = Jz{SL) + Jz{int) such that Jz = q^ /\q\Jo{SL), 
instead of Jz{SL), is used in the calculations of longitudinal cross sections. This simple pre- 
scription is identical to what has been commonly used in many nuclear physics calculations 
[203]. 

In short, Sato and Lee derived an effective Hamiltonian consisting of bare A nN, 
vertex interactions and energy-independent meson-exchange vriV, 7iV — > vriV transition poten- 
tial operators, by applying a unitary transformation to a model Lagrangian of iV, A, tt, p, cj, and 
7 fields in order to achieve consistency in describing 7*7r reactions and nN scatterings. Form 
factors are added to the resultant effective Hamiltonian and with the use of a Schroedinger 
propagator, the parameters are adjusted to give a good description of the existing data. 



3.2.3 Dubna-Mainz-Taipei model 

Dubna-Mainz-Taipei (DMT) model was developed in two stages. When the dynamical ap- 
proach was first proposed in Ref. [178], a phenomenological separable form for f^jv and the 
Schroedinger propagator {E — Hq)^^, together with a derived from a chiral effective La- 
grangian, were employed in Eq. (3.21). Efforts were then made to construct a meson-exchange 
model for the nN scattering [180,191] (called Taipei- Argonne nN model hereafter), in order to 
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achieve consistency in the treatment of ttN scattering and reactions. However, the resultant 
ttN interaction was used in Eq. (3.21) only to describe the pion threshold photoproduction in 
[180]. Putting together the meson-exchange models of Vj^^ and Vt^n and apply them to the photo- 
and electroproduction of pion in the resonance region was performed only recently in Ref. [91]. 

In the TA nN model, several three-dimensional reduction schemes of the Bethe-Salpeter 
equation for a model Lagrangian involving N, A, vr, a, and p fields were investigated. It was 
found that all of the resultant meson-exchange ttN models can give similar good description of 
the ttN phase shifts up to 400 MeV. However, they have significant differences in describing 
the ttNN and ttA'^A form factors and nN off-shell matrix elements. 

The TA model contains almost identical physics input as the SL nN model. Namely, the 
Lagrangian they considered both consist of tt, p, N, A fields except that the TA model also in- 
cludes the t— channel a exchange. The relevant interaction Lagrangian used in the TA model is 
given as follows: 



AO) (v) 

-g.NNNaN - g,NNN{j,p^ + ^^.Ad^'p" - d'^p'^)} ■ ^tTV 

(0) _ 1 

+{^A^[^^^ - (Z + -)YY]T^N ■ d^TT + h.c.}, (3.42) 

where A^ is the Rarita-Schwinger field operator. The parameters, namely, all the coupling con- 
stants, the bare nucleon and A masses and , and the "off-mass-shell coupling param- 
eters" Z for the ttA^A vertex [171], are adjusted such that the predictions of ^^at agree well 
with the experimental data, including nucleon mass, renormalized iiNN couplings constant 
fl^j^/An = 0.079 and the s— and p— wave ttN phase shifts. 

The effective chiral Lagrangian for the strong interaction is then extended to include the 
electromagnetic interaction in a gauge invariant way, i.e., with the minimal substitution of 
^ — ieAfj,. In addition, the channel vector-meson (p, (v) exchange is known to make 
non-negligible contribution. The relevant effective Lagrangian is given by 



777. 7r 

Unn = 9vnnN (^^^V^ - ^a^^d'^V^''^ N , (3.44) 

where N and V {= p,uj) denote the nucleon and vector meson field, respectively. Xy is the 
radiative coupling determined hy V ^ tt'j decay. The interaction described above would give 
rise to a contribution to of Eq. (3.21) which contains the Born terms in pseudoscalar cou- 
pling and the t— channel (p, u) vector-meson exchanges and is normally called the background 
transition potential and denoted by v^^. 
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However, different theoretical treatments lead to difference in f even for the same set of 
diagrams considered. For example, for the m— channel nucleon exchange diagram, TA model 
has (omitting the isospin indeces) 



k';p, k)=['-^) ^,hSN{p' - k) 75 H\ (3.45) 

where Sn{p) = {i> — + ie)~^ is the nucleon propagator. In SL model, however, with the 
use of the unitary transformation, Eq. (3.45) becomes, 

v:],'^{p',k';p,k) = {l^f^,/i[S^{p'-k) + SM{p-k'h,}^'h,}^', (3.46) 

m-^ 2 

instead. In addition, SL used go = 1 / {E — Hq) in solving the equation t^^r = v^^ + VT^Ngot-KN, 
while the TA model considered many different three-dimensional propagators. 

In the DMT model calculations of [91] for electromagnetic pion production, only the meson- 
exchange TiN model for v^.^, obtained with the Cooper- Jennings [196] three-dimensional re- 
duction scheme, has been employed. It was chosen because Cooper and Jennings [196] argued 
that the reduction scheme they developed has the advantage of preserving the chiral symme- 
try contained in the original Lagrangian. However, two remarks are in order here. First, even 
though the t^^N matrix elements used in the DMT calculations are obtained within the Cooper- 
Jennings' reduction scheme, namely, the time component of the relative momentum in Vt^n is 
fixed according to the Cooper- Jennings prescription, all the particles are put on mass-shell in 
the derivation of f of DMT. In addition, the Cooper- Jennings propagator was used in solving 
^TrAf = + VTrN9otTrN but thc Schrocdingcr propagator is used for go instead in the evaluation 
of tjjr = Vjjr + v^irgotTrN- With the differences in Vt^n and go as expounded above, it will not be 
a surprise if the off-shell behaviors of the i^jv obtained in SL and TA would differ substantially 
in some kinematical regions. 

The prescription used to preserve gauge invariance in the DMT model goes as follows. The 
aim is to continue k in v'^^{k, q), which is originally derived only for the on-energy-shell val- 
ues ks, to the off-energy-shell region as called for in Eq. (3.23). It starts from the expressions 
for the CGLN amplitudes F/s derived from the tree approximation to the effective chiral La- 
grangian used for the model. The F/s can be expressed as linear combination of the the invariant 
amplitudes A[s of Eq. (A. 6) as given in Eqs. (la-b) in the Appendix of Ref. [204], where the 
coefficients are functions of Q'^, the energies of photon, pion, initial and final nucleon, qo,ko,E, 
and E', respectively, and the total energy W. These energy variables are then replaced with 
expressions in three-momentum |q| and |k|, e.g., , = ^|qP — Q^, ko — rn^ -\- and E — 

\J Mjf + q2 etc, while W is always kept as the initial energy, namely, W = |q| -I- \J Mj^ + q^. 
Such substitutions have the property that the threshold behaviors of the F/s are retained, i.e., 

Fi oc const., F2 oc |q||k|, F3 oc |q||k|, 

F4 oc |k|^ Fr oc const., Fg oc |k|. (3.47) 
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In the invariant amplitudes, the Mandelstam invariants t and u appearing in the t- and li-channel 
exchange diagrams, are expressed in terms of the three-momentum |q| and |k| by : 



t^(p- = 2M^ - 2^M2 +q2^M2 +k2 + 2q • k, (3.48) 



u={p-ky = M^ + mi- 2^jMl + q^m^ + - 2q ■ k. (3.49) 

The resulting expressions of F/s, which are functions of |q|, |k| and x = cosO are then used to 
obtain, according to Eq. (If) in Ref. [204], the multipole amplitudes v^^{k, q) where k can now 
be off-energy-shell. A dipole form factor [(«^ + k|;) / (a^ + k^)]^ with a = 440 MeV is then 
multiplied to all the resultant multipoles to ensure the convergence in the integral in Eq. (3.23). 

There is another difference between SL and DMT dynamical model calculations. It lies in 
the scheme in the separation of background and resonance contributions. In the SL model, the 
full tjT, is decomposed as in Eqs. (3.32-3.40). In the DMT model, is instead decomposed as 
follows. 



t,,{E)^t^,{E)+t^,iE), (3.50) 
where 

t^^iE) = <^ + < goiE) tME) , (3.51) 

t^,{E) = + ^o(^) tME). (3.52) 

Here t^^, in contrast to i^^ of Eq. (3.33), includes the contributions not only from the nonreso- 
nant mechanisms but also some of the contributions from the A excitation as contained in the 
full tjrAfwhich lead to the renormalization of the vertex 7*A^A. However, all processes which 
start with the electromagnetic excitation of a bare resonance are summed up in t^^. One feature 
of such a decomposition is that each term in Eq. (3.50) would fulfill the condition imposed by 
Fermi-Watson theorem, i.e., their respective multipole amplitude in channel a would have the 
same nN phase shift 6a as in Eq. (3.24). 





+ 

B ,~A 

Fig. 29. Relation between dressed and bare 7A'^A vertices of Eq. (3.54). The shaded and solid circles 
represent the dressed 7A'^A and ttA'^A vertices, respectively. 

Multipole amplitude of t^^ can be evaluated in a straightforward manner as in Eq. (3.23) once 
tT^N is given. For the resonance contribution t^^, it is easy to see that it takes the form of. 



t^,{E) = hi^^{E)g^{E)h%. (3.53) 



69 



The relation between t^^ and i^^ of Eq. (3.34), the A resonance contribution with dressed and 
bare e.m. vertices as depicted in Fig. 29, is given by 

iU^) = C(^) + < 9o{E) i^AE) . (3.54) 

Note that gA{E), which depends on the physical mass Ma and total width Fa, appears in 
i^^{E),t^^ and so the resonance position and the total width can be extracted from nN 

scattering. The DMT model calculation takes advantage of this fact and writes the resonance 
multipole amplitude of Eq. (3.53) in the following form : 

where /^a is the usual Breit-Wigner factor describing the decay of a resonance with total width 
V i^iyV) and physical mass Ma. Namely, the well known Breit-Wigner form for the resonance 
contribution is assumed, in close analogy to the standard way of analysis of the experimental 
data as done in [92]. The expressions for the form factors /^Aj /ttA and total width Fa will be 
considered in the next subsection. The energy dependence in /^a(W^) is introduced to account 
for the effect of possible excitation of the A via inelastic channels like ttA. The phase 
in Eq. (3.55) is used to adjust the phase of the resonance contribution to be equal to the cor- 
responding ttA^ phase shift ^ttat in the (3,3) channel as required by the Fermi-Watson theorem. 
At = Ma = 1232 MeV = for any Q^. Therefore, this phase does not affect the 
dependence of the 7A^A vertex. 

In the SL and DMT models, the dependence for the 7A^A vertex is parametrized in the 
same way as in Ref. [92], i.e. 



K{Qi') = (0)— (1 + P. Q') e-^-"^' Gzp(g'), (3.56) 

with a = M,E,S referring to magnetic dipole, electric quadrupole, or Coulomb quadrupole 
transitions, where qa is defined in Eq. (2.6), and where qw is defined below Eq. (3.10). The 
form factor Gd{Q'^) = 1/(1 + Q^/O.Tl)^ is the usual nucleon dipole form factor. In the case of 
the magnetic transition the parameters P and 7 can be determined by fitting A^{Q^) to the data 
for the well known G%f form factor, which is obtained from the Mi+ pion electroproduction 
multipole in the (3,3) channel atW = Ma = 1232 MeV, as given by Eq. (2.13). 

We want to remind the reader that the definition of the A e.m amplitudes in different models 
can be different. As we have seen above in the SL model they can contain contributions from the 
excitation of the A resonance via the nonresonant mechanism and describe the dressed 7A^A 
vertex. In the DMT model such mechanism is included in the unitarized background t^^ and the 
electromagnetic vertices A-^(Q'^) describe the bare e.m. vertices. 

In general, in accordance with the considered above dynamical models, the form factor 
can be decomposed in three terms 

Gl,{Q^) = Gt'iQ'') + Gr"'='°"'(g') + G^i^'-'^'^iQ^) , (3.57) 
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Fig. 30. Real and imaginary parts of the and multipoles. The dashed (dotted) curves are 

the results obtained within the DMT model for t^^ including (excluding) the principal value integral 
contribution in Eq. (3.23). The solid curves are the full DMT results including also the bare A excitation. 
For the Ei^ multipole, the dashed and solid curves practically coincide due to the small value of the bare 
The open circles are the results from the Mainz dispersion relation analysis [?], whereas the solid 
circles are obtained from the VPI analysis [205]. 



where G^^p , Gtr^^"""^ and G^/' '"^'^' are the contributions from t^^, from the second term in rhs 
of Eq. (3.54) and from the t^^ Eq. (3.33), respectively. Therefore, different models, which have 
a different way for the separation of background and resonance contributions, the values for the 
parameters A^(0), (3a and 7^ in Eq. (3.56) can be different. The same arguments are also true 
for the form factors of the electric and Coulomb transitions. 

In Fig. 30, we show the resonant multipoles Mj^i^^^-* and E^jj,'^'^ as obtained in the DMT model. 
For MijI_^' , one sees a large effect of the pion off-shell rescattering (difference between dotted 
and dashed curves), which results from the principal value integral part of Eq. (3.23). The total 
pion rescattering (dashed curves) contributes for half of the M]^'],^^^ as seen in Fig. 30 for the 
DMT model, the remaining half originates from the bare 7A^A excitation. Furthermore, one 
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Fig. 31. (Color online) The s-, p-, and some of the d-wave vrA^ scattering phase-shifts in the DUO model 
(red solid curves). The data points are the results of the SAID single-energy solution SP06 [184]. 



sees that almost all of the E2 strength is generated by the vrA^ rescattering. 



3.2.4 Dynamical Utrecht-Ohio model 

The dynamical Utrecht-Ohio (DUO) model has been developed in Refs. [206,183]. The 
model is based on a rrN-'jN coupled-channel equation which when solved to the first order 
in the electromagnetic coupling e leads to the electroproduction amplitude, T^^y* = + 
T7r7rG7rKr7*, whcrc V^^* is an basic electroproduction potential, Gj^ is the pion-nucleon propa- 
gator and Ttt^ is the full vrA^ amplitude. Thus, the pion rescattering effects are included as the 
final state interaction. 

The ttN amplitude satisfies an integral equation on its own. The details on constructing tiN 
amplitude are presented in [206] . The corresponding fit of this model into the nN elastic scat- 
tering phase shifts is shown in Fig. 31. 

In Fig. 32 one can see the diagrammatic representation of the DUO model for the pion elec- 
troproduction. The model potential V-^-y includes the Bom term (using the pseudo-vector ttNN 
coupling), the t-channel exchange of p— and cu— mesons, and the A-isobar exchange. 

The nN final state interaction dresses the s-channel nucleon and resonance contributions, 
leading in particular to the mass, field and coupling constant renormalizations. Therefore, both 
N- and A-pole contributions in V^^y are included using the bare mass and coupling parameters 
obtained from the equation for the vrA^ amplitude. The renormalization conditions together with 
unitarity demand that the same propagators and wN vertices, including the cutoff functions, 
appear in both the nN and potentials. Thus, all these ingredients are fixed by the analysis 
of TcN scattering. 

On the other hand, the electromagnetic interaction is constrained by the electromagnetic 
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Fig. 32. The electromagnetic pion production amplitude in the DUO model. 



gauge invariance. At this point one is often concerned with the problem of how to introduce the 
electromagnetic form factors for nucleon and pion in a way consistent with gauge invariance. 
A common solution to this problem is to choose all of the electromagnetic form factors that go 
into the Bom term (i.e., nucleon, pion and axial form factors) to be the same. This prescription 
does enforce the current conservation, however the Ward-Takahashi (WT) identities cannot be 
satisfied in this way. Furthermore, it is clear that the requirement of gauge invariance should 
not be able to restrict the dependence of the electromagnetic interaction. Finally, it is known 
that these form factors need not to be and are not the same. 

The DUO model introduces the procedure where an arbitrary form factor F{Q'^) can be ac- 
commodated by the following replacement of the current: 



J'' ^ J"'{Q^) = J^ + [F{Q^) - 1] O^V, , (3.58) 

where O^'' = g^" — q^q" /q^, and q is the photon 4-momentum, = —q^. It is easy to see 
that the resulting current J'^ obeys exactly the same WT identities as J'^. Thus, as long as 
gauge-invariance is implemented at the real-photon point, the inclusion of the form factors via 
Eq. (3.58) will give the gauge-invariant current for ^ 0. For example, the bare NN^ and 
7r7r7 vertex functions and the KroU-Rudermann term are: 



^NN^y = er + e [Fi(Q^) - 1] O'^'^t. + F^iQ') ta^'q, , (3.59a) 

ZTri]\f 

r^,^ = e{k + k'r + e [F^iQ') - 1] O^^ik + k% , (3.59b) 

Jkr-^^ + [Fa{Q') - 1] O^'^^^} 75 . (3.59c) 

This procedure allows one to use the experimentally determined form factors in the Bom terms. 
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3.3 MAID - a phenomenological multipole solution 



MAID, developed by the Mainz group, is currently the most comprehensive unitary isobar 
model which includes contributions from all other known resonances up to second and third 
resonance regions, in addition to the A(1232). Its first version was published in 1999 [92] as 
commonly referred to as MAID98. Since then it is updated once every two years to accommo- 
date for the new data by readjusting the parameters and can be accessed to, together with DMT, 
at the website of Ref. [185]. 

The basic equations used for the MAID are similar to those of the DMT model [91]. It also 
starts from the equation t^T,{E) = v^^, + v^^^ go{E) t^j^{E) with go — {E — Hq)~^. For any 
resonant channel, the transition potential v^t,- consists of two terms 



v^^{E)^v^^ + v^^{E), (3.60) 

where v^^ is the background transition potential and v^^{E) corresponds to the contribution of 
the bare resonance excitation. The resulting i-matrix can be decomposed into two terms as in 
Eqs. (3.50-3.52), 



where 



t^^{E) = v^, + < go{E) tME), (3.62) 
t^,{E) = v^, + ^o(^) Un{E). (3.63) 

In MAID, the background potential y^^" is constructed in the same way as in the DMT model 
except now the Bom terms are calculated with a hybrid model (HM), i.e., an energy dependent 
mixing of pseudovector-pseudoscalar (PV-PS) ttNN coupling , 

tHM _ ^^m tPV I tPS (1. (^A\ 

where ko is the asymptotic pion momentum in the ttN cm. frame which depends only on W and 

(3/2) (3/2) 

is not an operator acting on the pion field. From the analysis of the Mi_ and £'0+ multipoles 
it was found that the most appropriate value for the mixing parameter is = 450 MeV [92]. 
Note that for the pion pole term, the pion couples with on- shell nucleons only, and HM, PV, 
and PS couplings are all equivalent. As a matter of fact, only the multipoles Eq^, Mi_, Lo+, 
and Li_ are affected by the use of such a mixed coupling scheme. In all, v^^"" depends on 7 
parameters: the PV-PS mixing parameter A^, 4 coupling constants and 2 cut-off parameters for 
the vector-meson exchange contributions. 

In MAID98, the background contribution was given as t:^;"(MAID98) = v!^f{W,Q^), a 
real and smooth function. The unitarization of the total amplitude was done by introducing an 
additional phase in the resonance contribution which ensures the phase of the total amplitude 
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(background plus resonance) be equal to the corresponding pion-nucleon scattering phase 6q. 
In the new extended version (MAID2005), all the the background contributions of s-, p-, d- and 
/-waves are complex numbers as prescribed in K-matrix approximation, namely, 



t^fiW, Q') = v^fiW, Q') [1 + K^{W)] , (3.65) 

where t^j^ — [r]aex.p(2i6a) — l]/2i is the pion-nucleon elastic scattering amplitude with the 
phase shift 5„ and the inelasticity parameter rja taken from the analysis of the VPI group (SAID 
program). Below the two-pion threshold where T]a = 1, Eq. (3.65) becomes t^^^iWjQ"^) = 
e^^°'cos6av!^^"{Q'^). Such a structure of the background appears naturally from Eq. (3.23) if 
the principal value integral contribution in the pion rescattering term is neglected. Part of such a 
"pion cloud" contribution (principal value integral) will be effectively included in the resonance 
sector with the use of dressed ^NR vertex as will be explained below. 

As in the DMT model, the resonance contributions Q"^) in MAID is parametrized as 

in Eq. (3.55), 

where A^{Q'^) should now be interpreted, in contrast to A^{Q'^) of Eq. (3.55), as corresponding 

to the excitation of a dressed resonance in order to partly account for the neglected principle 
value integral term of Eq. (3.23) in t!^fiW, Q^) of Eq. (3.65). f^NiW) is the usual Breit-Wigner 
factor describing the decay of the A^* resonance with total width Ttou partial ttA^- width F^tv and 
spin j, 

1 kw Mn r^jv 



- Ml 



_{2j + l)n\k\WRrlt 

The factor Cttat is ^3/2 and -1/^3 for the isospin 3/2 and isospin 1/2 multipoles respectively. 
The energy dependence of the partial width r^^N is given by 



B r ( l''^' ( X' + k% y Wn 




where X is a damping parameter, TtotiW) is the total width, kji is the pion cm. momentum at 
the resonance peak (W = Wr) and is the single-pion branching ratio. Expressions for the 
total width can be found in [92]. The W dependence of the -jNN* vertex is given by the form 
factor 

where the damping parameter X is the same as in Eq. (3.68)) and Qb, = qw at W = Wr. 
The parameter n is defined by the best fit of the experimental multipoles at the photon point 
((5^=0). At the resonance position f^N{WR) — 1. The phase (t){W) in Eq. (3.66) is introduced 
to adjust the phase of the total multipole to equal to the corresponding t:N phase shift 5^ below 
the two pion threshold. For the s- and p- waves such unitarization procedure is applied up to 
W = 1400 MeV. At higher energies for these waves the phase (f)(W) is taken as constants. Note 
that in the case of the A(1232) excitation, the phase = at W=1232 MeV. In this case, the 
Fermi-Watson theorem is applied up to 1^ < 1600 MeV, where the inelasticity parameter rja of 
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pion scattering amplitude is still close to 1 . For the resonances in the d- and /-waves resonance 
phase (f) assumed to be constant and determined from the best fit. 

While in the original version of MAID [92] only the 7 most important nucleon resonances 
were included with mostly only transverse e.m. couplings, in the new extended version MAID05 
all four star resonances below W — 2 GeV are included (with transverse electric, — Ei±, 
magnetic, = Mi±, and coulomb ,A^ = Si±) couplings. They are P33(1232), Pii(1440), 
Di3(1520), 5ii(1535), 53i(1620), 5n(1650), /^i5(1675), Fi5(1680), AssllTOO), Pi3(1720), 
^35(1905), P3i(1910) and P37(1950). For all these 13 resonances the total numbers of the e.m. 
couplings is 52 (34 for proton and 18 for neutron channels) They are independent of the total 
energy and depend only on Q^. They can be taken as constants in a single-Q^ analysis, e.g. in 
photoproduction, where = but also at any fixed Q^, where enough data with W and 9 
variation is available. Alternatively they can also be parametrized as functions of Q'^ 

3.4 Dispersion relation approach 

The dispersion relations approach, or "S-matrix theory" was developed mostly in the 50's 
and 60's of the last century. It was suggested as an alternative to the field theoretical approach 
to the strong interaction to avoid an expansion in the strong coupling constant. The underlying 
postulate of this theory is that the S-matrix is a Lorentz-invariant analytic function of all mo- 
mentum variables with only those singularities required by unitarity. It is generally believed that 
the analyticity arises from the causality property as in the case of Kramers-Kronig relations for 
the dielectric constant. 

The starting point of the dispersion relation approach is the fixed-t dispersion relations which 
follow from the assumed analytic properties, and from crossing symmetry Eq. (A. 17). They can 
be conveniently written in the following matrix notation as [208]: 



where Sthr — {t^-k + and A^°^^{s^ t, Q"^) is the Bom term evaluated in pseudoscalar cou- 

pling as can be found in [167,208]. It consists of s- and w-channel nucleon poles, and t-channel 
pion pole with residues proportional to the electromagnetic form factors of the nucleon and 
pion. 

The next step is to apply the multipole projection to the dispersion relations Eq. (3.70) 
which leads to the following set of coupled integral equations for the CGLN multipoles M.a — 



P r 

Re A{s, t, Q^) = AP'^^s, t, Q') + - / 



ds' + [elv^l ImA{s',t,Q''), (3.70) 



{Ei±, Mi±, Li^/u), 



Re Ma{W) = M^"^'{W) + M^'^'iW) 




(3.71) 



where a and (3 are the set of quantum numbers and 
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dW 



,ImMa{W')rc,iW') 
{W - W)ra{W) 




(3.72) 



The detailed expressions for the kernels K's and the kinematical factor r^iW) are given in 
Refs. [167,209]. 

One of the methods widely used to calculate the dispersion integrals in Eqs. (3.71)-(3.72) is 
based on the Fermi-Watson theorem [198] such that below the two-pion threshold, we can use 
the following relation between the real and imaginary parts of the amplitude: 



If one further makes an assumption about the high-energy behavior of the multipole phases, 
we obtain a system of coupled integral equations for ReM.a{W). This is the standard method 
to apply fixed-i dispersion relations to pion photoproduction at threshold and in the A(1232) 
resonance region, which was successfully used by many authors [165-168]. The reliability of 
this method at low energies (W < 1400 MeV) is mainly based on the finding that Eq. (3.73) can 
be applied to the important P33 multipole, dominated by the A(1232) resonance contribution, 
with good accuracy up to = 1600 MeV. 

Another method to calculate the dispersion integrals is based on isobaric models [210-213] 
which allow extending the use of fixed-i DR to higher energies. With this approach, the imagi- 
nary parts of the pion photo- and electroproduction multipoles are expressed in terms of back- 
ground (M.^) and resonance (M.^) contributions. 



In a recent work [169], both parts were modeled similar to MAID and good agreement with the 
data is found for the neutral pion photoproduction at threshold. 



ImMa{W,Q^) ^ ReMa{W,Q^) temSaiW) . 



(3.73) 



Im Ma{W, Q^) = Im M^{W, Q^) + Im M^{W, Q^) . 



(3.74) 
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4 Chiral effective-field tlieory in tlie A-resonance region 

In this section we review the recent extension of chiral perturbation theory into the A- 
resonance region. Such an extension yields the opportunity to study the pion-production pro- 
cesses in the A-resonance region within the systematic framework of effective field theory 
(EFT). Some of the applications that are relevant to the 7A^A transition will be discussed here 
as well. 

In Sect. 4. 1 we remind the reader of some basic facts about the chiral Lagrangians for the pion 
and nucleon fields. In Sect. 4.2, we introduce the spin- 3/2 formalism of Rarita and Schwinger 

and discuss its consistency with respect to the description of a spin-3/2 particle. In Sect. 4.3, the 
inclusion of the spin-3/2 A-isobar field in the chiral Lagrangian is discussed with the emphasis 
on the consistency of the spin-3/2 theory. In Sect. 4.4, we discuss the power counting in the pres- 
ence of the A-resonance, and thus introduce the "5-expansion". In Sect. 4.5, we demonstrate 
several applications of the 5-expansion in the processes involving the A-resonance excitation. 

4. 1 Effective chiral Lagrangians 

The relevant effective Lagrangian of low-energy QCD will in our case include the pion, the 
nucleon, and the A-isobar fields. The color degrees of freedom, as well as heavier mesons and 
baryons are assumed to be "integrated out", thus setting an upper limit on the energy range 
where this theory is valid. To begin with we need to write down the most general Lagrangian 
involving these fields and consistent with the underlying symmetries of QCD. In particular, the 
chiral symmetry is known to govern the interaction of hadrons at low energies. 

Let us briefly recall here that chiral symmetry is a symmetry of massless quark Lagrangian, 



where qf is the quark field, D is the covariant QCD derivative, and the summation over the flavor 
index / is understood. This Lagrangian is invariant under the following rotations of quark fields 
in the flavor space 



where index a — l^rij — 1, with n/ being the number of flavors, ^'s are 2 x {rij — 1) independent 
parameters and r's are the SU(n/) Pauli matrices. Since the rotations are done independently 
for the left-handed and right-handed quarks (9l Or), we have a global SUL(n/)xSUij(n/) 
symmetry, known as chiral symmetry. 

In QCD chiral symmetry is broken both spontaneously, by non-perturbative effects of the 
QCD vacuum, and explicitly, by the quark masses. Both mechanisms break chiral symmetry 
down to the SU(n/) symmetry under rotation Eq. (4.2) with 9l — Or, the isospin symmetry. 
The spontaneous chiral-symmetry breaking (xSB) leads to the generation of the chiral quark 

More details and insights can be found in Refs. [214—216] and references therein. 



(4.1) 



i(l - 75) exp {lOlT^f,) + i(l + 75) exp (^^^r;^OJ , 



(4.2) 
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condensate, {qq) ~ —(230 MeV)^, and to the appearance of — 1 massless modes — Gold- 
stone bosons (GBs). The expUcit xSB mechanism is responsible for giving the mass to the GBs. 
Accordingly, the effective low-energy Lagrangian of QCD should contain GB fields, and pre- 
serve explicit chiral symmetry up to the terms vanishing in the limit of massless GBs masses — 
the "chiral Umit". 

In what follows we restrict our consideration to QCD with two flavors (uf = 2), up and down 
quarks only. In this case the three Goldstone bosons are the pions, described by a (pseudo)scalar- 
isovector field tt". The chiral Lagrangian can conveniently be written in terms of the unimodular 
non-linear representation of the pion field: 



U{x) = e2'"(^)/^% with TT = iTrV^ = - 

2 \^/27r- -7r° J 



(4.3) 



where n^'^ represent the charge eigenstates of the pion. Under the chiral transformation, U — > 
RUL\ with L and R elements of SUl(2) and SUi^(2), respectively. Therefore, the Lagrangian 
containing an even number of ?7's is chirally symmetric. It is also obvious that only derivatives 
of U can enter such chirally-symmetric terms. (Terms where U enters without derivatives re- 
duce trivially to terms with derivatives, e.g., Tx\U'^ duUd^W U] = Tr[dij,Ud^W].) This shows 
that the chirally-symmetric interactions are proportional to the momentum of Goldstone bosons, 
and therefore is weak at low energies — the fact that plays a crucial role in the construction of 
the perturbative expansion in powers of momenta, which is utilized in xPT. 

In applications that furthermore will be considered, the isospin-breaking effects are negligi- 
ble. We therefore shall assume exact isospin symmetry, and in particular take the masses of u 
and d quarks to be equal (m^ — ma = iriq). The lowest-order chiral Lagrangian of pion fields 
is then given by 

4^) = lf^TT[d,Ud^U'' + 2Bm, {U + U% (4.4) 

where, at this lowest order, — 92.4 MeV is the pion decay constant and B is related to the 
scalar quark condensate as S = — {qq) / f^. The explicit xSB, linear in U, term gives rise to 
the pion mass. Expanding in the pion field, 

= ^d^Ti'^dV - im^TT^ + 0(71^), (4.5) 

one can identify = 2Bmg, which is the celebrated Gell-Mann-Oaks-Rennner relation. 

Next we include the nucleon, described by an isodoublet spinor field, N = {p, n)^ . Its left- 
and right-handed components, N^^ji = i(l =i=75)A^, transform under SUL(2)xSU/i;(2) as — > 
LNl and Nji RNji. In this representation it is not easy to write down a chirally symmetric 
Lagrangian. In particular, the nucleon mass term, NN — Mn{NlNji -\- NrNi), breaks the 
symmetry, while should preserve it exactly (because does not vanish in the chiral limit). 
The trick is to redefine the nucleon field such that the its left- and right-handed components 
transform in the same way, which is achieved by redefining N'i^ = uNl and N'j^ = Nr, with 



In the notation £^'^\ the superscript indicates the order of the chiral Lagrangian, given here by the 
number of derivatives of Goldstone-boson fields and insertions of their mass. 
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u = VU. Then, N'^ j^ K N'^ j^, where 

K ^^/Ru^/Lu^ = y/Zu^y/Ru (4.6) 

is an SUy(2) matrix which depends on the pion field. The latter fact demands a little more work 
for the derivative of the nucleon field, which now transforms as (omitting from now on the 
prime on the redefined nucleon field): di^,N Kd^N+{d^K)N. Obviously for the Lagrangian 
construction it is desirable to have a derivative which transforms in the same way as the field, 
i.e., D^N — > KD^N. Thus, one is led to consider 

= \fR {d^u) v) y/W + Ku (di^u^) . (4.7) 

Noting that {d^j,u)v} + u{d^j,v}) — 0, it is useful to introduce the SU(2) vector and axial-vector 
currents, 

= i rXi^) ^^.{u d^u'' + u^'d^u) , (4.8a) 
= i T«a;(x) = i (ut d^u - u d^v)) , (4.8b) 
and observe that, under SU^ xSUr, they transform as 



v^^Kv^K^ + i{d^K)K\ (4.9a) 
a^^Ka^Kl (4.9b) 

Therefore, the "chiral covariant derivative" of the nucleon field can be defined as follows: 

D^N ^ d^N + iv^N ^ K D^N. (4. 10) 

The derivative of the axial-vector field can also be defined in a covariant fashion: 

D^a^ = df,a^^i[v^,, aj\ KV^^a^K^. (4.11) 



We are now in position to write down a chiral Lagrangian with nucleon fields. Any hermitian 
Lagrangian built from a combination of the nucleon and axial-vector fields, as well as their 
covariant derivatives, will be chirally symmetric. The lowest-order such Lagrangian is given 
by: 

£W ^N{ip-Mr, + gA4 Ts) N , (4.12) 
where qa — 1.267 is the nucleon axial-coupling constant. 

In practice these Lagrangians need to be expanded in the pion field. To do that conveniently 
we write 

u^ui + i{7[/fn)u2, (4.13) 
where ui^2 are real functions of tt^ = 47r^: 
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-i = -sg = E7^ -175 > (4.14) 



1 / ^2 



oo 2 

In terms of these functions, the vector and axial currents are given by 



«2 = ^sin^ = V-- TTT -TT^ . (4.15) 



Mi9^Mi + (l//2)7r«2a^(7r«2) = -^T«£"^^7r''(a^7r^)u 



1 

47^ 



2 



4/.^ 



TT 



3! 2/2 5! 3/^4 7! 4/6 



+ . . . 



(4.16) 



at,^Uid^{7TU2) - {d^Ui)TtU2 = — T^a^TT*' [(5"'' -^1-^2 + (ttV/tT^) (1 - U1U2) 

2/7r 



2/. " 



TT TT 

+ 



^3!/2 5!/4 7!/6 



+ . . . 



(4.17) 



The currents, obviously, are of the first order in the derivatives of the pion field. 

Finally let us note that in the presence of the electromagnetic field (AfJ, the electric charge of 
the pions is accounted for by making the "minimal substitution": S^tt" d^n"' + e s"**^ A^tt'', 
in the above expressions. Similarly, the proton charge is included by the minimal substitution 
in the chiral derivative as: D^A^ D^N - ie\{l + t^)A^,N. 



4.2 Inclusion of the spin-3/2 fields 



The A(1232) is a spin-3/2 resonance. Therefore its spin content can conveniently be described 
in terms of a Rarita-Schwinger (RS) field [217]: '0^'^'', where ji is the vector and a the spinor 
index; the latter index is omitted in the following. The free Lagrangian of the massive RS field 
is given by 

Cks = {tY"'' da - , (4.18) 

where M is the mass, and the totally-antisymmetric products of 7-matrices are defined as: 
7^^^^ = 5[7'',7l, Y""" = UY^n"} = -^£'"''"^7/375 (using the convention £0123 = +!)• The 
corresponding Euler-Lagrange field equations are: 



iY'^'^do^ip^ - MY^-^u = 0, (4.19a) 

d^{iY'"'da - MYn^u = = Y^'d^^,, (4.19b) 

-^^.{iY"°'^a - MY")^. = = -C^^Y'd^ + 3MY)^., (4.19c) 

which can equivalently be written as: 
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{i^ -M)i,^ = 0, 

a- v = o, 

7 • -0 = 0. 



(4.20) 



Thus, the RS field obeys the Dirac equation, supplemented with the auxiliary conditions, or 
constraints. The constraints ensure that the number of independent components of the vector- 
spinor field is reduced to the physical number of spin degrees of freedom (sDOF). 

Note that the constraints are built into the Lagrangian. This is achieved by making the La- 
grangian to be symmetric under a certain local transformation of the RS field. To exhibit this 
local symmetry, observe that the massless RS Lagrangian is, up to a total derivative, symmetric 
under 

^P^{x)^^^{x) + d^e{x), (4.21) 

where e is a spinor field. This gauge symmetry is a fermionic analog of the gauge symmetry of 
the electromagnetic field, and just as in that case, it leads to a reduction of the number of sDOF 
to 2, as is required for a massless field with a spin. The mass term breaks (partially) this gauge 
symmetry to raise the number of sDOF to 4, as is required for a massive field with spin 3/2. 

Clearly, the coupling of the RS field must be compatible with the free theory construct, in 
order to preserve the physical sDOF content of the theory. However this fact has largely been 
ignored in the literature on the field-theoretic description of the A. There are plenty of examples 
of the so-called "inconsistent" couplings, i.e., couplings that violate the free-theory constraints. 
Besides involving the unphysical sDOF, such couplings lead to fascinating pathologies, such as 
negative-norm states [223,224] and superluminal (acausal) modes [225,226]. 

The quest for "consistent" spin-3/2 couplings was raised from time to time, see, e.g., [227- 
232]. One of the most viable proposals up to date is the one of gauge-invariant couplings [222,233], 
suggesting that the couplings invariant under the gauge transformation (4.21) are consistent. In- 
deed, in the case of gauge-invariant couplings, only the mass term breaks the gauge symmetry, 
hence changing the sDOF content, and it is known to do that in a correct way. 

A notable feature of the gauge-invariant couplings is that the corresponding vertices satisfy a 
trans versality condition: 

p^r'^(p,...) = 0, (4.22) 

where p is the four-momentum and ji is the vector index of a spin-3/2 leg. The RS propagator, 
obtained by inverting the operator in Eq. (4.18), can be written as the following anticommutator : 

= -\ {{i> - M)-\ {g^, - -^p^p, - 17^,)} . (4.23) 



'A canonical method for determination of constraints is due to Dirac [218]. See Refs. [219-222] for 
applications of Dirac's method to the spin-3/2 case. 

In this case the degrees-of-freedom counting goes as follows. The vector-spinor has 16 components, 
the field equations show that there are 8 conditions on them, and the Dirac equation. The latter halves 
the number of independent components, and thus in total we have: (16-8)72=4, the number equal to the 
number of different spin polarizations of massive spin-3/2 particle. 
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It apparently contains a spin- 1/2 sector, that can be made explicit by writing the propagator in 
terms of the co variant spin projection operators: 

^-(f' = ^' + + M)P^:i - {piSl - pS!^1) . (4.24) 
where 

^S^'^ = 9^.u - - ^{i>l^.Pv + P^.l.i>) (4.25) 
projects onto the pure spin-3/2 states, while 



P2l''S^=P^lPv/P^^ 

P{Sl=P'P.l,JiV3p'), (4.26) 
<1^=pX7../(V3/) 

are projection operators onto the spin 1/2 states. It is easy to see that, in combination with gauge- 
invariant couplings satisfying Eq. (4.22), the spin- 1/2 contributions decouple from observables, 
e.g., 

r^ip, . . .) SM r^ip, . . .) = r^ip, . . .) Plli'\p) r^{p, . . .) (4.27) 

Such a decoupling of the lower-spin contribution is of course a desirable effect and often had 
been implemented in the literature "by hand". Namely, one would either drop the spin- 1/2 terms 
in the non-local decomposition of the propagator, Eq. (4.23), see, e.g. [234,235,101]), or use a 
local decomposition, e.g., 

Siiv{p) — — ^ _ ^J^^y ^ ~3M2~ (^'^'^ ~ |7/i7i^) + {liJ,Pv — lixPv) , (4.28) 

and retain only the P^^/^) term therein [236,237,103]. Both methods share a common problem: 
the ad hoc deletion of momentum-dependent terms may destroy the symmetries, such as chiral 
and electromagnetic-gauge invariances. For example, while the full RS propagator is guaran- 
teed to obey a Ward-Takahashi identity with an electromagnetic coupling obtained by 'minimal 
substitution' into Eq. (4.18), the truncated propagator does not even have an inverse. In other 
words, since the above-mentioned procedures are not based on a Lagrangian it is not clear how 
to implement the symmetries. In contrast, the spin-3/2 gauge-invariant couplings ensure the 
spin- 1/2 decoupling automatically. The only question is then how to implement the spin-3/2 
gauge and other symmetries in the same Lagrangian. Prior to attempting to answer this, we 
would like to make one more remark. 

While the spin-3/2 gauge-invariant couplings make sense, perhaps equally consistent seems 
the idea of having couplings which break the gauge symmetry in the same way as the mass 
term. However, such couplings should then be proportional to the mass, in order to provide 
consistency in the massless limit. And the couplings proportional to the mass can be rewritten 
in a gauge-invariant way by using the free-field equation, Eq. (4.19a), or equivalently by a field 
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redefinition. To give an example, consider a coupling of the RS field to a spinor ^ and a scalar 
0: 

Ant = 9M ^Y^'i^i. + H.c, (4.29) 
where is a dimension [mass]~^ coupling constant. This coupling is known to affect the con- 
straints in the same way as the mass term [238,222]. Upon the field redefinition, 

^^,^^^, + 9^d^(f>, (4.30) 

we have C-rs + ^int ^rs + ^Lv ^^^h the new coupling being 

jC/^, = ig-^-f^^'^ido.i;^) + H.c, (4.31) 

which evidently is a gauge-invariant coupling. 

One should emphasize here that the above example is quite exceptional, because in gen- 
eral field redefinitions lead to the appearance of higher-order couplings, the so-called "contact 
terms". Their appearance, however, is not troublesome as long as we deal with effective theories, 
where all possible higher-order term are present anyway, unless forbidden by other symmetries. 
The latter condition brings us again to the question of how the spin- 3/2 gauge symmetry will 
co-exist with other local or global symmetries of the effective theory, and in particular the chiral 
symmetry. 

Certainly it is not easy, in many cases impossible, to incorporate several different symme- 
tries in a given interaction. For example, in the description of electrically charged RS field, the 
only possible way to have both the spin- 3/2 and the electromagnetic gauge symmetries, in a 
closed form, is to allow for general covariance in a de-Sitter geometry, resulting in an extended 
supergravity [239]. In many other cases, e.g., spin-5/2, even such formidable possibilities are 
unavailable. Nevertheless, within the EFT framework one can envision a following method for 
inclusion of the spin-3/2 (and other higher-spin) gauge symmetries. 

To start with, one may construct the effective Lagrangian disregarding the higher-spin gauge 
symmetry. Then, in case of a massive spin-3/2 field, the gauge-invariant couplings can be ob- 
tained by the following substitution [240]: 

^ = {Mjn-'Y''"idai^, (4.32) 

everywhere in the interaction Lagrangian. This step ensures the spin-3/2 gauge symmetry while 
providing an on-shell equivalence of the new and old couplings. Unfortunately in doing so, the 
other symmetries are likely to be violated since, while the field transforms covariantly under 
other symmetries, its derivative does not. To restore the other symmetries one needs to replace 
the derivatives by covariant derivatives ('minimal substitution'): 

d^ipu ^ D^iP,, (4.33) 

hence violating again the spin-3/2 gauge symmetry. 

To break out of this loop, we need to note that generically = + T^, where the "connec- 
tion" is at least one order higher (in the EFT expansion) than d^. Then, if the first substitution 
restores the spin-3/2 symmetry to, say, order n, the other symmetries will be violated only at or- 
der n -I- 1 or higher. So to order n all the symmetries are satisfied. When going to the next order. 
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we would need to make the minimal substitution to restore the other symmetries to that order 
and the first substitution to restore the spin-3/2 symmetry. One can continue this procedure to 
establish all the symmetries to any given order in the effective expansion. 

We emphasize that substitution (4.32) is certainly not unique, other free equations of the RS 
field can be used, e.g.. 



= {i/M) Y {d.ij, - d.iP^), (4.34) 
r = (1/M) £'^^"^^375 (4.35) 

all leading to gauge-invariant couplings, equivalent at the order where the field equations are 
used. Further on we will apply this method to write down chiral Lagrangians involving the 
A-isobar field. 



4.3 Chiral Lagrangians with A 's 



The spin-3/2 isospin-3/2 A(1232) can be represented by a vector-spinor isoquartet field, 
Afj^ — {A^'^, A'Y, where the four components correspond to the charge states of 

the A-isobar. Inclusion of an isoquartet into the chiral Lagrangian can be done by generalizing 
slightly the formalism of Subsect. 4. 1 for the isodoublet field to the case of isospin 3/2. It ba- 
sically amounts to replacing the SU(2) generators in the fundamental representation (t"/2) by 
the generators in the isospin-3/2 representation. The latter have the following form: 



^3 



2i 
3 



( ^ 

VS/2 1 

1 V3/2 

V3/2 ) 

-VS/2 

V3/2 0-10 
1 -V3/2 

y V3/2 



:diag(l, 



and satisfy T'^T" = 5/3. 
The chiral transformation of the isoquartet field is then given as 



(4.36a) 



(4.36b) 



(4.36c) 



(4.37) 
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where K4 is a four-dimensional SUy(2) matrix. Similarly, we introduce the isospin-3/2 vector 
and axial- vector currents: 



^(3/2) ^ ^a^a(^) ^ Tr(r'^^;^), (4.38a) 

a(f/^) = r«a;(x) = r«Tr(T«a^), (4.38b) 
and the chiral covariant derivative of the A-field: 

D^A, = {d^ + zv^^/^^)A,. (4.39) 



Before writing down the effective Lagrangians using these ingredients, let us remark on the 
other frequently used representation. The so-called isospurion representation is based on an 
isovector-isodoublet field [84,241,242]: = T^A^, with T the isospin-l/2-to-3/2 transition 
matrices defined as 



^ \^ -1 V3 

( Vsol o\ 

1 V3 



rp2 



—I 



7 



rp3 



10 
10 



l^a. b 

Q I ' » 



(4.40a) 



(4.40b) 



(4.40c) 



0. Under a chiral rotation the 



and satisfying T^T^'t = 
isospurion transforms as 

A"^ ^ K"^ K A^, (4.41) 

where K, an SUy(2) matrix given by Eq. (4.6), acts on the doublet components, while K"^^ — 
^Tr{T'^KT^K^) transforms the isovector components. Note that the latter object has the follow- 
ing important properties: K'^^t^ = Kt"-K\ K"-'^K^^ = 5'"^, where 5 is the Kronecker symbol. 
The transformation properties of the isospin-3/2 field in the two different representations are 
related via: 

K4 ^ T"^ K T'' K^K (4.42) 

Knowing these relations, one can easily go from one representation to another. They both, of 
course, are equivalent at the level of observables. 

In writing down the chiral Lagrangian we adopt the isoquartet representation. The first-order 
Lagrangian of the A is given by: 



= A, {ir"' D, - Ma rn - \Ha A, 4^'/'^ 75 ^^ 



(4.43) 



where Ma is the mass of the A-isobar, Ha is the axial coupling constant of the A given, in the 
large- A'c limit, by Ha — {9/5)gA- Note that the electric charge of the A can be accounted for 
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by the following minimal substitution: —^D^~ iei(l + 3T^)A^. 

The chiral interactions in Eq. (4.43) obviously do not have the spin-3/2 gauge symmetry, see 
Eq. (4.21). We follow the program outlined in the previous section to incorporate this constraint. 
First, we write out the spin-3/2 gauge-invariant couplings which are on-shell equivalent to the 
ones in Eq. (4.43), e.g.. 

The 'minimal substitution' into the second term, which is required to restore the chiral symme- 
try, generates higher-order contributions, such as: 



4'^ = -^e'^'"''' idoA,) [^^'/'^ 75 + ^Ha vf'^ + H.C. (4.45) 
The spin-3/2 gauge-symmetry is manifest in the following on-shell equivalent Lagrangian: 



= -J^^'"'"'' [^'^'^ 75 + 4^'^'^] vf^^ Yid,A, - d.A,) + H.C. (4.46) 

This of course is not the complete second-order Lagrangian, only the term required by chiral 
symmetry when the first-order Lagrangian (4.44) is used. 

The terms required by chiral symmetry are not only of higher order in the derivatives of 
the pion field, they are are of higher order in pion field itself, and hence in many cases may 
appear only in multi-loop corrections. Here we shall focus on a single-pion production with no 
more than one-loop contributions, thus it suffices to consider the couplings with no more than 
three pion fields. The relevant terms of the chirally and gauge symmetric Lagrangian of the A, 
expanded to third order in the pion field, read 

= W^^"""'^^'^'' i^A'^) daTl" + 0(71% (4.47) 

Similarly, we write a few relevant couplings of the A^A-transition Lagrangian in the form 
which manifests the spin-3/2 gauge symmetry and is expanded to leading order in the pion 
field: 

^NA = T^^NT'^ Y"' id,A,) ^att" + H.C., (4.48a) 
£gk = JJJ^'^T'Y-' (d./^w") (a„A„) + H.C, (4.48b) 



Ma 



F'"' + H4:.48d) 
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where F'^^ and F^'^ are the electromagnetic field strength and its dual. Note that the electric and 
the Coulomb 7A^A couplings are of one order higher than the magnetic one, because of the 75 
which involves the "small components" of the fermion fields and thus introduces an extra power 
of the 3 -momentum. 

Finally, for the applications below it will also be useful to write out the magnetic moment 
coupling for the nucleon and A fields: 

^N=^N i (4^^ + n) a,^ N F^^ , (4.49a) 
= ^ ^[1 + «:Sf) + 3(1 + K^P) Ta] A. F'^'^ , (4.49b) 

where k,^^\ k,^^^ correspond in both the N and A case with the isoscalar and iso vector anoma- 
lous magnetic moments. Also in both cases the anomalous magnetic moments are defined as 
the deviation from the gyromagnetic ratio (g = jji/ s) from 2, the natural value for an elementary 
particle of any spin s [243-245]. In this notation the magnetic moment of the A corresponds 
with 



/^A = ^ (3eA + \nf + T3 ) (4.50) 

where eA = (1 + 37^) /2 is the charge of the A in units of e. 

We should piont out that such a choice of the anomalous magnetic moment of the A is not 
(yet) widely used. One usually defines it as a deviation of the magnetic moment from the mag- 
neton value: e eA/(2MA). Namely, in Eq. (2.20a), (0) corresponds precisely with the conven- 
tional definition of the anomalous magnetic moment. The relation between the two conventions 
is obvious: K^'""^ = f;^^'''\q) - 2 . 



4.4 Power counting, renormalization and naturalness 



The fact that the strength of the chiral interactions goes with derivatives of pion fields allows 
one to organize a perturbative expansion in powers of pion momentum and mass — the chiral 
perturbation theory [7,8]. The small expansion parameter is p/K^sB, where p is the momentum 
and A^5B ~ 47r/^ ^ 1 GeV stands for the scale of spontaneous chiral symmetry breaking. 
Based on this expansion, one should be able to systematically compute the pion-mass depen- 
dence of static quantities, such as nucleon mass, magnetic moments, as well as the momentum 
dependence of scattering processes, such as pion-pion and pion-nucleon scattering. It generi- 
cally is an effective-field theory (EFT) expansion, in this case a low-energy expansion of QCD. 
One expects to obtain exactly the same answers as from QCD directly, provided the low-energy 
constants (LECs) — the parameters of the effective Lagrangian — are known, either from ex- 
periment or from QCD itself. 

One of the principal ingredients of an EFT expansion is power counting. The power counting 
scheme assigns an order to Feynman graphs arising in loopwise expansion of the amplitudes. 
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(a) (b) 
Fig. 33. Nucleon self-energy and electromagnetic-moment corrections to one loop. 

and thus defines which graphs need to be computed to a given order in the expansion. In a way, 
it simply is a tool to estimate the size of different contributions without doing explicit calcula- 
tions. Of course, the main requirement on a power-counting scheme is that it should estimate 
the relative size of various contributions correctly. So, if a given graph is power counted to be 
of subleading order, but in explicit calculations gives a dominant effect, the power-counting 
scheme fails. 

In xPT with pions and nucleons alone, the power counting for a graph with L loops, 
(Nn) internal pion (nucleon) lines, and Vk vertices from /cth-order Lagrangian, estimates its 
contribution to go as ]9"xpt^ ^ith the power given by [246]: 

n^PT = AL-2K-Nr, + J2 kVk • (4.51) 

k 

For example, the one-loop correction to the nucleon mass given in Fig. 33(a) is characterized 
by L = 1, A^^ = 1, Nn = 1, Vi =2, while the correction to the nucleon electromagnetic vertex 
Fig. 33(b) has L = 1, A^^^ = 2, A^^ = 1, Vi = 2, and V2 = I ( two vrA^iV vertices from and 
the 77r7r vertex from >C(^^). Therefore, both of these graphs count as order p^. 

Another important ingredient is the renormalization program. First of all, although the in- 
teractions are non-renormalizeable in the usual sense, one always should be able to absorb the 
infinities into the available LECs. This is ensured simply by the fact that the effective Lagrangian 
contains all the possible terms allowed by symmetries and entering with arbitrary coefficients. 
Second, the renormalization program should be compatible with power counting. This issue is 
related to the arbitrariness of renormalization, namely the arbitrariness of the finite part that 
remains after the renormalization. In the EFT framework, it is clear that only the parts which 
satisfy power counting should remain. That is, if power counting estimates the graph to be, say, 
of order p^, but its unrenormalized contribution contains a piece which goes as p^, that piece 
must be completely absorbed into the p^ counter term. Again, symmetries and the generality of 
the effective Lagrangian ensure that this can always be done [247]. 

In the above example of the one-loop correction to the nucleon mass, an explicit calculation 
shows [246] that the loop produces 0(m°) and O(m^) terms, both of which are large. In fact, 
they are ultravioletly divergent, hence are infinitely large. However, the appears of such terms 
is not in violation of the power counting, because there are two LECs: the nucleon mass in the 
chiral limit, M(°\ and c^v, which enter at 0(m°) and O(m^), respectively, and renormalize 
away the large contributions coming from the loop. The renormalized result (up to p"^ terms) is 
given in Eq. (2.47), such that the loop contribution begins at 0{m^) in agreement with power 
counting. 

Similarly, in the case of the chiral correction to the electromagnetic interaction, an explicit 
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calculation will reveal (see, e.g., [248,249]) that the magnetic contribution goes as 



Fig. 33(b) = 



2M, 



e 




( 



9aMjX\ 
47rA ) [ 



Co + Ci 



m. 



'TT 



+ 0{ml) 



(4.52) 



where e and q are the four-vectors of photon polarization and momentum, respectively, and 
where cq j are real coefficients which depend only on isospin. Recalling that both the charge e 
and the photon momentum count as one power of p, we find that the cq term is of order p^, hence 
superficially the power counting is violated. However, there is a LEC coming at this order from 
in Eq. (4.49). The cq term, would it be finite or infinite, can be absorbed by this LEC, thus 
providing correct power-counting for the renormalized loop graph. 

A cornerstone principle of effective field theories in general is naturalness, meaning that 
the (dimensionless) LECs must be of "natural size", i.e., of order of unity. Any significant fine- 
tuning of even a single LEC leads, obviously, to a break-down of the EFT expansion. Therefore, 
if an EFT describes well the experimental data, but at the expense of fine-tuned LECs, the result 
is negative: EFT fails in the description of those data. Such an EFT can still be useful for getting 
insights into the physics beyond the EFT itself. Namely, by looking at the form of the fine-tuned 
operators, one might be able to deduce which contributions are missing. 

For instance, it is well known that the NLO xPT description of the pion-nucleon elastic scat- 
tering, near threshold, requires relatively large values for some of the LECs, see e.g., Ref. [9], It 
is not difficult to see that the operators corresponding with those unnatural LECs can be matched 
to the "integrated out" A-resonance contributions. The problem is that the A is relatively light, 
its excitation energy, A = Ma — M^v ~ 0.3 GeV, is still quite small compared to A^^^ ~ 1 
GeV. Integrating out the A-isobar degrees of freedom corresponds to an expansion in powers 
ofp/A, with p ~ m^, which certainly is not as good of an expansion as the one in the meson 
sector, in powers of p/A^sB- 

The fine-tuning of the "Deltaless" xPT seems to be lifted by the inclusion of an explicit A- 
isobar. Also, the limit of applicability of the EFT expansion is then extended to momenta of 
order of the resonance excitation energy, p ^ A. Such momenta can still be considered as soft, 
as long as A/A^sb can be treated as small. The resulting xPT with pion, nucleon, and A-isobar 
degrees of freedom has two distinct light scales: and A. Perhaps the most straightforward 
way to proceed is to organize a simultaneous expansion in two different small parameters: 
e = m^/A^5B and 6 = A/ A^-sb- However, for power counting purposes, it is certainly more 
convenient to have a single small parameter, and thus a relation between e and 6 is usually im- 
posed. We emphasize that the relation is established only at the level of power counting and 
not in the actual calculations of graphs. In the literature up to date two such relations between 
e and 5 are used: (i) e ~ 5, see [84,86,241,242,104], which we will commonly refer to as the 
"e-expansion", (ii) e ~ 5"^, the "5-expansion" of Ref . [90]. The table below (Table 2) sum- 
marizes the counting of momenta in the three expansions: Deltaless (^xPT), e-expansion, and 
5-expansion. 

An unsatisfactory feature of the e-expansion is that the A-resonance contributions are always 
estimated to be of the same size as the nucleon contributions. In reality (revealed by actually 



The same counting was independently developed by Hanhart and Kaiser [250] in application to the 
pion production in nucleon-nucleon collisions. 
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0{p) 


0{l) 


e-expansion 


0{e) 


0{e) 


(^-expansion 


0(52) 


0{6) 



Table 2 

The counting of momenta in the three different xEFT expansions discussed in the text. 

computing these contributions), they are suppressed at low energies and dominate in the the 
A-resonance region. Thus, apparently the power-counting in the e-expansion overestimates the 
A-contributions at lower energies and underestimates them at the resonance energies. The 5- 
expansion improves on this aspect, as is briefly described in what follows. 

In the ^-expansion, the power counting depends on the energy domain, since in the low-energy 
region (p ~ m^) and the resonance region (p ~ A), the momentum counts differently, see Ta- 
ble 2. This dependence most significantly affects the power counting of the direct resonance 
exchanges, the so-called one-Delta-reducible (ODR) graphs. Figure 34 illustrates examples of 
the ODR graphs for the case of Compton scattering on the nucleon. These graphs are all char- 
acterized by having a number of ODR propagators, each going as 



SoDR 



s-Ml 2Ma p-A' 



(4.53) 



where p is the soft momentum, in this case given by the photon energy. In contrast the nucleon 
propagator in analogous graphs would go simply as ^at ~ 1 /p. 








Fig. 34. Examples of the one-Delta-reducible (1st row) and the one-Delta-irreducible (2nd row) graphs 
in Compton scattering on the nucleon. 

Therefore, in the low-energy region, the A and nucleon propagators would count as 0{l/5) 
and 0{\/5), the A being suppressed by one power of the small parameter as compared to the 
nucleon. In the resonance region, the ODR graphs obviously all become large. Fortunately they 
all can be subsumed, leading to "dressed" ODR graphs with a definite power-counting index. 
Namely, it is not difficult to see that the resummation of the classes of ODR graphs results in 
ODR graphs with only a single ODR propagator of the form 



1 



^ODR 



^ODR ^ 



p-A-S' 



(4.54) 
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where S is the A self-energy. The expansion of the self-energy begins with p^, and hence in 
the low-energy region does not affect the counting of the A contributions. However, in the 
resonance region the self-energy not only ameliorates the divergence of the ODR propagator at 
s — M\ but also determines power-counting index of the propagator. Defining the A-resonance 
region formally as the region of p where 

\p-A\<5^K^SB, (4.55) 

we deduce that an ODR propagator, in this region, counts as 0{l/5^). Note that the nucleon 
propagator in this region counts as 0{\/5), hence is suppressed by two powers as compared to 
ODR propagators. Thus, within the power-counting scheme we have the mechanism for esti- 
mating correctly the relative size of the nucleon and A contributions in the two energy domains. 
In Table 3 we summarize the counting of the nucleon, ODR, and one-Delta-irreducible (ODI) 
propagators in both the e- and (5-expansion. 





e-expansion 


5-expansion 
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1/6 


1/6^ 


SoDI 


1/6 
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Table 3 

The counting for the nucleon, one-Delta-reducible (ODR), and one-Delta-irreducible (ODI) propagators 
in the two different expansion schemes. The counting in the (5-expansion depends on the energy domain. 

We conclude this discussion by giving the general formula for the power-counting index in 
the (5-expansion. The power-counting index, n, of a given graph simply tells us that the graph is 
of the size of For a graph with L loops, Vk vertices of dimension k, pion propaga- 

tors, Npj nucleon propagators, A^a Delta propagators, Nqdr ODR propagators and Nqdi ODI 
propagators (such that A^a = Nqdr + Nqdi) the index is 




2%PT - Na, P ~ m^; 
n^PT - ^NoDR - Nodi , p ^ A, 



where n^pT, given by Eq. (4.51), is the index of the graph in ^PT with no A's. For further 
details on the S counting we refer to Ref. [90]. The rest of this section is devoted to applications 
of xEFT to several processes that are relevant for the — > A transition. 

4.5 In practice: the next-to-leading order calculations 

4.5.1 Pion-nucleon scattering 

The pion-nucleon (ttN) scattering amplitude at leading order in the ^-expansion in the reso- 
nance region, is given by the graph (LO) in Fig. 35. This is an example of an ODR graph and 
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(LO) (NLO) 
Fig. 35. The leading and next-to-leading order graphs of the TrAT-scattering amplitude. 



(a) (b) (c) (d) 

Fig. 36. The leading and next-to-leading order graphs of the A self-energy. 

thus the A-propagator counts as 6~^. The leading-order vertices are from and since p ^ 6, 
the whole graph is 0{6~^). 

At the NLO, the graphs labeled (NLO) Fig. 35 begin to contribute. The ttA^^A vertices de- 

(2) (2) 

noted by dots stand for the Ha coupling from Z^Xta and the circles for the hi coupling from £)vA' 
see Eq. (4.48). The NLO graphs are thus 0{6'^). It is not difficult to see that graphs containing 
the loop correction to the vertex, as well as the nucleon-exchange graphs, begin to contribute at 
N^LO [0{6)]. 

The ODR graphs contribute only to the P33 and D^^ partial waves. The D^s contribution is 
due to the "negative-energy states" contribution and is suppressed by 5^, as compared to the 
"positive-energy states" contributions to P33. The D33 will therefore be omitted from our con- 
siderations. 

The P33 contribution can conveniently be written in terms of the following partial- wave 'K- 
matrix' : 

1 T(W) 

Kps3 = , (4.56) 

where W = ^/s is the total energy and T is an energy-dependent width, which arises from the A 
self-energy. At this stage it is already taken into account that the real part of the self-energy will 
lead to the mass and field renormalization and otherwise are of N^LO. Thus, only the imaginary 
part of the self-energy affects the NLO calculation. 

In the ODR graphs of Fig. 35, the A-propagator is dressed by the self-energy given to NLO 
by the graphs in Fig. 36. The Lorentz-covariant self-energy of the spin-3/2 A is a rank-2 tensor, 
has the 7-matrix structure and is dependent on one four-momentum. As such, this object can 
in general be decomposed onto a basis of 10 independent scalar functions. The procedure for 
computing these functions is cumbersome [251,252]. A consistent procedure of "dressing" the 
Rarita-Schwinger propagator with such a self-energy does not exist, because the propagator is 
derived from the Lagrangian with only two independent structures: the kinetic and the mass 
term. As discussed above, the form of the Rarita-Schwinger Lagrangian is constrained by the 
requirement of physical spin degrees-of -freedom count, and hence the expression for the self- 
energy must be constrained as well. The spin-3/2 gauge-invariant couplings, introduced in the 
previous two subsections, automatically ensure these constraints. In that case, the most general 
form of the self-energy can be written as 
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where p(3/2) is the covariant spin-3/2 projection operator, and has the spin-1/2 Lorentz 
form, thus 2 independent Lorentz structures as needed. 
The "dressed" propagator of the A to NLO is then given as follows: 



{j, - Ma)[1 - Z2 ^'(Ma)] - {Mf^ -M^ + Z2 S{Ma) ) 



-PiT\p), (4.58) 



where is the "bare" mass, Z2 is the field renormalization constant, and 



Our 'on-shell renormalization' conditions read: 



(4.59) 



Z2 = l- Z2ReS'{MA), 
= M^^^ -MA + Z2 RqS{Ma), 



and thus the renormalized NLO propagator is given by 

Sapip) = 7 



Ma) [1 - i Imr'(MA)] - i Imr(MA) 
The energy-dependent width in Eq. (4.56) is then given by 

r(W) = -2Im[i;(MA) + (W- Ma) S'(Ma)] 
and therefore the expression for the K-matrix becomes 

Imr(MA) 



K 



P33 



W-Ma 



+ Imi;'(MA) 



(4.60) 
(4.61) 



(4.62) 



(4.63) 



(4.64) 



An elementary calculation (see, e.g., [109]) of the graphs in Fig. 36 yields, in the region 

W e [Mn + myr, Ma + mj^], the following result: 



Imr(MA) = -7r 
Im U'{Ma) = 



/,2 I 9/, 1, A 

^ + Ma)' - mlf/\A' - mlf' , 



24Mi (SttA) 



^Mi(87r^): 



:[(M^ + MaY - mlf'sjA-^-ml 



X 



(4.65a) 



Mi - {M% - ml)'' - \{A' - ml){Ml + MaM^ - ml) . (4.65b) 



The -kN scattering phase-shift is related to the partial-wave K-matrix simply as 

5i — arctan Ki , 



(4.66) 
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where / stands for the conserved quantum numbers: spin (J), isospin (J) and parity (P). We 
emphasize that the P33 phase (corresponding to J = 3/2 = /, P = +) is the only nonvanishing 
one at NLO in the resonance region, and is computed by substituting the NLO expressions for 
the self-energy [Eq. (4.65)] into Eq. (4.64). We can then fix the LECs Ha and hi by fitting the 
result to the well-established empirical information about this phase-shift. 

Thus, in Fig. 37 the red solid curve shows the NLO description of the empirical P33 phase- 
shift represented by the data points. The curve is obtained by taking Ha = 2.85 and hi = 0, 
and is characterized by x^/point ~ 1, where we assume 0.5 degree uncertainty in the empirical 
values. (The best fit is obtained by slightly decreasing Ha and increasing hi to about 0.6, how- 
ever the improvement is very small and we prefer to neglect hi for simplicity). The blue dashed 
line in Fig. 37 shows the LO result, obtained by neglecting ImZ" and hi, and taking hA = 2.85. 
This corresponds with the so-called "constant width approximation". At both LO and NLO, the 
resonance width takes the value P(Ma) ^115 MeV. 




1.17 1.19 1.21 1.2321.25 1.27 1.29 

W [GeV] 



Fig. 37. (Color online) The energy-dependence of the P33 phase-shift of elastic pion-nucleon scattering 
in the A-resonance region. The red solid (blue dashed) curve represents the NLO (LO) result. The data 
points are from the SP06 SAID analysis [93]. 

Of interest is the position of the resonance in the complex plane, = M^°^^^ — {i /2) pjf . 
In this NLO calculation, the "complex pole" parameters come out to be 



Mr^=M^ - ^ ^ ^^^^^^^ ^ 1.211GeV, (4.67a) 

All these numbers are within the range quoted by the Particle Data Group [15]. 
Note that the calculations presented here satisfy (the two-body -nN) unitarity exactly. Indeed, 
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Fig. 38. Diagrams for the eN —>■ eirN reaction to LO and NLO in the 5-expansion. The dots denote the 
vertices from the Ist-order Lagrangian, while the circles are the vertices from the 2nd order Lagrangian 
{e.g., the 7A^A-vertex in the first two graphs is the gM coupling from C^"^^). 
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Fig. 39. The ^yNA vertex at 0{S^). The sliced vertex (1) stands for the gE and gc couplings from the 
3rd order Lagrangian. The wiggly double line in (2) stands for the vector-meson propagator. 

the partial-wave S-matrix obtained from the graphs in Fig. 35 is given by 

1+tK, 



Si 



(4.68) 



l-iKi 

where Ki and 6i, given respectively by Eqs. (4.56) and (4.66), are real numbers, hence l^^l = 1. 



4.5.2 Pion photo- and electroproduction 

We now turn to the analysis of the pion electroproduction process. Since we are using the 
one-photon-exchange approximation (until Sect. 6), the pion photoproduction can be viewed as 
the particular case of electroproduction at = 0. 

The pion electroproduction amplitude to NLO in the 5-expansion, in the resonance region, is 
given by the graphs in Fig. 38, where the shaded blob in the 3rd graph denotes the NLO 7A'^A 
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vertex, given by the graphs in Fig. 39. The 1st graph in Fig. 38 enters at the LO, which here is 
0{5~^). The A self-energy in this graph is included to NLO, see Fig. 36. All the other graphs 
in Fig. 38 are of NLO= 0(5°). Note that the A-resonance contribution at NLO is obtained by 
going to NLO in either the vrA^A vertex (2nd graph) or the 7A^A vertex (3rd graph). Accord- 
ingly, the A self-energy in these graphs is included, respectively, to NLO (Fig. 36) and to LO 
[Fig. 36(a)] in the vrA^A coupling. 

The vector-meson diagram. Fig. 39(2), contributes to NLO for ~ AA. One includes it 
effectively by giving the (7jv/-term a dipole Q^-dependence (in analogy to how it is usually done 
for the nucleon isovector form factor): 

The analogous effect for the gE and gc couplings begins at N^LO. 

An important observation is that at = only the imaginary part (unitarity cut) of the 
loop graphs in Fig. 39 contributes to the NLO amplitude. Their real-part contributions, after the 
renormalization of the LECs, begin to contribute at N^LO, for <C AA^sb- At present we 
will consider only the NLO calculation where the vrA-loop contributions to the 7 A- vertex 
are omitted since they do not give the imaginary contributions in the A-resonance region. We 
emphasize that such loops might become important at this order for ~ AA^sb ~ 0.3 GeV^ 
and should be included for the complete NLO result. The present calculation is thus restricted 
to values < 0.3 GeV^. 




W (GeV) 



Fig. 40. (Color online) The M^^ multipole for pion photoproduction as function of the invariant en- 
ergy. Green dashed curves: A contribution without the 7A^A-vertex loop corrections, [i.e., only the first 
three graphs in Fig. 38 with Fig. 39(1) contribution are taken into account]. Blue dotted curves: adding 
the Born contributions, 2nd line in Fig. 38, to the dashed curves. Black solid curves: the NLO calcula- 
tion, includes all graphs in Fig. 38 as well as the loop corrections. The data point are from the SAID 
analysis (FA04K) [93] (red circles), and from the MAID 2003 analysis [92] (blue squares). 



97 




1.18 1.2 1.22 1.24 1.26 1.28 

W (GeV) 



Fig. 41 . (Color online) The El_^ multipole for pion photoproduction. Curve conventions and data points 
are the same as in Fig. 40. 
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Fig. 42. (Color online) The ^j,,: multipole at = as function of the invariant mass W of the vrA^ 
system. Curve conventions are the same as in Fig. 40. 

In Figs. 40-42 we show the results for the pion electroproduction resonant multipoles m[^'^\ 
e[^/'^\ and Si^/"^^ as function of the invariant energy W = ^/s around the resonance posi- 
tion, at g2 = 0. The M^^^^^ and Ef^"^^ multipoles are well established by the MAID [92] and 
SAID [93] partial-wave solutions, thus allowing one to fit two of the three 7A^A LECs at this 
order as: qm = 2.97, gE = —1.0. The third LEC is adjusted to for a best description of the pion 
electroproduction data at low Q'^ (see Sect. 5), yielding gc = —2.6. The latter values translate 
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into G*M = 3.04, 0% = 0.07, and G*c = 1.00 for the Jones-Scadron forai-factors at = 0. 
As is seen from the figure, the NLO results (solid curves) give a good description of the energy 
dependence of the resonant multipoles in a window of 100 MeV around the A-resonance posi- 
tion. Also, these values yield Rem = —2.2 % and Rsm = —3.4 %. 

The dashed curves in these figures show the contribution of the A-resonant diagram of Fig. 38 
without the NLO loop corrections in Fig. 39. For the Mi+ multipole this is the LO and part of 
the NLO contributions. For the Ei+ and 5*1+ multipole the LO contribution is absent (recall that 
gE and gc coupling are of one order higher than the gM coupling). Hence, the dashed curve 
represents a partial NLO contribution to Ei^ and Si+. 

Note that such a purely resonant contribution without the loop corrections satisfies unitarity 
in the sense of the Fermi-Watson theorem [198], which states that the phase of a pion elec- 
troproduction amplitude A4i is given by the corresponding pion-nucleon phase-shift: Aii = 
\Aii\ exp(i6i). As a direct consequence of this theorem, the real-part of the resonant multipoles 
must vanish at the resonance position, where the phase-shift crosses 90 degrees. 

Upon adding the non-resonant Born graphs (2nd line in Fig. 38) to the dashed curves, one 
obtains the dotted curves. The non-resonant contributions are purely real at this order and hence 
the imaginary part of the multipoles do not change. While this is consistent with unitarity for the 
non-resonant multipoles (recall that the non-resonant phase-shifts are zero at NLO), the Fermi- 
Watson theorem in the resonant channels is violated. In particular, one sees that the real parts 
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Fig. 43. (Color online) Non-resonant multipole amplitudes (in units /m.,^) of pion photoproduction 
as function of the invariant mass W of the vrA^ system. The solid curves result from our NLO calcula- 
tion. The data points are from the SAID analysis (FA04K) [93] (red circles), and from the MAID 2003 
analysis [92] (blue squares). 
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Fig. 44. (Color online) The resonant multipoles of pion electroproduction as function at the 
A-resonance position. The green dotted curve is the A-contribution alone. The solid curves are results 
of the NLO calculation of Ref. [89]. Also shown are results of the SAID analysis (FA04K) [93] (black 
dashed-dotted curves), and the MAID 2003 analysis [92] (blue dashed curves). The data points are from 
BATES at = o.l27 [80] andMAMI: = [23] , = 0.06 [79], = o.l21 [78], = o.2 [77]. 

of the resonant multipoles now fail to cross zero at the resonance position. The complete NLO 
calculation, shown by the solid curves in the figure includes in addition the TrA^-loop corrections 
in Fig. 39, which obviously restore unitarity. The Fermi-Watson theorem is satisfied exactly in 
this calculation. 

Next we examine the results for the non-resonant multipoles, which all receive contributions 
of the Bom graphs only. In Fig. 43, we show the NLO calculations for the non-resonant s-, 
p- and (i-wave pion photoproduction multipoles in the A (1232) region in comparison with the 
two state-of-the-art phenomenological multipole solutions IMAID and SAID. At this order the 
non-resonant multipoles are purely real. The multipole solutions show indeed that the imagi- 
nary parts of non-resonant multipoles, around the A resonance, are negligible in comparison to 
their real parts. 

From the figure we conclude that for most of the non-resonant multipoles, the parameter-free 
NLO results, agree fairly well with the phenomenological solutions. The largest differences are 
observed for M^^"^^^ multipole. This multipole corresponds with nucleon quantum numbers. 
The cause of the appreciable difference in this channel is largely due to the nucleon anomalous 
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magnetic moment contributions, which are not included in this calculation (since they appear at 
N^LO), but which are included in the phenomenological solutions. 

Finally, in Fig. 44 we show the NLO xEFT results for the dependence of the resonant 
multipoles, at the resonance position. The red solid curve is with and the green dotted curve 
without the chiral-loop corrections of Fig. 39. The blue dashed curves and the black dashed- 
dotted curves represent the results of MAID and SAID, respectively. The data point are from 
the recent MIT-Bates and MAMI experiments. We observe from the figure that the chiral loops 
play a crucial role in the low momentum-transfer dependence of the Rem ratio. The effect of 
the "pion cloud" is most pronounced in the E2 ^NA transition. 



4.5.3 Compton scattering 

The 7iVA transition plays an important role in the 'yN process in the A-resonance 

region. Already well below the resonance the effects of the A excitation are appreciable, and 
mostly due to its significant contributions to the magnetic and the backward spin polarizabilities 
of the nucleon (cf. [72] and references therein): 



^ ^"'^"^^^ 6 [x 10-^ fm3] , (4.70) 



7f ^ - "r-^/As.. fl + #^ ) - 9 [X fm^] . (4.71) 



{Mn + Ma) a 

2 

These numbers can for instance be compared with the known empirical values for the proton 
case [253,254]: = 3.2 ± 1.2 [x 10"^ fm^], 7(p-^^p) = -38.7 ± 1.8 [xlQ-^ fm^], to 

convince oneself that the effects are significant. 

Note also that the leading chiral-loop contribution to magnetic polarizability is /J^P"^*^) ~ 1.2 
in the heavy-baryon xPT [255] and I3^^~^^'^^ ~ —2 in the manifestly covariant xPT [256,257]. 
The latter number is certainly preferred, if we are to reconcile the large A contribution to 
the magnetic polarizability with experiment. More discussion of this issue can be found in 
Ref. [258]. 

The calculations of Compton scattering in xPT without A's are shown to be applicable only 
to energies not far above the pion-production threshold, W = M^ + m^r, see Refs. [259,260]. 
Traditionally calculations of this reaction in the resonance region are done using isobar-type 
of models [175,176,261] or dispersion-relation approaches [262,263]. Only relatively recently, 
first attempts to compute the A-resonance region in a xEFT framework have been completed to 
some degree of success [90,264,265]. 

In Fig. 45 we show the results of Ref. [90] for the Compton-scattering observables, namely, 
differential cross-section and the linear beam asymmetry. The results are shown as a function 
of the photon lab energy, u — {s — M^)/(2M;v), and the scattering angle in the center-of- 
mass system 6 cm. ■ The complete calculation (red solid curves) represents the NLO result in the 
5-expansion over both the low-energy and the resonance energy region. It is instructive to com- 
pare it to the NLO calculation in the Deltaless HBxPT (blue dashed lines). The latter clearly 
breaks down at energies above the pion-production threshold, uj^hi — ttIt^- These results also 
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Fig. 45. The xEFT calculations of the 7p— > 7p differential cross section (da / dQ^.m.) and the pho- 
ton beam asymmetry (S) compared with recent experimental data: SAL93 [267], LEGS97 [268], 
LEGSOl [25]. Left panel: energy dependence at a fixed scattering angle. Right panel: angular depen- 
dence at a fixed energy. The long-dashed orange line represents the sum of nucleon and pion Born 
graphs, the blue dashed line gives the NLO .^HBxPT prediction, and the red solid line is the full result 
at NLO in the J-expansion [90]. 

demonstrate that the such computed A contributions (inferred by the difference between the 
red solid and blue dashed curves), while being small at low energies, are dominating the res- 
onance region (lu around 340 JVIeV), in qualitative agreement with the power-counting of the 
(5-expansion scheme. 

In all fairness one should note that the NLO results shown in Fig. 45 are obtained using 
gM = 2.6 and = —6.0, the values of LECs which are inconsistent with the photoproduction 
analysis in the same expansion scheme [89] (see also the previous subsection). The source of 
the discrepancy is most likely to be the use of the heavy-baryon expansion in the Compton- 
scattering work [90]. Relativistic effects are expected to be important in the resonance region. 
As mentioned above, they already make a significant impact on such low-energy quantities as 
polarizabilities. A future study of Compton scattering in a manifestly covariant xEFT frame- 
work is called for to clarify this issue. 



4.5.4 Radiative pion photoproduction 

The radiative pion photoproduction (7A^ — > nNy) in the A-resonance region is used to ac- 
cess experimentally the magnetic dipole moment (]V1D]VI) of the A [269]. On the theory side this 
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Fig. 46. Diagrams for the jp ir^pj' reaction at NLO in the 5-expansion. 

reaction had extensively been studied within the isobar type of models (e.g., [270-273]). Here, 
however, we only discuss the more recent study [109] performed within the xEFT framework. 

Consider the amplitude for the 7p — > 7r°p7' reaction. The optimal sensitivity to the MDM 
term is achieved when the incident photon energy is in the vicinity of A, while the outgo- 
ing photon energy is of order of m^. In this case the 7p — > n^pj' amplitude to NLO in the 
(5-expansion is given by the diagrams Fig. 46(a),( b), and (c), where the shaded blobs, in ad- 
dition to the couplings from the chiral Lagrangian, contain the one-loop corrections shown in 
Fig. 46(e), (f). 

The general form of the 7A A vertex is given in Eq. (2.20). To NLO it suffice to keep only the 
F* and F2 terms, which both receive corrections from chiral loops Fig. 46(e), (f). It is important 
to emphasize that the Ward-Takahashi identity, 

g,.r^"^(p',p) = e [{S-'r^{p') - (S-'rf'ip)] , (4.72) 

to NLO leads to the relation F*{0) = 1 - S'{Ma), where E' is given in Eq. (4.65). This 
condition is verified exactly in the NLO calculation [109]. 

Figure 47 shows the pion mass dependence of real and imaginary parts of the A+ and A++ 
MDMs, according to the calculation of Ref. [109]. Each of the two solid curves has a free 
parameter, a counterterm ka from O^' , adjusted to agree with the lattice data at larger values 
of ttItt. As can be seen from Fig. 47, the A MDM develops an imaginary part when m-,^ < A, 
whereas the real part has a pronounced cusp at = A. The dashed-dotted curve in Fig. 47 
shows the result [276] for the magnetic moment of the proton. One can see that /^a+ and fip, 
while having very distinct behavior for small m^, are approximately equal for larger values of 

We next consider the NLO results for the 7p 7r°p7' observables. The NLO calculation 
outlined above, completely fixes the imaginary part of the 7 A A vertex ^° . The expansion for 
the real part of the 7A A begins with LECs from which represent the isoscalar and isovector 
MDM couplings: /t^^ and /t^^ in Eq. (4.49). A linear combination of these parameters, = 
[3+{K^A^ + K^P)/2]{e/2MA), is to be extracted from the 7P TT^p'-f' observables. Several such 
observables are shown in Fig. 48 for an incoming photon energy E!^'^ = 400 MeV as function 
of the emitted photon energy E'^'^-"^-. 

In the soft-photon limit {E'^'^-'^- — > 0), the 7p — > 7r°p7' reaction is completely determined 
from the bremsstrahlung off the initial and final protons. The deviations of the 7p — > 'k'^p^' 
observables, away from the soft-photon limit, will then allow to study the sensitivity to /xa+- It 
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For an alternative recent calculation of the imaginary part of the A MDM see Ref. [266] 
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Fig. 47. Pion mass dependence of the real (solid curves) and imaginary (dashed curves) parts of A++ and 
A"*" MDMs [in nuclear magnetons]. Dashed-dotted curve is the result for the proton magnetic moment 
from Ref. [276]. The experimental data points for A++ and A+ (circles) are the values quoted by the 
PDG [15]. Quenched lattice data are from Refs. [274] (squares) and from Ref. [275] (triangles). 

is therefore very useful to introduce the ratio [273]: 



R 



]_ da 



(4.73) 



7 



where da /dE'^ is the 7p ix^p^' cross section integrated over the pion and photon angles, 
and o"vr is the angular integrated cross section for the 7^ ix^p process weighted with the 
bremsstrahlung factor, as detailed in [273]. 

This ratio R has the property that in the soft-photon limit, the low energy theorem predicts 
that i? — > 1. One firstly sees from Fig. 48 that the EFT calculation exactly satisfies this low- 
energy theorem. Furthermore, the EFT result for R shows clear deviations from unity at higher 
outgoing photon energies, which are in good agreement with the first data for this process [269]. 
The sensitivity of the EFT calculation to the /^a is a very promising setting for the dedicated 
second-generation experiment by the Crystal Ball Coll. at IVIAIVII, which is currently under 
analysis (for first results, see Ref. [277]). It improves upon the statistics of the first experiment 
(Fig. 48) by at least two orders of magnitude and will allow for a reliable extraction of 
using the EFT calculation presented here. 

Besides the cross section, the asymmetries for linearly and circularly polarized incident pho- 
tons are also displayed in Fig. 48. The photon asymmetry for linearly polarized photons, S, at 
E"^ = exactly reduces to the 7]? — > vr^p asymmetry. From Fig. 48 one sees that in the soft- 
photon limit, the NLO calculation is in good agreement with the experimental data point, and 
predicts a nearly constant energy-dependence. It also predicts a weak sensitivity on the IVIDIVI 
value. The linear beam asymmetry is therefore an excellent observable for a consistency check 
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Fig. 48. The outgoing photon energy dependence of the —>■ n^pj' observables for different values 
of fi^+ (in units e/2M/^). Top panel: the ratio of 7p -K^p^' to 7p — > -K^p cross-sections Eq. (4.73). 
Data points are from [269]. Middle panel: the linear-polarization photon asymmetry of the ^p —>■ -K^p^' 
cross-sections differential w.r.t. the outgoing photon energy and pion cm. angle. The data point at E'^ = 
corresponds with the jp — > ir^p photon asymmetry from [23]. Lower panel: the circular-polarization 
photon asymmetry (as defined in [273]), where the outgoing photon angles have been integrated over the 
indicated range. 

of the EFT calculation. 

The asymmetry for circularly polarized photons, Scire, (which is exactly zero for a two body 
process due to reflection symmetry w.r.t. the reaction plane) has been proposed [273] as a 
unique observable to enhance the sensitivity to fx a- Indeed, in the soft-photon limit, where the 
7P rr^py process reduces to a two-body process. Scire is exactly zero. Therefore, its value 
at higher outgoing photon energies is directly proportional to fi^. One sees from Fig. 48 (lower 
panel) that the EFT calculation supports this observation, and shows sizably different asymme- 
tries for different values of /ia- A combined fit of all three observables shown in Fig. 48 will 
therefore allow for a very stringent test of the EFT calculation, which can then be used to extract 
the A+ MDM. 



4.5.5 Errors due to neglect of higher-order effects 

An advantange of the xEFT framework over, say, dynamical models, is that one can make a 
resonable estimate of the theoretical uncertainty of the calculations due to the neglect of higher- 
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Fig. 49. Examples of N^LO contributions to the eN eirN reaction. Open circles denote the e.m. 
coupling to the anomalous magnetic moment of the nucleon. 

order effects. Consider the pion electroproduction case which we have computed above to NLO. 
Some of the neglected next-next-to-leading order (N^LO) contributions are shown in Fig. 49. 
Of course, there is no substitute for an actual calculation of those effects, but at present we only 
know that they must be suppressed by at least one power of 6 (— A/A^sb) as compared to the 
NLO and two powers of 5 as compared to the LO contributions. Therefore, we can estimate 
the size of the N^LO contribution to an amplitude A as: A^i^qS, or Ai^qS"^. The theoretical 
uncertainty of a calculation up to and including NLO can thus be estimated as: 

A^rr = I^Lo + ^NLol (4.74a) 
In cases where the amplitude does not receive any LO contributions, we have simply: 

A^r^\A^i.o\5. (4.74b) 



There are a few caveats in the implementation of such a naive estimate. First of all, as we 
increase the light scales, the expansion parameter must increase too. And this does not concern 
A alone but, in our case, also and Q'^. Therefore, it is more realistic to estimate the error 
using an expansion parameter such as: 



A 



+ 



AT 



M: 



N 



1/2 



+ 



1/2 



(4.75) 



where all the light scales are treated on equal footing and hence are averaged over. 

Secondly, what if the amplitude happens to vanish at some kinematical point. According to 
Eq. (4.74) the theoretical calculation at that point would be perfect, which is of course unlikely 
to be true in reality. So, when considering dependencies on kinematical variable(s), we shall 
take an average of the error over some appropriate region of that variable. 

Given these two points, we are led to the following formula for the theoretical uncertainty of 
the NLO calculation for an amplitude A, 



A. 



\AlJ\ LO^O 
\A\aJ, LO = 0, 



(4.76) 
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and the subscript "av" indicates that the appropriate averaging is performed. 
The theoretical uncertainty of the NLO calculation of an observable O is: 

2\0\aJ, LO = 0, 

where the factor of 2 takes into account that an observable is a product of two amplitudes. 
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5 Results for observables and discussion 



In this Section, we review the current status of pion photo- and electroproduction data in 
the A-resonance region, and discuss the results of the theoretical approaches described above, 
namely dynamical models and xEFT. An extensive review of meson productions from the nucle- 
ons, including single pion, two pions and r] up to second resonance region, has been published 
recently [12]. Another recent review [1 1] on meson productions also goes up to the second res- 
onance region, but is restricted only to photoproduction. 

After some introductory remarks in Sect. 5.1, we shall consider pion photoproduction in 
Sect. 5.2. In Sect. 5.3, the focus will be on the electroproduction case for both the low and 
larger Q"^ regions. 



5.1 Introductory remarks 

The first experimental work on meson photoproduction dates as early as 1949. However, most 
of the present precision data became available thanks to the progress made in accelerator and 
detector technology during the last two decades. At the present time, the new generation of elec- 
tron accelerators at JLab, ELS A, and MAMI are dedicated to the study of nuclear and hadronic 
structure and are equipped with state-of-the-art detector systems. The high intensity electron 
beams, besides being used for electron induced reactions, also provide good photon sources for 
photonuclear reactions via bremsstrahlung. Another way to produce high energy photon beams 
is through Compton backscattering. Laser backscattering systems have been installed at BNL 
(LEGS), ESRF (GRAAL), and at SPringS (LEPS) to produce photon beams to investigate nu- 
cleon resonances and meson photoproduction. 

Turning to the theory, let us recall the key ingredients which enter into the dynamical mod- 
els and the xEFT calculations. Chiral perturbation theory (xPT) is now widely accepted as the 
"basic" theory to describe the low energy interactions of Goldstone bosons among themselves 
and with other hadrons. There is generally good agreement between the xPT predictions and 
experiments [278]. For example, state-of-the-art xPT calculation on n^p at threshold 

[279,280], which has gone up to order p'^ and found large contributions from one-loop charge- 
exchange rescattering, agree beautifully with the new precise measurements [282-287]. 

The xEFT 5-expansion presented in Sect. 4 provides an extension of xPT into the A-resonance 
region. In such scheme, for reactions such as nN scattering and electromagnetic pion produc- 
tion, diagrams with a dressed A in the intermediate states (the ODR graphs) as depicted in 
Fig. 35a and Fig. 38 (upper left diagram) are the leading diagrams . In the NLO, nN scatter- 
ing is still dominated by the A— excitation diagrams as depicted in Fig. 35. For pion electro- 
magnetic production, only the Bom terms with electric ^NN coupling, in addition to the loop 
correction to the 7A^A vertex as given by graphs in Fig. 39 , enter in the NLO. The only ques- 
tion left then is whether the expansion converges. Namely, whether NNLO terms are smaller 



It is interesting to note that the LO xEFT calculation provides the field theoretic justification for the 
previously employed A— saturation model (see e.g. [288,178]). 
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than NLO. At NNLO, for example, the background nN interaction corresponding to the crossed 
Bom diagram would begin to contribute to the A self-energy in the two-loop diagrams which is 
the main driving mechanism in generating the A— resonance in the Chew-Low theory. Conver- 
gence problem aside, EFT in the resonance region does provide a consistent expansion scheme 
which contains the fundamental features of QCD, namely chiral symmetry, gauge invariance, 
covariance, and crossing symmetry, as in the case for xPT for low-energy hadron phenomena. 

As in xPT, dynamical models also start from an effective chiral Lagrangian. However, the cal- 
culation of the ^-matrix in dynamical models is not based on field theory. Rather, the effective 
Lagrangian is used to construct a potential used in a Lippmann-Schwinger type of scattering 
equation. The solution of the scattering equation includes the rescattering effects to all orders 
and thereby unitarity is ensured while the crossing and the other symmetries are not. In dynam- 
ical models applied to pion electromagnetic production, considered in Sect. 3, gauge invariance 
is problematic. The only hope is that the ad hoc schemes which have been proposed provide 
good approximations. The success of DMT for the 7p — > 7r°p at threshold [91] offers a hope in 
this direction. 

Dynamical models [289,191,236,14,290,206] have been able to provide a good description 
of nN scattering lengths and the phase shifts in S, P, and D waves up to 11^ ~ 1.5 GeV, see 
for instance Fig. 31. It should be pointed out here that within the dynamical models, the A res- 
onance is a coherent contribution of an elementary A and Chew-Low graphs as obtained in the 
cloudy bag model [291]. One advantage of dynamical model calculations over xEFT is that they 
provide a unitary framework up to higher energies because they do not rely on a perturbative 
expansion in small energy scales. 



5.2 Pion photoproduction 

The majority of the photopion production data comes from measurements of the differential 
cross section by unpolarized photons on unpolarized nucleons. All other photoproduction ex- 
periments involve the use of polarized photon beams, polarized nucleon targets, a combination 
of both polarization or the measurement of the polarization of the recoil nucleon. 

The earliest of the polarization experiments used unpolarized photon beams and unpolarized 
targets and measured the recoil polarization of the nucleon, P{0), by a scattering of this recoil- 
ing nucleon on a secondary target. High-energy polarized photon beams have been produced 
either by laser backscattering or coherent bremsstrahlung in a crystal. They have been used to 
measure the beam asymmetry, S(^), from an unpolarized target. 

The first data on the target asymmetry, T{6), was obtained in 1972 for the 7 + p — tt"*" + n 
process. Since then an impressive set of measurements of da/dQ, P{9),T{9), and T.{9) has 
been assembled on the final states 7r+n, 7r°p, and 7r~p. The data coverage on the first two chan- 
nels spans the whole resonance region with good angular coverage. The n'p channel data, on 
the other hand, is more scarce since it involves the neutron (deuterium) target. 

Double polarization observables have also been measured. They fall into three categories: 
beam-target, beam-recoil, and target-recoil. Polarized beam-polarized target experiments yield 
G{e) and H{e) and have been measured at Kharkov [301,302] and at MAMI [294,295,303]. 
Beam-recoil and target-recoil measurements are experimentally more difficult as they involve a 
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-2.7 


XEFT 


-133 


-253 


-0.077 


3.48 


-2.2 



Table 4 

Comparison of the values for the helicity amplitudes, Qn^/^, I^'N^a, and Rem extracted from ex- 
periment, with their values in dynamical models (DMT [91], SL [14], DUO [183]) and xEFT calcula- 
tions [88,89]. The numbers within the parenthesis, in the cases of DMT and SL, correspond to the bare 
values. 

second scattering process to analyze the recoil nucleon polarization. No data have been reported 
for them. The most precise determination to date of the Rem ratio in the photo-excitation of 
the A(1232) has come from the simultaneous measurements of p(7,p)7r° and ^(7, n)7r+ with 
polarized photon beams [22,24]. 

A comprehensive account of the status of photoproduction experiments can be found in [93] 
and the full database can be accessed at the SAID website [184]. 

Let us turn to the results of the theoretical approaches considered above. In the xEFT [88,89] 
and in the dynamical model calculations of SL [14], DMT [91], and DUO [183], the jNA cou- 
pling constants are all obtained by fitting to the experiments. We present in Table 4 the resulting 
values for the helicity amplitudes ^41/2, ^3/2. transition electric quadrupole moment Qp^A+, 
transition magnetic dipole moment, iip^A+, and Rem, together with the PDG [15] values and 
experimental results measured at MAMI [22,23] and LEGS [24,25]. The values for Qp^A+ and 
transition magnetic dipole moment Hp^A+ listed in the same row with the PDG values are the 
estimates of Eqs. (2.25) and (2.24) respectively. 

In the DMT and DUO calculations, the best fit to the multipoles obtained in the recent analy- 
ses of the Mainz [167] and VPI groups [205] is shown in Fig. 50. One sees that both models are 
able to describe the energy dependence of the resonant multipoles very accurately throughout 
the A-resonance region. 

Interpreting the results of the dynamical models, one finds that a sizeable part of the Ml 
^NA strength and almost all of the E2 strength is generated by the nN rescattering. One sees 
from Table 4 that the bare values for the helicity amplitudes obtained in both SL and DMT, 
as given within the parentheses, amount to only about 60% of the corresponding dressed val- 
ues, and are close to the predictions of the constituent quark model, as was first pointed out 
in Ref. [14]. The large reduction of the helicity amplitudes from the dressed to the bare ones 
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Fig. 50. Real and imaginary parts of the Mj^_^ and El__^ multipoles. The solid curves are the full 
DMT results, whereas the dashed curves are the results of the DUO model. The blue circles are the 
results from the Mainz dispersion relation analysis [167], whereas the red squares are obtained from the 
VPI analysis [205]. 

results from the fact that the principal value integral part of Eq. (3.23), which represents the 
effects of the off-shell pion rescattering, contributes approximately for half of the Mi_^ ^ , as 
seen in Fig. 30 for the DIVIT model. The large pion rescattering effects obtained in the dynam- 
ical models are generally interpreted as an important manifestation of the presence of the pion 
cloud in the nucleon. In xEFT, the size of the vrA^ contribution to the Ml and E2 transitions 
is quite different from the dynamical models, due to the renormalization of the LECs gM and 
at the resonance position for = 0. The renormalized nN loop contributions then satisfy 
power counting and as a consequence contribute substantially less than in dynamical models. 

(3/2) (3/2) 

The resonant multipoles and i^i^ obtained in the xEFT calculation at NLO are 

shown in Fig. 40 and Fig. 41, respectively. The imaginary part of these multipoles arise from 
purely the A contibutions, see the upper graphs in Fig. 38. The Bom graphs (2nd line of Fig. 38) 
contribute to the real part of the resonant multipoles, as well as to the non-resonant multipoles 
as shown in Fig. 43. Up to NLO in the 5-expansion, the non-resonant multipoles are thus purely 
real, whereas the corresponding nN phase shifts are zero, in accordance with the Fermi- Watson. 

The predictions of xEFT and the dynamical model calculations of SL, DlVlT, and DUO for the 
total and differential cross sections, single polarizations including beam and target asymmetry 
and recoil polarization, and double polarizations G and H for both 'jp — > vr^p and 'jp vr+n 
are shown in Figs. 51-59. 

For the 7p n^p and 7p ix^n processes, one sees from Figs. 51-59 that below the 
resonance the dynamical models are in overall agreement with the available data. Above the 
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Fig. 51. The total cross sections for -jp it"p (left panel) and ■yp — > vr+n (right panel) through the 
A(1232) resonance region as function of the total cm. energy W. Dashed-dotted (black) curves : DMTOl 
model. Dashed (red) curves : SL model. Dotted (green) curves : DUO model. Solid (blue) curves : NLO 
XEFT calculation. The data points are from Refs. [292] (inverted blue triangles), [293] (green triangles), 
and [294] (red circles). 
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Fig. 52. The helicity difference total cross section (J3/2 — cri/2 for 7P ~^ '^^P (left panels) and ^p vr+n 
(right panels) through the A (1232) resonance region. Curve conventions as in Fig. 51. The data points 
are from MAMI [294] (red circles). 

resonance, the SL model starts to fall increasingly below the n^p and in particular the vr+n 
cross section data, see Figs. 51, 53, and 54. For the DUO model at higher energies, one ob- 
serves more structure in the n^p P, T, G, and H observables for forward and backward angles, 
where however the data base is scarce. The DMT model shows good agreement with all 7r°p 
and TT+n observables throughout the A-energy range. 
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Fig. 53. The differential cross section (left panels) and photon asymmetry (right panels) for —>■ n^p 
through the A(1232) resonance region. Curve conventions as in Fig. 51. The data points are from MAMI 
: Refs. [23,296] (red circles) and [295] (red squares). 



The presented NLO xEFT calculations (solid blue curves) are done in the scheme designed 
to work in a small energy domain around the resonance, as discussed in Sect. 4. It is therefore 
of interest to see how far away from the resonance such calculations actually work. One sees 
from Figs. 51-59 that the NLO calculation reproduces well the pion angular dependencies of 
the 7r°p and vr+n observables around resonance. Below the resonance, the xEFT results for the 
7r°p cross sections, and to lesser extent for the vr+n cross sections, start to fall below the data. It 
was checked that this is mainly due to the absence of the anomalous magnetic couplings in the 
Born diagrams and in the irN loop corrections to the 7A^A vertex. Such effects are coming in 
at NNLO in the (5-expansion. For energies about 50 MeV or higher above resonance, the NLO 
begins to overestimate the n^p cross sections, and to a lesser extent the n^n cross sections, 
which also shows up in larger deviations from the data for S (for 7r+n) and T and P (for n^p). 
Nevertheless, given the simplicity of the NLO calculations in the (5-expansion and its potential 
for systematic improvement, one can be pleased with the results. Overall, the xEFT to NLO 
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Fig. 54. The differential cross section (left panels) and photon asymmetry (right panels) for 7p ^ vr+n 
through the A(1232) resonance region. Curve conventions as in Fig. 51. The data points are from MAMI 
: Refs. [23] (red circles) and [295] (red squares). 

in the 5-expansion provides a simple and elegant description of pion photoproduction observ- 
ables in a 100 MeV window around the A-resonance. The precision of the available data can 
be exploited to study the higher order effects in the 7A^A transition using calculations beyond 
NLO. 



5.3 Pion electroproduction 



The pion production by electrons was first observed by Panofsky and his collaborators at 
Stanford in 1955. Even though it was first attempted at Cornell in 1965, it was not until in 
the beginning of 1970's that double-arm or coincidence experiments, in which the scattered 
electron and one of the recoiling hadrons are measured in coincidence, were actively pursued. 
A summary of the coincidence experimental results on pion electroproduction, up until early 
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Fig. 55. The helicity difference differential cross section ~ ^c'"i/2 for IP ~^ P (1^^ panels) and 

7P — > TT^n (right panels) through the A(1232) resonance region. Curve conventions as in Fig. 51. The 
data points are from MAMI [294] (red circles). 



1980's can be found in the extensive review of Ref. [306]. 

Extensive sets of high statistics data of pion electroproduction have been collected at BATES, 
MAMI, and JLab in the last several years. The majority of them are on the ep e'pn^ reaction. 
They cover a large range of invariant mass W, from threshold up to 2.5 GeV and a wide range 
in momentum-transfer squared = 0.05 — 6 GeV^, and the full coverage of pion polar and 
azimuthal angles. Single and double polarization measurements have also been performed. Sin- 
gle polarization experiments are carried out with polarized electron beams to measure the beam 
helicity asymmetry. Double polarization experiments are performed with a polarized electron 
beam together with a polarized target or with recoil nucleon polarization measurements. Data 
on the threshold tt*^ production are measured mostly at low Q'^ region in order to test the chiral 
perturbation theory. They have been obtained prominently at MAMI, see Ref. [307] and refer- 
ences contained therein. 

Data on the 7r+n channel were much less abundant in the past. Most of them are in the mass 
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Fig. 56. The target asymmetry T for —>■ n^p (left panels) and 7p — > vr+n (right panels) through the 
A(1232) resonance region. Curve conventions as in Fig. 51. The data points for 7p — > vrV are from 
Refs. [297] (blue squares) and [298] (red circles). The data points for jp vr+n are from Refs. [299] 
(blue squares) and [300] (red circles). 



range of the first and second resonances and at forward pion cm. angles, with some at backward 
angles. As in the case of photoproduction, data on the n^p channel from a neutron (deuterium) 
target is even more scarce. Currently, CLAS at JLab has an active programs to measure cross 
sections for the vr+n [308] and n^p [309] channels. 

Measurements made before the early 1980's are summarized in Tables 5-7 in Ref. [306]. The 
new electroproduction data in the resonance region, measured recently at BATES, ELS A, JLab, 
and IMAIVII, up until 2003 are given in Table 2 in Ref. [12]. 

Electroproduction of neutral pions on a proton target has been employed most extensively for 
the purpose of determining the Q"^ evolution of the Rem and Rsm ratios. Many experiments 
have been performed at BATES, IMAIVII, and JLab. At BATES, the experiments have been set 
up for central invariant mass W = 1.232 GeV and Q'^ = 0. 127 GeV^. With the use of either un- 
polarized or polarized electron beams and choice of specific kinematics, quantities like induced 
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Fig. 57. The recoil asymmetry P for 'yp —>■ Tr^p (left panels) and 'yp —>■ vr+n (right panels) through the 
A(1232) resonance region. Curve conventions as in Fig. 51. The data points for yp — > ir^p are from 
Ref. [297]. The data points for jp ir^n are from Ref. [299]. 



polarization, and a^, a^x, cttl', and have been measured [310,95,315,80]. Both unpolar- 
ized and polarized experiments are done at IVIAJVII [78,311,77,79] with ranges from a low 
value of 0.06 GeV^ [79] to 0.2 GeV^. In [78], the recoil proton polarization was measured in the 
double-polarized p{e, e'p)n^ reaction and all three proton polarization components were mea- 
sured simultaneously. All three halls at JLab, Hall A [312,142], Hall B (CLAS) [81,313,131], 
and Hall C [130], have programs on 7r° electroproduction with either unpolarized or polarized 
electrons and a large variety of observables have been measured. The experiment of Ref. [131] 
measured values as high as 6.0 GeV^. The aim of such large Q'^ measurements is to map out 
the transition towards the pQCD regime. The large amount of precision data have been used to 
extract the A (1232) transition form factors from = to 6 GeV^. They also provide 

very stringent tests for theoretical models. 

In Figs. 60-62, the different virtual photon absorption cross sections around the resonance 
position are displayed at three different Q"^ values: Q'^ = 0.06, 0.127, 0.20 GeV^, where recent 
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Fig. 58. The double polarization asymmetry G for — > -K^p (left panel) and 7p — > 7r"'"n (right panel) 
atW = 1.232 GeV. Curve conventions as in Fig. 51. The data points are from Kharkov [301] (black 
squares) and from MAMI [303] (red circles) . 
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Fig. 59. The double polarization asymmetry H for 7p 7:^p (left panel) and 'yp —>■ ■K'^n (right panel) 
atW = 1.232 GeV. Curve conventions as in Fig. 51. The data points are from Kharkov [301] (black 
squares). 

precision data are available. We compare these data with the predictions made within the xEFT 
framework [88,89] and using the SL, DJVIT, and DUO dynamical models. 

In the xEFT calculations, the low-energy constants and qe, were fixed from the resonant 
pion photoproduction multipoles. Therefore, the only other low-energy constant from the chiral 
Lagrangian entering the NLO calculation is gc- The main sensitivity on gc enters in axL- A 
best description of the a^L data at low is obtained by choosing gc = —2.6. In Figs. 60-62, 
we also include as an illustration the theoretical uncertainty of the NLO xEFT result, estimated 
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Fig. 60. The pion angular dependence of the 7*p —>■ ir^p cross sections at W = 1.221 GeV and 
= 0.06 GeV2. Dashed-dotted (black) curves : DMTOl model. Dashed (red) curves : SL model. Dotted 
(green) curves : DUO model. Solid (blue) curves : xEFT results of Refs. [88,89]. The bands provide an 
estimate of the theoretical error for the xEFT calculations. Data points are from the MAMI experiment 
of Ref. [79]. 

using Eq. (4.77), where the average is taken over the range of Q^^ . From Figs. 60-62, one sees 
that the NLO xEFT calculation, within its accuracy, is consistent with the experimental data for 
these observables at low Q^. Only for cttl', the xEFT calculation overpredicts the data with 
increasing Q^. This is because this observable is sensitive to the background multipoles i?o+ 
and 5*0+ which are not very accurately reproduced in the NLO xEFT calculation, see Fig. 43. 

The dynamical models are in basic agreement with each other and the data for the transverse 
cross sections. Differences between the models do show up in the axL and axi' cross sections 
which involve the longitudinal amplitude. In particular for utl' the differences reflect to large 
extent how the non-resonant So+ multipole is described in the models. 

In Fig. 63, we compare the data for the resonant multipoles up to the largest measured values 
of with the results obtained for the three dynamical models (SL, DIVIT, DUO) reviewed in 
this work. One sees that all three models give a very accurate account of the dominant m['^'^^ 
multipole. They are also in good agreement with the Rsm ratio, which grows to a value of 
around —25 % around — 6 GeV^. Only for the Rem ratio, the three dynamical models 
show clear differences. While in the DIVIT and DUO models the Rem ratio shows a zero cross- 
ing around ~ 3 GeV^, the SL result is in good agreement with the data, which give small 



Note that cjtl and axu do not receive any LO contributions and therefore the L0= case in Eq. (4.77) 
must be apphed in the estimate. 
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Fig. 61. The pion angular dependence of the 7*p ir^p cross sections at = 1.232 GeV and = 
0.127 GeV^. Curve conventions as in Fig. 60. The bands provide an estimate of the theoretical error for 
the xEFT calculations. Data points are from BATES experiments [95,315,80]. 

negative values for Rem (around —4 %) at the largest measured values of — 6 GeV^. As 
remarked in Sect. 2.8, the data are very far from the pQCD prediction for Rem which should 
approach +100 % at asymptotically large values of Q^. The DMT calculation [91] finds nearly 
as large values for the subdominant (according to pQCD) helicity amplitude A3/2 as for the 
pQCD dominant amplitude A1/2, at accessible values of Q^. In both the calculations and the 
data, one therefore sees no onset of hadron helicity conservation in the 7A^A transition up to 
Q2 ~ 6 GeV2. 



120 



e p -> e p / @ W=1.221 GeV, 2^=0.2 GeV^ 




10 



I I I I I I I I I I I I I I I I I 




I I I I I I I I I I I I I I I I I 



I I I I I I I I I I I I I I I I I 




30 60 90 120 150 180 
0,(deg) 



^ 
C -2 





1 1 1 1 








.V/ • ■ X. 




.■'//' ^ 


■ ■ 


^ ✓ / : /: : 


— V-:v 


^ y.</. -./. y. 






~, ,111 


1 1 1 1 1 1 1 1 1 r 




30 60 90 120 150 180 
0, (deg) 



Fig. 62. The pion angular dependence of the 7*p ir^p cross sections at = 1.221 GeV and = 
0.20 GeV^. Curve conventions as in Fig. 60. The bands provide an estimate of the theoretical error for 
the xEFT calculations. Data points are from the MAMI experiment of Ref. [141]. 
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Fig. 63. Q dependence of the resonant multipoles atW = 1.232 GeV. Upper panel : imaginary part of 

(3/2) 

M{_^ ; middle panel : Rem', lower panel : Rsm- The curves are the results obtained in different dynam- 
ical models : DMTOl (dashed-dotted black curves); SL (dashed red curves); DUO (solid green curves). 
The data are from BATES [80] (blue squares), MAMI [23,79,141] (red circles), JLab/Hall C [130] (open 
circles), JLab/CLAS [81] (black squares) and [131] (inverted triangles), and JLab/Hall A [142] (blue 
stars). 
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6 Beyond the one-photon exchange approximation 



In all the up-to-date studies of the A excitation by electrons {eN eA), as discussed in 
the previous sections, the electromagnetic interaction between the electron and the nucleon 
is assumed to be mediated by a single photon exchange (I7). In this Section, we discuss ef- 
fects beyond this approximation, in particular the effects due to two-photon exchange (27). 
The spectacular discrepancy between the polarization-transfer [152,154,155] and Rosenbluth- 
separation measurements [316] of the nucleon elastic form factors is largely due to 27-exchange 
effects [317-319]. Subsequently, the 27 corrections to the eN — > eA process were studied in 
Ref. [320], where these 27 corrections were evaluated in a parton model. Recently, the 27 cor- 
rections to eA^ eA were studied in a hadronic model [321]. In what follows, we will review 
the general formalism for the 27 corrections to eA^ — > eA and discuss their effect within the 
partonic calculation of Ref. [320]. 



6.1 General formalism 



When describing the eN — > e A process, the total number of helicity amplitudes is 32, and is 
reduced to 16 by using parity invariance. Furthermore, in a gauge theory lepton helicity is con- 
served to all orders in perturbation theory when the lepton mass is zero. Neglecting the lepton 
mass, reduces the number of helicity amplitudes to 8. Recall that in the I7 approximation, there 
are only three independent helicity amplitudes which can be expressed in terms of three form 
factors of the 7*A'^A transition as discussed in Sect. 2.1. 

The 7*A^A transition form factors are usually studied in the eN eirN process in the A- 
resonance region. In the I7 approximation, the 7*A^ txN virtual photon absorption cross 
section da^/dVtT^ was presented in Sect. 3. Beyond the I7 approximation, one can still formally 
define a cross section da^/ dVL^, from the experimentally measured eN — > ei^N five-fold differ- 
ential cross sections using Eq. (3.7), with momentum transfer q defined asq = k — k'. However, 
beyond the I7 approximation F does not have the physical interpretation any more of a virtual 
photon flux factor. The expression of the thus defined cross section da^/d^lj, will then also be 
modified from the I7 expression of Eq. (3.8). For unpolarized nucleons, the cross section at 
W — Ma can, in general, be parametrized as (using the notations introduced in Sect. 3) : 



= (To + e cos(20^) cttt + he sin(20^) aipj^ 

as 

+ V2ee+ cos (pT^aTL + hV2ee- sin (pT^a'j^j;^, (6.1) 



with h = ±1 the lepton helicity and e± = \/l ± e. 

It is convenient to multipole expand Eq. (6. 1) for the ep — > eA+ — > evriV process (i.e. keeping 
only the A multipoles) as : 
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(70 = ^0 + ^(3 COS^ - 1) A2, 

gtt = sin^ 6^^ Co, cttl = 5 sm(26'^) D^, 

g'tt = sin^ Co, (7^L = 5 sin(2e^) D^, (6.2) 

where ^40 can be written as: 

where X denotes the isospin factor which depends on the final state in the A"*" ttN decay 
as : X(7r°p) = 2/3 and X(7r+n) = 1/3, and where an denotes the reduced cross section aR, 
including 27 corrections. 

We are now in position to discuss the corrections to Rem and Rsm as extracted from arr 
and axL- Experimentally, these ratios have been extracted at = Ma using : 



J _ 3A2 - 2Co 17 4MiQ^ 2 , 



^Tm — TTIT^FT" ~ ^SM — Rsm Rem + ■ ■ ■ (6.4) 
C^Ma oAq 

where the omitted terms involve cubic products of Rem and Rsm- These formulas are usually 
applied [81] by neglecting the smaller quantities Rgj^^ and Rem ■ Rsm- Ref- [320], however, 
kept these quadratic terms, and proposed a second method of extracting Rem- 

Tfexp,II _ -(^0 - A2) -2 Co I7 

- 3(Ao-A2)-2Co~^^^- ^^-^^ 

This method apparently avoids corrections at the one-photon level, which may prove to be 
significant at larger momentum transfer due to an appreciable contribution of the Rsm term. 
We denote the corrections to Rem and Rsm by : 

R ~ + SR'^ + 5R^\ (6.6) 
where 5R^'^ denotes the corrections due to the quadratic terms in Eqs. (6.4),(6.5), which are : 

O^EM — ^ q2 q2 ^SMi 

SR'g^^'^O, (6.7) 
<^-^5M — Rem ■ RsM- 



In Fig. 64, we show the effect of the I7 corrections on Rem and Rsm- We see that its effect on 
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Fig. 64. dependence of the I7 corrections, 5R^^ , to Rem (upper panel) and Rsm (lower panel) at 
W = 1.232 GeV. The uncorrected data are from JLab/Hall C [130] (open circles), JLab/CLAS [81] 
(solid squares) and [131] (inverted triangles). The error bars reflect the statistical error only. The red 
circles include the corrections 5R^^j^.^. The triangle at = is the real photon point from [23], whereas 
the star at = 1 GeV^ is the extraction from the JLab/Hall A recoil polarization experiment [142], 
which provides an independent way to extract Rem and Rsm- 



Rsm is negligible, whereas it yields a systematic downward shift of the Rem result. This shift 
becomes more pronounced at larger e, and it is found that at the highest Q"^ values it decreases 
Rem by around 2 %. To avoid such a correction, it calls for extracting Rem according to the 
procedure // as we outlined above. Alternatively, double polarization experiments provide a 
very useful cross check on the extraction of the Rem and Rsm ratios. In Ref. [142], the angular 
distributions of 14 recoil polarization response functions and two Rosenbluth combinations have 
been measured for ep epn^ at = 1 GeV^. The analysis of Ref. [142] avoids mailing a 
multipole truncation. In particular for Rem, one sees from Fig. 64 that it yields a value which 
is lower than the one of Ref. [81] based on a multipole truncation, in agreement with the trend 
observed when applying the correction of Eq. (6.4). 

In the following, we briefly describe the partonic estimate of Ref. [320] for the 27 contribution 
to the iV ^ A electroproduction amplitudes at large Q'^. 
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Fig. 65. Handbag approximation for the eA^ eA process at large momentum transfers. In the par- 
tonic scattering process (H), the lepton scatters from quarks with momenta pq and p^. The lower blob 
represents the ^ A GPDs. 



6.2 Partonic calculation of two-photon exchange effects 



In the partonic framework, illustrated in Fig. 65, the 27 amplitude is related to the ^ A 
GPDs, as discussed in Sect. 2.7. This calculation involves a hard scattering subprocess on a 
quark, which is embedded in the proton by means of the — ^ A GPDs. Besides the domi- 
nant vector GPD Hm, corresponding with the magnetic dipole 7A^A transition considered in 
Sect. 2.7, it also involves the leading axial- vector GPD, denoted by Ci, which are the leading 
GPDs in the large Nc limit [121] contributing to this transition. Analogously to the relation of 
Eq. (2.57) relating the GPD Hm to the 7*A^A magnetic dipole form factor G^j, one can relate 
the GPD Ci to the 7*A^A axial transition form factor C^, introduced by Adler [322]. 

The eA^ — > eA handbag amplitude can then be specified in terms of the following two char- 
acteristic integrals [320] : 



J-l X 



C* 



1 dx 
-1 X 



Q2 9m "t 9 a 

S - W (27) hard 

9 A + 9m 



sgm(a;) -Ci, 



(6.8) 
(6.9) 



where all hard scattering quantities in the square brackets are given in Ref. [319]. 

In terms of the above integrals A* and C*, the reduced cross section of Eq. (6.3) including 27 
corrections, can be expressed as [320] : 



, _i, 002 I . leMjg^ 



QIQ' 



SM 



+ 



G*M V 3 



N 



2 g+g_ qi^-Mn + m^ 



(6.10) 



Furthermore, the 27 exchange corrections to Rem and Rsm, are : 
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Fig. 66. dependence of the 27 corrections to Rem (upper panel) and Rsm (lower panel) at 
W = 1.232 GeV and for different values of e. Figure from Ref. [320]. 
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M 



To provide numerical estimates for the 27 corrections, one needs a model for the N ^ A 
GPDs which appear in the integrals A* and C*. In Ref. [320], the large Nc relations discussed in 
Sect. 2.7 were used. For the GPD Hm, the model of Sect. 2.7.2 is used, whereas the axial GPD 
Ci is expressed through the isovector combination Ci{x, 0, Q'^) = ^/3[H'^ — H'^]{x, 0, Q"^) of 
the nucleon axial GPDs, see Ref. [121]. 

In Fig. 66 we show the 27 corrections to Rem and Rsm estimated using the modified Regge 
GPD model. We see that the 27 effects are mainly pronounced at small e and larger Q^. For 
Rem they are well below 1 %, whereas they yield a negative correction to Rsm by around 1 %, 
when Rsm is extracted from cttl according to Eq. (6.4). 

Another way to obtain Rsm is from a Rosenbluth-like analysis of the cross section ctq for 
the eA^ eA reaction. In Fig. 67, we show the reduced cross section from Eq. (6.10). One 
sees how a different extraction of Rsm can yield a significantly different result. Starting from 
a value of Rsm (e-g-^ -10 %) extracted from axL, adding in the 27-exchange corrections gives 
a sizable change in the slope of the Rosenbluth plot. When fitting the total result by a straight 
line in an intermediate e range, one extracts a value of Rsm around 3 percentage units lower 
than its value as extracted from ari- This is sizable, as it corresponds with a 30% correction on 
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Fig. 67. Rosenbluth plot of the reduced cross section aji of Eq. (6.10) for the eA^ eA reaction at 
= 3 GeV^. The dashed curve corresponding with Rsm = —10% is the I7 result, whereas the 
solid curve represents the result including 27 corrections for the same values of Rsm- The dotted curve 
corresponding with Rsm = —13% corresponds with a linear fit to the total result in an intermediate e 
range. Figure from Ref. [320]. 

the absolute value of Rsm- This is similar using the Rosenbluth method to extract the proton 
elastic form factor ratio G^j/G^j, where the 27-corrections were found to be large [317]. 

Note that obtaining -Rsm from the Di/Aq ratio depends on isolating the cos0^ term in the 
cross section. This term implicitly uses the decay A Nn to give us some information about 
the polarization of the A. The "Di/Aq" method is thus the analog of the polarization transfer 
measurement of Ge/Gm in the elastic case. Lilce the elastic case, the calculations of Ref. [320] 
show that there are sizeable corrections to the Rosenbluth-like determination of Rsm, and sig- 
nificantly smaller corrections if Rsm is obtained from Di/Aq- It will therefore be interesting to 
confront the Di/Aq extraction of Rsm with the new Rosenbluth separation data in the A region 
up to ~ 5 GeV^ which are presently under analysis [323,324]. 
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7 Conclusions and outlook 



In this work we have been reviewing the recent progress in understanding the nature of the 
A-resonance and of its electromagnetic excitation. 

On the experimental side, high-precision data for low-Q^ pion electroproduction in the A- 
resonance region have become available from BATES, MAMI and JLab. As a result of these 
measurements, a consistent picture has emerged for the small quadrupole ratios Rem and 
RsM- The electric quadrupole over magnetic dipole ratio Rem is around — 2 %. The Coulomb 
quadrupole over magnetic dipole ratio Rsm is around —6 % in the range up to about 1 GeV^, 
and becomes more negative with increasing values of Q^. The measurements at large show 
no trend towards the pQCD predictions. Up to the largest measured values around 6 GeV^, 
the dominant magnetic dipole 7*A'^A form factor decreases faster than the pQCD scaling pre- 
diction ~ 1/ Q S and the Rem ratio shows no trend towards a zero crossing to reach the pQCD 
prediction of +100 %. 

The small value of the Rem ratio has been understood in the large Nc limit where the Rem 
ratio is of order ~ ^/N^- The large A^'c limit leads to the spin-flavor symmetry which is used 
in many quark models. In such quark models the small non-zero values of Rem and Rsm are 
interpreted as indicators of hadron deformation, because they require d-wave components in the 
and/or A-wave functions. The constituent quark models yield however too small (absolute) 
values for the Rem and Rsm ratios, and underpredict the Ml iV — > A transition by about 25 
%. A natural way to improve on such quark models is to include pion degrees of freedom. A 
variety of models have shown that pionic effects to the 7A^A transition are indeed sizeable. In 
particular the strong tensor interaction due to one-pion exchange between quarks may yield to 
a simultaneous spin flip of two quarks in the nucleon and yield non-zero E2 and C2 amplitudes 
even if the quark wave functions have no rf- state admixture. 

The effect of pionic degrees of freedom in the 7A^A transition has been discussed within the 
context of dynamical models. We have reviewed in detail three examples of such models and 
compared their results. These dynamical models yield a good description of experimental data 
and provide a model-dependent framework for extracting the 7A'^A form factors. 

A systematic framework to account for the pion loop effects in the ^ A transition has 
emerged in recent years in terms of the chiral effective field theory. Namely, the extensions of 
chiral perturbation theory to the A(1232)-resonance energy region have been developed. The 
5-expansion is based on the small parameter 5 equal to the excitation energy of the resonance 
over the chiral symmetry-breaking scale. The other low-energy scale of the theory, the pion 
mass, counts as 5'^, which is crucial for an adequate counting of the A-resonance contributions 
in both the low-energy and the resonance energy regions. This framework has been applied to 
the process of pion electroproduction. A complete calculation of this process in the resonance 
region up to, and including, next-to-leading order in the 5-expansion satisfies gauge and chiral 
symmetries perturbatively, and Lorentz-covariance, analyticity, unitarity (Watson's theorem) 
exactly. The low-energy constants (LECs) entering at this order are the 7A^A couplings qm, 
9e, gc characterizing the Ml, E2, C2 transitions, respectively. They are extracted from a fit to 
the pion photoproduction data, and for gc to the pion electroproduction at low Q'^. Once these 
LECs are fitted, both the energy (W) and Q'^ dependencies follow as predictions. The NLO re- 
sults were found to give a good description of both the energy-dependence around A-resonance 
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and the low Q'^ dependence of the pion electroproduction data. 

We have also reviewed the recent progress achieved in the lattice QCD simulations of the 
77VA transition. They now clearly establish a non-zero negative value for the Rem and Rsm 
ratios. However the present lattice calculations have been obtained for pion masses well above 
300 MeV. A naive linear extrapolation in the quark mass to its physical value, yields a value for 
Rsm which is however in disagreement with experiment. It has been reviewed that the xEFT 
framework plays a dual role in that it allows for an extraction of resonance parameters from ob- 
servables and predicts their pion-mass dependence. It was seen that the opening of the A — ^ ttA^ 
decay channel at = Ma — Mjv induces a pronounced non-analytic behavior of the Rem 
and Rsm ratios. While the linearly-extrapolated lattice QCD results for Rsm are in disagree- 
ment with experimental data, the xEFT prediction of the non-analytic dependencies suggests 
that these results are in fact consistent with experiment. 

The quark structure of the electromagnetic N ^ A transition can be accessed in hard scat- 
tering processes such as the deeply virtual Compton scattering processes 7*A^ — > 7A, where 
the virtual photon 7* carries a large virtuality. The non-perturbative information in this pro- 
cess can be parametrized in terms of — > A generalized parton distributions. The ^ A 
transition form factors are obtained as first moments of these GPDs. In the large Nc limit, the 
dominant magnetic GPD is related to an isovector combination of nucleon GPDs. Using this 
large Nc relation combined with a Regge parametrization for the nucleon GPDs, fitted to the 
elastic nucleon form factors, the magnetic dipole GPD was calculated. Its first moment yields a 
prediction for the dependence for the Ml ^*NA form factor. This prediction is very well in 
agreement with the recent data up to the highest measured value of around 6 GeV^. This cal- 
culation also shows that the Ml 7*A^A form factor follows approximately the same fall-off as 
the nucleon (isovector) Pauli form factor, yielding a faster fall-off at large Q"^ compared with the 
pQCD prediction. Likewise, we have used a large Nc relation between the low Rem 
ratio and the neutron electric form factor to calculate the electric quadrupole N A GPD. The 
resulting prediction for Rem is in good agreement with the data up to around 1.5 GeV^. 

At large momentum transfer, we discussed the 27 contribution to the eA^ — * eA(1232) 
evriV process. The 27 amplitude has been related in a partonic model to the iV — > A GPDs. For 
Rem, the 27 corrections were found to be small, below the 1 % level. However, the neglect of 
a quadratic term in a truncated multipole expansion to extract Rem yields corrections at 
the 1-2 % level. The 27 corrections to Rsm on the other hand are substantially different when 
extracting this quantity from an interference cross section or from unpolarized cross sections, as 
has been observed before for the elastic eA^ eN process. They turn out to be important when 
extracting Rsm at large from the unpolarized longitudinal cross section using a Rosenbluth 
separation. 

We end this review by spelling out a few open problems and challenges (both theoretical and 
experimental) in this field : 

(1) Issues in xEFT: 

a) To interpret the precise data on pion photo- and electroproduction at low Q'^ with mini- 
mal model assumptions, requires xEFT calculations to go beyond making predictions of 
7*A^A form factors and directly perform the comparison at the level of the observables. 
The xEFT to NLO in the (5-expansion power counting scheme have shown to yield an 
overall satisfying description of a large number of pion photo- and electroproduction 
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observables. 

b) The xEFT calculations in the (5-expansion have been obtained at NLO for the 7*iV — 
ttN process. The present theoretical uncertainties of these calculations are larger than the 
precision of present data. The next step for the theory is to go beyond a NLO calculation. 

c) The xEFT calculation also allows to make predictions for Compton and virtual Compton 
scattering processes in the A region. A comparison of the NLO predictions for 7*A^ 
TiN , real and virtual Compton scattering may yield a powerful cross-check of the theory. 

d) In the low Q"^ region, it will be very worthwhile to measure the complementary channels 
in particular virtual Compton scattering in the A region. 

e) The Ml excitation of the A gives a large (para)magnetic contribution to the magnetic 
polarizability of the nucleon {/^ Q ■ 10~^ fm^ ). The challenge is to reconcile this large 
contribution with the small value observed in experiment. 

f) At present, xEFT calculations are performed in different power counting schemes. A 
challenge it to combine the low-energy chiral expansions around threshold with the 5- 
expansion around A-resonance. 

(2) Resonances on the lattice: 

Lattice QCD calculations with dynamical quarks have just started for the electromag- 
netic — > A form factors. They are currently performed for pion masses above 300 
MeV where the A is stable. The xEFT shows a strong non-analytic behavior in the 
quark mass for the 7*A^A form factors when the pion mass becomes smaller than the 
N — A mass difference. It will be a challenge to implement the unstable nature of the 
A on the lattice and to extend these calculations to pion mass values sizeably below 300 
MeV. 

(3) Development of dynamical models: 

Building in symmetries in a consistent way is a non-trivial challenge for the dynamical 
models. Future challenges in this field are to find ways to systematically improve on 
approximations, identify the limitations, and quantify the theoretical uncertainties. 

(4) TheN ^ A GPDs: 

a) The parametrization of the iV — * A GPD for the C2 transition has not been addressed 
until now. Model calculations of — > A GPDs as well as lattice simulations of its 
moments will be very useful in order to make a comparison with the quark distributions 
in a nucleon, as parametrized through the nucleon GPDs. 

b) Of high interest is the measurement of the deeply virtual Compton scattering 7*A^ — > 7A 
process to directly access the N ^ A GPDs. 

(5) Large momentum transfer regime: 

a) At presently accessible momentum transfers (Q^ < 6 GeV^), the pQCD prediction 
(Rem — 1, RsM const) seems to be nowhere in sight. A challenge is to map out 
the transition from xEFT at low Q"^ to pQCD by measuring the 7*A^A form factors to 
larger Q'^ (e.g. in the near future at the planned JLab 12 GeV facility). 

b) At larger Q^, constrain the model dependence in the extraction of multipoles. In this 
respect, double polarization observables have already been shown to be a very valuable 



131 



cross-check and should be further investigated, 
c) Check the importance of 27 exchange corrections by Rosenbluth measurements (via the 
gl cross section) of Rsm- 



In retrospect, the impressive level of precision obtained in the measurements of the electro- 
magnetic excitation of the A-resonance during the past few years has challenged the theory to 
arrive at a QCD based description of the A-resonance. Several new theoretical developments 
are under way and have shown promise and potential. We hope that the present work will stim- 
ulate further efforts in this direction and contribute to the very fruitful interplay between theory 
and experiment in this field. 
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A General analysis of the electromagnetic pion production on the nucleon 



In this Appendix we present the general properties of the pion photo- and electroproduction 
reactions: 



7{q) + N{p)^N{p')+7i{k), 

e{k) + N{p)^e'{lf) + N{p')+7r{k) 



(A.l) 
(A.2) 



Namely, we consider the kinematics, various decompositions into invariant amplitudes, the mul- 
tipole expansion, and the isospin structure. 



A.1 Kinematical details 



In both cases, the invariant T-matrix can be expressed as 



(A.3) 



where is the nucleon electromagnetic current matrix element 



J/* p'k\f\p >, 



(A.4) 



and is the photon polarization vector for photoproduction and 



_ eu{qf)^^u{qi) 

for electroproduction, with Q'^ = —q^ = —{If — li^. 
The T-matrix can be expanded in terms of Lorentz-covariant operators as 



(A.5) 



u{p') 



u=i 



u{p), 



(A.6) 



where s = (p + g)^, t = [p' ~ p)^ are the Mandelstam variables, n = 4(6) for photoproduction 
(electroproduction), and A's are scalar amplitudes. The expressions for the operators Mj can be 
found in Refs. [163,164,204]. 
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A. 2 CGLN decomposition 



Traditionally the Lorentz structure of the photopion amplitude is expressed in terms of the 
Pauli matrices and (two-component) spinors [163]: 

Tf,{-f*N ^ TTiV) = — -4 ^Xi, (A.7) 
M]si 

with JT given by [10,188]: 

— ia ■ eFi-\- (T ■\^a ■ (q x e) F2 + i<7 • q k • e F3 + i<T • k k • e F4 
+^<T • qq • e F5 + i<T • k q ■ e Fg — zcr • keo — icr • qeo Fg, (A. 8) 

where cr is spin operator of the nucleon and k = k/|k|, q = q/|q|. Note that with the use of 
cr = cr — (cr • q)q and k = k — (k • q)q in the first term of Eq. (A. 8), the F5 and Fq terms differ 
from the ones often used in the literature, e.g., [204,165]. 

Obviously, the first four terms in Eq. (A. 8) arise due the transverse polarizations of the pho- 
ton, while amplitudes F5 and Fq come from the longitudinal polarization of the virtual photon. 
The longitudinal component of the current is related to the scalar component via current con- 
servation and hence: 

\<l\F, = qoFs, \<l\F^ = qoF^. (A.9) 



Note also that the scalar amplitudes in Eq. (A. 8) (commonly knows the 'CGLN amplitudes') 
are complex functions of three independent kinematical variables, e.g., Fi = 0^, Q^), 

with the total energy W, the pion scattering angle 6^, and the photon virtuality Q'^. 



A. 3 Helicity amplitudes 

Cross sections and polarization observables of pion photo- and electroproduction are also 
often expressed in terms of the helicity amplitudes (see, e.g., [208]): 

ff^'M =< /^'l-^l V >, (A.10) 

where ji {jj!) denotes the helicity of the intitial (final) nucleon and A the photon helicity. 

There are 6 (4) independent helicity amplitudes for the case pion electro- (photo-) production. 
They are defined Table A.3 provided ij! — \ and the phase convention of Jacob and Wick [325] 
is used. 

The six helicity amplitudes Hi,...,H^ can be expressed as the following linear combination 
of the six independent CGLN amplitudes Fi, Fg: 
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Table A. 1 

HeUcity amplitudes if^i^i xix- 

Hi^-^sinO cos(^/2) (F3 + F4) , 

i/2 = A/2cos(^/2) [(F2 - Fi) + i(l - cos^) (F3 - F4)] , 

= sin^ sin(^/2) (F3 - F4) , (A. 11) 

i/4 = V2sin(^/2) [(Fi + F2) + i(l + cos^) (F3 + F4)] , 

He = sm{e/2) (Fe-F,). 

A.4 Multipole decomposition 

Multipole amplitudes are obtained with the use of eigenstates of parity and angular momen- 
tum, instead of plane wave states. For electroproduction, there are six types of transitions pos- 
sible to a ttN final state with angular momentum / and parity P, which are classified according 
to the character of the photon, transverse or scalar (or alternatively, longitudinal) and the total 
angular momentum J = Z ± | of the final state. In addition, the transverse photon states can 
either be electric, with P = (—1)^, or magnetic, with P — (—1)^+^, where L is the total orbital 
angular momentum of the photon. 

The expansion of the CGLN amplitudes in terms of the multipole amplitudes is given as 
follows [204]: 

Fi = J2{ilMi+ + Ei+)PU, + [{I + 1)M,_ + Ei_]PU}, 

l>0 

F2 = J2[{l + l)Mi+ + lMi_]Pl, 

i>i 

F3 = J2[iEi+ - Mi+)Pl'^, + {Ei_ + Mi_)Pi,l 
i>i 

F4 = - Ei+ - Mi_ - Ei_)Pl', (A.12) 

l>2 

F, = J2[{l + l)Li^P;^,-lLi_Pi,], 

Fe = Y.[lLi--{l + l)Li+]Pl 
i>i 

where P/s are the derivatives of the Legendre polynomials of argument x — cosO. Note that in 
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the literature the scalar transitions are sometimes employed and described by Si± multipoles, 
which correspond to the multipole decomposition of the ampUtudes F7 and Fg. They are con- 
nected with the longitudinal ones by Li± = {ojq/l q |)'S'/±. 



A.5 Isospin decomposition 



Assuming isospin conservation, there are three independent isospin amplitudes for the single 
pion production off a nucleon by a photon. These are: for the isoscalar photon, and A^^^"^^ 
and A*^'^/^) for the isovector photon. The latter two correspond to the nN system with total 
isospin 1=1/2 and / = 3/2, respectively. The isospin structure of any photopion amplitude 
{e.g., the multipoles) can be written as: 



(A.13) 



where a is the pion isospin index and Tq, the nucleon isospin matrices. The isospin amplitudes 
and are related to the amplitudes of Eq. (A.13) as : 



^(3/2) ^ ^(+) _ ^C- 



^(1/2) ^ ^{+) + 2A(- 



(A.14) 



It will also be useful to define the proton pA'^^Z'^'' and neutron „^^^/^^ amplitudes with total 
isospin 1/2, 



/I (1/2) _ /l{0) , 1/1(1/2) .(1/2) _ 4(0) _ 1/1(1/2) 

pj± — J± -t- ^yi , nJi — J± ^Ji 



(A.15) 



With this convention the physical amplitudes for the four physical pion photo- and electropro- 
duction processes are 



A(7*p 7171+) = a/2 



^A(V2) _ ^^(3/2) 



yl(7*p^p7r°)=p/l(^/2)^f^(3/2)^ 

3 

A(7*n p7r~) = a/2 



4(1/2) , 1^(3/2) 
3 



(A.16) 



Under the interchange of s and u (crossing), the functions A\ ' are either even or odd, i.e., 
and this crossing symmetry can readily be derived from the assumption of charge conjugation C 
invariance for photo- and electroproduction. These crossing symmetry properties can be com- 
pactly expressed in matrix form with a six-vector A with elements Ai,.., Aq, 
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A{s,t,u,Q^) = [^]A{u,t,s,Q^), (A.17) 
where [^] is a 6 x 6 diagonal matrix 

[el = e diag{l, 1, -1, 1, -1, -1}, (A.18) 

and the parameter ^ is defined by ^ = 1 for isospin index (+, 0), and ^ = — 1 for isospin index 
(-)• 
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